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Thermal  quench  residual  stresses "" 
are  typically  compressive  at  the 
surface  and  tension  in  the  interior. 
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Residual  stress  distribution  of  A1 7075 
cylinder  quenched  in  cold  water  spray 
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Quench  rate  can  appreciably  impact 
residual  stress  magnitude  and  range. 
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Stress  range:  32  ksi  20  ksi  3  ksi 

Faster  quench 

.  Effect  of  quench  water  temperature  on  7XXX  A1  forging  residual  stress 


Advanced  stress  relief  technology  is  now 
commercialized  to  alleviate  forging  machining 
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Check  point  number 

Contrast  in  machining  performance  of  aluminum  7050-T74xx  die-forging 
Alcoa  quench,  water  quench,  compressive  stress  relief  in  finish  dies 
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Illustration  of  mechanical  stretch 
associated  with  cold  restrike 


Residual  stress  bias  can  be 
masked  as  a  geometry  effect. 
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Crack  initiates  where 
stress  normal  to  crack 
plane  is  tensile 
(high  FCGR) 


Crack  initiates  where 
stress  normal  to  crack 
plane  is  compressive 
(slow  FCGR) 


Effect  of  specimen  type  on  FCGR  measurement 
from  non-stress  relieved  7xxx  Al  extruded  rod. 
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Fracture  toughness  test  results  from 
partially  stress  relieved  product  can  be 
misleading  without  proper  interpretation. 
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Specimen  width,  W  (in.) 

nteraction  of  internal  stress  state  and  specimen  size  on  fracture 
toughness  measurements  from  comparable  4-in.  thick  7050-T74 
hand  forged  billets,  one  stress  relieved  and  the  other  not. 
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Effect  of  through-thickness  internal  stress  distribution  (normal  to  the 
notch  plane)  on  fracture  toughness  specimen  precrack  shape. 


Residual  stress  induced-  J 
bending  can  measurably  impact 
crack  tip  stress  intensity  factor. 


ALCOA 


(a)  Specimen  location 
within  parent  slab. 


(b)  isolated  specimen 
longitudinal  residual 
stress  distribution. 


(c)  Clamping  moment 
developed  after 
machining  crack 
starter  slot. 


Portrayal  of  crack  tip  clamping  associated  with 
residual  stress  distribution  aligned  paraiiel  to  a 
compact  specimen  notch  plane. 
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Classical  LEFM  similitude  breaks  down 
when  residual  stress  is  appreciable  &  ignored. 

Part  Coupon 
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Residual  stress  induced  bias  is  magnified  in  the  coupon  test 


r. 


A  simplified  correction  practice 
has  been  devised  for  Kic  testing. 


ALCOA 


A8  =  5l-82 

Measure  the  specimen  height  before  and 
after  machining  the  crack  starter  notch. 


Load 
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5%  secant  offset 
line  from  new  origin 


Analyze  test  record  with 
new  origin  displaced  by  A6 


Crack  opening  displacement 

Kic  test  residual  stress  correction  schematic 
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New  quench  technology  and  better 
test  methods  are  available  to 
reduce  the  residual  stress  effect. 
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Quench  &  age,  **  Quench,  compression  stress  relieve  &  age 


protocol  to  purge  residual  stress  bias 
)m  FCGR  data  has  been  established. 

Sued,  ASTM  STP  743,  1981  &  Bush  et.  al,  ASTM  STP  1189,  1993) 
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FCGR  data  from  two  incompletely  stress 
relieved  7050-T7452  forgings 

(S-L  orientation,  R  =  0.33,  high  humidity  air) 


Closure-based  correction  of 
suspect  FCGR  data  is  derivable 
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Crack  opening  displacement  (v) 

Effect  of  residual  stress  on  load-COD  trace 


Residual  stress  bias  can  be 
deduced  from  decay  in  closure 
load  with  crack  advance. 


ALCDA 
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Load-COD  traces  from  FCG  tests  of  high  strength  Al 
product  evaluated  in  a  non-stress  relieved  (top)  and 
_ stress  relieved  (bottom)  condition. 
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Alloy  7050  plate  and  forging  display 
similar  FCGR  behaviors  when  the 
residual  stress  effect  is  purged. 
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Closure  corrected  FCGR  data  for  7050-T745X 
plate  and  non-stress  relieved  forgings 

(S-L  orientation,  R  =  0.33,  high  humidity  air) 


Suspect  FCGR  data  on  varied  products 
is  readily  found  in  the  open  literature. 
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AK(ksiVin)  AK,  ksi  VTrT 

Grade  A/B  D357-T6  Casting  5xxx  A1  Alloy  Butt  Weld  Joint 


Conclusions/recommendations 

Understanding  of  heat  treat  induced  residual  stress  cause/effects 
on  toughness  &  crack  growth  characterizations  has  been 
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Absuact 

A  study  was  performed  using  two  aluminum  alloys,  2024-T351  and  7075- 
T651,  subjected  to  prior  corrosion  and  then  fatigue.  The  resulting  crack  nucleation 
and  the  "short"  fatigue  crack  growth  were  observed  using  an  "in-situ"  scanning 
electron  microscope.  This  unique  machine  allowed  visual  observation  of  a  crack  on 
the  surface  of  the  material  during  fatigue  cycling  at  high  magnifications.  Baseline 
tests  without  corrosion  were  conducted  to  compare  the  "short"  fatigue  crack  growth 
rates  with  those  of  prior  corroded  specimens. 

Four  tests  were  conducted.  One  specimen  was  subjected  to  corrosion  prior 
to  fatigue  and  one  specimen  was  run  as  a  basehne  fatigue  test  for  each  of  the  two 
aluminum  alloys.  Specimens  were  pitted  in  3.5%  salt  water  for  28  hours.  It  was 
hypothesized  that  the  prior  corrosion  would  accelerate  the  fatigue  crack  growth  rate 
in  the  "short"  crack  region  and  that  cracks  would  form  from  discontinuities  on  the 
surface,  such  as  corrosion  pits. 

Based  on  the  results,  it  was  found  that  prior  corrosion  does  influence  the 
fatigue  crack  growth  rate  and  that  the  two  materials  exhibit  different  crack  growth 
behavior.  The  specimens  subjected  to  prior  corrosion  exhibited  accelerated  crack 
growth  rates  compared  to  the  baseline  tests. 
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The  7075  specimens  had  faster  crack  growth  rates  than  the  2024  specimens. 
The  basehne  tests  and  the  prior  corroded  2024  specimen  displayed  crack  growth 

rates  that  increased,  then  decreased,  then  increased  again  as  AK  increased.  The 

7075  specimens  did  not  show  a  dramatic  drop-off  in  crack  growth  rate  as  seen  in 
the  2024  tests. 

This  work  is  unique  in  that  combining  prior  corrosion  and  short  crack 
behavior  has  not  been  greatly  researched.  These  are  relatively  new  fields  of 
research  but  they  are  important  because  they  model  real-world  behavior  of 
materials.  Further  work  should  be  conducted  in  this  area  to  statistically  substantiate 
the  results  presented  here. 


Introduction 

Research  at  BatteUe  has  found  that  corrosion  costs  the  United  States  300 
billion  dollars  annually,  one  third  of  which  is  avoidable  [1].  Corrosion  acting  with 
fatigue  is  a  major  maintenance  concern  in  the  aircraft  industry.  As  the  aircraft  fleet 
ages,  these  time  dependent  failure  modes  are  becoming  even  greater  concerns.  In  a 
recent  survey  performed  by  the  authors  and  their  colleagues,  corrosion  and/or 
fretting  were  found  to  be  a  contributing  factor  in  at  least  687  incidents  and  accidents 
from  1974  to  1994  in  military  and  commercial  aircraft  in  the  United  States  [2]. 

Often  aircraft  are  left  to  sit  while  not  in  service  at  which  time  corrosion  may 
occur.  Because  inspections  are  currently  based  on  flight  hours  and  not  calendar 
hours,  the  corrosion  may  not  be  discovered  until  the  aircraft  has  been  put  into 
service.  During  this  time,  however,  cracks  may  have  nucleated  from  the  corrosion 
damage,  which  can  grow  during  subsequent  flight  of  the  aircraft.  These  cracks  are 
cause  for  alarm  as  they  may  not  be  discovered  before  they  have  become  critical. 
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Because  of  this,  prior  corrosion  damage,  such  as  pitting  damage,  is  cause  for 
concern  in  aviation  safety. 

"Short"  fatigue  crack  growth  added  to  prior  corrosion  damage  may  lead  to 
reduced  component  lifetime.  It  has  been  shown  in  the  past  two  decades  that  "short" 
crack  growth  rates  vary  from  long  crack  growth  rates  in  that  "short"  cracks  grow 
faster  below  the  threshold  stress  intensity  value.  Therefore,  life  predictions 
according  to  long  crack  growth  data  may  be  unconservative  when  a  crack  is  in  the 
"short"  crack  region.  Because  of  this  the  effort  in  "short"  crack  research  has  been 
increased  but  Utde  has  been  done  concerning  prior  corrosion  damage  in  the  "short" 
crack  regime.  Coupled  together,  prior  corrosion  and  "short"  crack  growth  may 
reduce  the  expected  lifetime  or  time/cycles  between  inspection  intervals  of  an 
aircraft  significantly. 

Prior  Corrosion 

Much  research  has  been  done  on  simultaneous  corrosion  and  fatigue 
damage  and  it  has  been  found  that  corrosion  accelerates  the  rate  of  fatigue  and 
failure.  However,  often  the  corrosion  and  fatigue  act  sequentially.  When  an 
aircraft  is  in  flight  it  experiences  cyclic  loading.  When  it  is  on  the  ground,  moisture 
builds  up  on  the  aircraft.  As  an  aircraft  sits  (and  when  it  is  in  the  air),  the  corrosion 
process  takes  place.  Thus,  the  sequential  action  of  these  time  dependent  failure 
modes  may  produce  different  results  than  the  simultaneous  action. 

It  has  been  found  that  pit  characteristics  vary  according  to  loading, 
environment,  and  material  [3].  When  comparing  the  sizes  and  shapes  of  pits 
formed  in  7075-T6  from  zero  load,  sustained  load,  and  fatigue  load,  Lorie  Grimes 
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found  that  the  corrosion  pits  formed  by  the  fatigue  load  were  larger  than  the  others 
and  that  they  were  more  spread  across  the  surface.  Therefore,  it  can  be  seen  that 
the  expected  fatigue  lives  due  to  each  will  vary  from  each  other. 

Du,  Chiang,  Kagwade,  and  Clayton  found  that  when  aluminum  alloy  2024- 
T3  was  fatigued  and  then  exposed  to  a  corrosive  solution  and  then  fatigued  again 
that  the  fatigue  life  was  greater  than  had  the  material  been  corroded  and  then 
fatigued  [4]. 

Short  Cracks 

The  idea  that  the  behavior  of  “short”  cracks  varies  from  long  cracks  did  not 
occur  until  the  mid  1970s.  All  work  on  crack  growth  was  done  assuming  long 
crack  behavior.  Scientists  discovered  that  in  the  short  crack  region  a  crack  exhibits 
characteristics  which  are  not  valid  using  linear  elastic  fracture  mechanics  (LEFM). 
Over  the  past  decade  it  has  been  realized  that  short  crack  growth  may  change  the 
estimated  life  of  a  component.  Short  cracks  may  make  the  estimated  life  non¬ 
conservative.  Therefore,  it  is  necessary  to  analyze  the  behavior  of  short  cracks  in 
order  to  learn  their  effects  on  crack  propagation  rates  and  how  to  take  these  potential 
effects  into  account. 

To  distinguish  between  a  long  and  short  crack  would  perhaps  simply  be  to 
say  where  LEFM  is  valid  and  where  it  is  not,  respectively.  R.O.  Ritchie  and  S. 
Suresh  have  discussed  the  types  of  short  cracks  and  what  mechanisms  lead  to 
accelerated  short  crack  growth  [5].  They  relate  crack  driving  force,  local  plasticity, 
microstructure,  crack  shape,  crack  extension,  premature  crack  closure,  and  local 
crack  tip  environment  to  the  differences  between  long  and  short  crack  growth.  A 


604 


crack  may  be  short  or  small  on  several  bases  but  there  is  no  quantitative  value  that  is 
yet  agreed  upon  to  say  when  a  crack  is  short  and  when  it  is  long.  According  to 
Ritchie  and  Suresh,  a  crack  may  be  short  with  respect  to  the  microstructure,  the 
section  size  (physically  short),  or  the  scale  of  local  plasticity  [6].  Physically  short 
cracks  were  observed  in  this  study  where  the  difference  in  minimum  and  maximum 
stress  intensity  factors  below  the  threshold  value. 

The  Combination  of  Aluminum.  Prior  Corrosion,  and  Short  Crack  Growth 

A  crack  spends  much  of  its  life  in  the  nucleation  and  formation  stages  where 
crack  growth  begins  in  the  short  crack  region.  To  understand  crack  behavior,  short 
crack  studies  need  to  be  performed  to  prevent  unexpected  failure.  It  is  obvious  that 
much  more  study  needs  to  be  done  involving  the  synergistic  effects  of  corrosion 
and  crack  formation  and  propagation  on  common  aircraft  materials,  such  as  2024- 
T3  and  7075-T6  aluminum  alloys. 

Test  Methodology 


Test  Apparatus 

The  fatigue  tests  were  conducted  in  an  electro-hydrauhc  servo-controUed  25 
kN  “in-situ”  fatigue  machine  coupled  with  MTS  components.  This  machine  is 
unique  in  that  the  specimen  is  attached  in  the  grips  inside  a  scanning  electron 
microscope  (SEM).  Thus,  the  specimen  may  be  monitored  at  high  magnifications 
while  the  fatigue  load  is  apphed  and  cracks  may  be  identified  as  they  form  and 
grow. 
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specimen  Description 

Because  of  their  widespread  use  in  aircraft,  which  are  susceptible  to  both 
fatigue  loading  and  corrosion  pitting,  2024-T3  and  7075-T6  aluminum  alloys  were 
studied.  The  material  was  obtained  in  rod  form  (0.500  in  OD)  from  Affiliated 
Metals.  A  groove  was  machined  into  the  rectangular  cross-section  of  the  specimens 
to  act  as  a  point  of  stress  concentration.  This  was  done  to  ensure  crack  growth  on 
the  side  of  the  specimen  facing  the  SEM  column.  It  was  found  that  the  stress 
concentration  factor  in  the  groove  was  1.409  according  to  [7].  Two  baseline 
specimens,  one  of  each  material,  were  tested  without  the  prior  corrosion.  These 
tests  were  compared  with  the  two  specimens  run  with  the  prior  corrosion  to 
determine  if  there  was  a  change  in  short  crack  growth  rate  due  to  the  prior 
corrosion. 

Pits  were  introduced  to  the  surface  prior  to  the  fatigue  loading.  This  was 
done  to  model  prior  corrosion  damage  on  structures.  The  specimens  were  covered 
with  paraffin  wax  all  over  to  protect  the  metal  except  for  the  small  region  of  the 
machined  groove  where  the  pitting  would  take  place.  The  specimens  were  pitted  in 
a  2  liter  bath  of  3.5%  NaCl  and  distilled  water  at  25  °C  for  28  hours.  The 
specimens  were  rinsed  with  distilled  water  and  acetone  and  tested  immediately  after 
being  corroded. 

Test  Conditions 

Fatigue  tests  were  conducted  at  a  frequency  of  10  Hz.  The  R  value 
(Pmii/PmJ  was  0.3.  Using  a  spreadsheet  and  an  initial  crack  length  of  0.01 
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mm,  the  load  was  calculated  iteratively  to  obtain  a  value  of  AK  around  1.0.  The 


solution  of  AK,  the  change  in  stress  intensity  factor,  was  derived  for  a  part-circular 
crack  in  a  simple  stress  field  on  the  surface  of  the  material.  From  A.C.  Pickard’s 
solution  for  a  part-circular  crack,  for  this  experiment  AK  becomes: 

A  K  =  MqMbMj  —A(j-\fm 
n 

where:  K  =  stress  intensity, 

Mq  =  general  correction  factor. 

Mg  =  back-face  correction  factor. 

Mg  =  side-face  correction  factor, 

Aa  =  change  in  stress, 
c  =  1/2  crack  length. 


The  experimental  matrix  is  shown  in  the  following  table. 


Table  1.  Experimental  Matrix  and  Testing  Loads 


Specimen 

Pre-crack  Load, 
kN  (%S,J 

Test  Load,  kN 
(%SJ 

Xl-2024 

2.90  (85%) 

2.90  (85%) 

P4-2024 

2.39  (70%) 

2.39  (70%) 

P2-7075 

3.84  (75%) 

3.07  (60%) 

P3-7075 

3.07  (60%) 

3.07  (60%) 

A.K.  Vasudevan  and  S.  Suresh  found  the  value  of  AKjjj  for  alloys  2024-T3  and 
7075-T6  exposed  to  a  corrosive  environment  [85].  For  2024-T3,  AK^j,  =  3.2 


MPaVm  and  for  7075-T6,  AK^  =  2.1  MPa'^m. 
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Crack  Growth  Readings 


Crack  growth  was  measured  using  photographs  taken  with  the  SEM.  A 
micron  marker  is  given  on  a  photograph.  This  allows  the  length  of  a  crack  to  be 
found  by  measuring  the  crack  at  a  known  magnification  and  multiplying  it  by  the 
micron  marker  value.  For  example,  at  lOOOx  magnification,  1  cm  equals  10  mm. 

A  crack  that  measured  3  cm  on  a  photograph  would  actually  be  30  pm  long. 

To  find  a  crack  initially,  the  surface  was  monitored  every  5,000  cycles. 
Once  a  crack  had  formed  and  was  detected,  measurements  were  taken  when  crack 
growth  was  visible  on  the  SEM.  The  specimen  was  loaded  as  crack  readings  were 
taken  so  that  the  crack  would  be  visible.  Crack  readings  were  taken  until  a  crack 
length  of  1-2  mm  was  measured  at  which  point  long  crack  growth  behavior  would 
take  over. 


Discussion  of  Results 


Baseline  Test  Results 

The  first  test  was  specimen  XI -2024.  Many  crack-like  features  were  visible 
on  the  surface  of  this  specimen  but  were  aligned  in  the  direction  of  loading  and  after 

numerous  cycles  did  not  appear  to  be  growing.  The  first  crack  detected  was  48  pm 

long.  Six  more  cracks  were  found  and  data  were  also  taken  on  those. 

The  second  test  was  P2-7075.  The  surface  of  this  specimen  contained  more 
discontinuities  than  the  2024  specimen.  It  had  a  cleaner  surface  in  that  it  was  freer 
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of  polishing  residue.  The  first  crack  was  detected  at  a  length  of  27  pm.  Eight  more 

cracks  were  discovered  and  all  but  one  propagated. 

For  each  of  the  specimens,  the  multiple  cracks  formed  at  constituent 
particles  (inclusions  and  intrinsic  discontinuities)  on  the  surface  that  appeared  to 
have  been  filled  with  polishing  compound  that  did  not  rinse  off  in  the  acetone  bath. 
The  cracks  grew  from  one  constituent  particle  to  another.  The  line  of  the  crack 
growth  was  not  always  nomal  to  the  direction  of  loading. 

These  two  specimens  displayed  similar  fatigue  crack  growth  data.  The 
crack  growth  rate  increased,  decreased,  and  then  tended  to  increase  again.  This  is 
common  behavior  for  short  crack  data.  Crack  #1  on  specimen  XI -2024  displays 

typical  crack  growth  data  for  this  specimen,  shown  in  the  da/dN  versus  AK  plot  in 

Figure  1.  The  initial  AK  was  2.64  MPaVm.  The  data  are  bunched  around  a  crack 
growth  rate  of  l.OOE-13.  It  sharply  decreases  and  then  sharply  increases. 

Figure  2  shows  a  low  initial  AK  of  2.01  MPaVm  of  crack  #2  on  specimen 

P2-7075.  The  decrease  in  crack  growth  rate  is  not  as  dramatic  on  the  7075 
specimen  as  on  the  2024  specimen.  But  P2-7075  does  display  the  same  increase, 
decrease,  and  slight  increase  again. 

The  crack  growth  rates  of  the  two  baseline  tests  were  similar.  It  appears  as 
though  the  7075  specimen  had  a  marginally  increased  crack  growth  rate  compared 

to  the  2024.  The  values  of  AK  were  lower  on  P2-7075  as  would  be  expected 
because  the  AK^  value  of  7075  is  lower  than  2024. 
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Figure  1 .  da/dN  versus  AK  for  XI -202' 


Stress  lotcDsity  Range,  aK  (MPavm) 

Figure  2.  da/dN  versus  AK  for  P2-707; 
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Prior  Corroded  Test  Results 

The  prior  corroded  surfaces  of  the  two  aluminum  alloys  showed  damage  in 
the  form  of  pits  with  some  exfoliation  of  the  surface  before  testing  began.  Figure  3 
shows  the  damage  on  the  surface  of  specimen  P3-7075.  Two  specimens,  one  of 
each  material,  were  tested  with  the  prior  corrosion  damage. 

A  crack  of  62  pm  was  first  detected  on  specimen  P4-2024.  Four  additional 

cracks  were  monitored  on  the  surface.  This  specimen  was  highly  damaged  due  to 
corrosion,  more  so  than  the  7075  specimen.  Multiple  pits  were  bunched  together 
and  large  areas  of  exfoliation  were  present. 

The  last  specimen,  P3-7075,  was  tested  at  loads  that  were  too  high  to  keep 

the  AK  low  enough  to  remain  in  the  short  crack  region.  The  first  crack  detected  on 
this  specimen  was  160  pm  long  at  a  AK  of  7.40  MPaVm.  A  shorter  crack  was 
found  that  was  64  pm  and  a  AK  of  4.61  MPaVm.  Five  cracks  were  monitored  on 
this  specimen. 

The  surfaces  of  the  prior  corroded  specimens  were  more  difficult  to  scan 
because  of  the  corrosion  induced  damage  on  them.  While  scanning  the  surface  it 
was  difficult  to  determine  exactly  what  was  a  crack.  Mudcracking  was  visible  on 
both  specimens.  This  made  it  difficult  to  know  if  a  crack  was  growing  through  the 
material  or  whether  it  was  a  layer  of  material  peeled  from  the  surface  that  was 
cracked. 

Specimen  P4-2024  displayed  similar  crack  growth  behavior  as  the  two 
baseline  tests  but  with  an  increased  crack  growth  rate.  The  increase,  decrease,  and 
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Figure  3.  Damage  on  the  surface  of  specimen  P3-7075 
before  fatigue  load  was  applied. 
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then  increase  was  present  in  the  crack  growth  data.  The  da/dN  versus  AK  plot  for 

crack  #1  is  shown  in  Figure  4.  The  initial  AK  was  2.48  MPaVm. 

The  data  for  crack  #1  on  specimen  P3-7075  are  given  in  Figure  5.  This  plot 
shows  how  the  crack  growth  rate  for  the  cracks  on  this  specimen  increased  but  did 
not  decrease  as  seen  in  the  other  tests.  The  crack  growth  rates  for  this  specimen 
were  faster  than  the  baseline  tests  and  the  prior  corroded  P4-2024  test.  This  may  be 

caused  by  the  increased  initial  value  of  AK. 

Overall,  the  7075  specimens  had  a  faster  short  crack  growth  rate  than  the 
2024  specimens.  The  prior  corrosion  did  increase  the  fatigue  crack  growth  rate  in 
comparison  to  the  baseline  tests.  The  damage  from  the  prior  corrosion  provided 
areas  of  higher  stress  concentration  and  reduced  area  from  which  cracks  formed.  It 
was  more  difficult  to  examine  the  surfaces  of  these  specimens  due  to  the  damage  as 
would  be  the  case  in  the  field  when  inspecting  components  subjected  to  corrosive 
conditions. 

Fractographv 

Visual  examination  of  the  fracture  surfaces  showed  that  crack  growth  took 
place  only  on  the  top  face  of  the  specimen  as  was  desired.  The  shape  of  the  fatigue 
region  of  the  crack  was  not  truly  half-penny  shaped.  This  may  have  an  effect  on 
the  validity  of  the  stress  intensity  factor  solution. 

A  fractographic  analysis  verified  that  crack  nucleation  occurred  at 
constituent  particles  on  the  base  specimens  and  from  pits  on  the  prior  corroded 
specimens.  This  is  most  likely  due  to  the  increased  stress  concentration  factor  at 
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these  points  and,  on  the  pitted  specimens,  the  reduction  in  area  may  also  be  a  factor. 
Figure  6  shows  a  crack  nucleation  site  on  specimen  P4-2024.  The  radial  marks 
point  back  to  the  point  of  origin  of  the  crack.  Other  cracks  on  the  surface  are  also 
visible  on  the  surface  which  did  not  cause  failure. 

Figure  7  shows  striation-like  markings  on  specimen  XI -2024.  Striation 
marks  were  difficult  to  find  on  these  specimens  because  of  the  small  fatigue  region. 
However,  they  were  present  on  both  types  of  aluminum.  A  secondary  crack  is 
visible  in  Figure  8  from  specimen  P3-7075.  Cleavage  characteristics  in  the  fatigue 
region  are  also  shown  on  this  figure.  The  fast  fracture  region  shows  ductile 
dimples. 

Because  testing  was  conducted  inside  of  an  SEM  considerable  fractographic 
characteristics  were  observed  during  testing.  Fractography  of  the  fracture  surface 
did  confirm  the  sites  of  crack  nucleation  at  pits  and  constituent  particles  and  showed 
the  shape  of  the  fatigue  region. 


Conclusions 

Based  on  the  results  presented,  the  following  conclusions  were  drawn: 

1 .  The  prior  corrosion  did  have  an  effect  on  the  short  crack  fatigue  crack 
growth  rate.  It  was  increased  from  the  base  tests. 

2 .  The  7075  specimens  had  faster  crack  growth  rates  than  the  2024  specimens. 
The  baseline  tests  and  the  prior  corroded  2024  specimen  displayed  crack 

growth  rates  that  increased,  then  decreased,  then  increased  again  as  AK 
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Figure  7.  Possible  striation  marks  in  specimen  Xl-2024. 


Figure  8.  Secondary  crack  and  ductile  dimples  on  the 
fracture  surface  of  specimen  P3-7075. 
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increased.  The  base  test  of  7075  material  did  not  show  as  dramatic  a  drop¬ 
off  in  crack  growth  rate  as  the  2024  tests.  The  prior  corroded  7075 
specimen  data  increased  but  did  not  decrease  as  the  other  specimen  did. 

3 .  The  damage  on  the  surface  of  the  prior  corroded  specimens  was  in  the  form 
of  corrosion  pits  and  exfoliation.  This  made  it  more  difficult  to  scan  the 
surface  of  the  material.  The  pitting  and  exfoliation  made  it  harder  to 
distinguish  cracks  from  mudcracking. 

4 .  Results  from  one  test  at  each  test  condition  are  not  conclusive  statistically. 
They  give  an  idea  of  the  behavior  but  should  be  backed  up  by  subsequent 
tests  that  are  statistically  planned. 

5 .  The  method  of  crack  measuring  provided  a  good  estimate  of  what  was 
taking  place. 
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SHORT  &  LONG  FATIGUE  CRACKS: 

Analysis  and  Tmnlications  to  Life  Prediction 
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Life  Prediction 


FATIGUE  DAMAGE  EVALUATION 
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RESIDUAL  STRESSES  STRESS  CONCENTRATIONS 
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Crack  Growth  Rate  Equation  with  Two  Parameters 
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Short  Crack  to  Long  Crack 
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Integration  with  Sensors  Data 
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Abstract 

Personnel  in  the  Quality  and  Integrity  Design  Engineering  Center  (QIDEC)  at  the 
University  of  Utah,  investigated  the  fatigue  characteristics  of  8090  T7  and  7075 
T7351  aluminum  alloys  under  fretting  conditions  in  air  and  vacuum  environments. 
The  purpose  of  this  investigation  was  to  provide  insights  into  the  effects  of 
corrosion,  specifically  oxidation  in  air,  on  the  fretting  fatigue  process.  As  a  result 
of  this  investigation,  it  is  concluded  that  crack  nucleation  and  early  propagation  are 
influenced  by  material,  environment,  and  local  conditions  in  the  area  of  relatively 
great  damage  where  the  crack  is  formed.  It  is  recommended  that  where  the  effects 
of  fretting  fatigue  are  an  issue,  the  local  conditions  that  influence  the  fretting 
fatigue  process  should  be  investigated.  Understanding  these  conditions  for  the 
materials  being  investigated  may  help  in  developing  procedures  for  alleviating  the 
effects  of  fretting  fatigue. 
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Introduction 


Life  reduction  due  to  fretting  is  having  significant  impact  on  Air  Force 
systems.  It  is  a  major  factor  causing  the  high  cycle  fatigue  problems  currently 
being  investigated  by  the  Air  Force.  This  paper  discusses  some  of  the  factors  that 
influence  fretting  and  the  fretting  fatigue  process. 

ASTM  defines  fretting  as  "a  wear  phenomenon  occurring  between  two 
surfaces  having  oscillatory  relative  motion  of  small  amplitude."^  In  air  the  process 
includes  oxidation  of  the  wear  debris  which  is  trapped  between  the  two  surfaces. 
When  fatigue  loading  is  present  the  combined  fretting  fatigue  process  can  result  in 
a  reduction  in  fatigue  life  of  more  than  an  order  of  magnitude  compared  to  fatigue 
without  fretting. 

The  fretting  fatigue  test  system  used  in  this  experimentation  was  developed 
within  QIDEC  to  enable  testing  either  in  air  or  a  scanning  electron  microscope 
(SEM)  vacuum  environment.  This  allows  investigation  of  fretted  surfaces  of 
specimens  tested  in  the  SEM  without  the  SEM  vacuum  being  broken.  Therefore, 
fretting  fatigue  damage  of  specimens  tested  in  the  SEM  can  be  considered  as 
resulting  from  wear  only,  whereas  damage  from  tests  in  air  results  from  the 
concurrent  wear  and  oxidation  mechanisms.  A  top-view  picture  of  the  fretting 
fatigue  load  frame  is  shown  in  Figure  1,  and  the  load  frame  mounted  in  the  SEM  is 
shown  in  Figure  2.  Additional  information  concerning  the  test  system  is  in 
references  2  and  3. 

The  experimentation  was  conducted  to  develop  cycles-to-failure  fretting 
fatigue  data  for  the  8090  T7  and  7075  T7351  aluminum  alloys  under  various 
conditions  in  air  and  vacuum  environments  and  to  gain  insights  into  the  effects  of 
corrosion  on  the  fretting  fatigue  process.  To  address  the  second  objective,  that  of 
gaining  insights,  fractographic  investigations  were  conducted  of  specimens  that 
had  failed  during  testing.  Additionally,  tests  of  eighteen  specimens  were  stopped 
prior  to  failure  so  that  the  fretted  surfaces  could  be  investigated  in  the  SEM  and  the 
specimens  could  be  sectioned  for  metallographic  investigation  to  gain  insights  into 
the  crack  nucleation  process.  This  paper  focuses  on  the  SEM,  metallographic,  and 
fractographic  investigations.  The  cycles-to-failure  data  are  summarized  in  this 
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paper  to  enable  the  reader  to  put  the  pictographic  results  in  perspective.  They  are 
reported  in  reference  2  and  in  the  paper  submitted  for  publication  as  noted  in 
references  4. 


Figure  1.  A  top-view  picture  of  the  fretting  fatigue  load  frame.  The  vertical  plate 
in  the  center  acts  as  the  SEM  lid.  The  part  to  the  left  of  this  plate  fits  into  the  SEM. 


Results  of  the  Investigation 


The  cycles-to-failure  data  from  this  program  are  summarized  in  Table  1  for 
each  material.  As  expected,  cycles-to-failure  decrease  when  the  wear  mechanism 
of  fretting  in  a  vacuum  is  added  to  fatigue,  with  a  further  decrease  when  the  testing 
is  conducted  in  air  so  that  oxidation  also  is  a  factor.  Note  the  influence  of  the 
materials.  For  the  8090-T7  alloy  the  wear  mechanism  is  relatively  important 
compared  to  the  corrosion  mechanism.  However,  for  the  7075-T7351  alloy  the 
corrosion  mechanism  is  dominant  in  reducing  cycles-to-failure  and  the  overall 
reduction  is  greater. 


Table  1.  Material  Average  Cycles-to-failure  Date 


Conditions  of  testing 

8090-T7 

7075-T7351 

Fatigue  without  fretting 

1,759,527 

2,537,695 

Fretting  fatigue  in  a  vacuum 

403,223 

2,377,757 

(assumed  wear  only) 

Fretting  fatigue  in  air 

341,980 

106,706 

(oxidation  and  wear) 

An  example  of  the  fretted  surface  of  a  specimen  tested  in  air  is  shown  in 
Figure  3.  Prior  to  testing,  the  surfaces  of  this  and  all  other  specimens  were  sanded 
and  polished  with  the  final  polish  being  1  micron  diamond  compound.  The 
direction  of  the  fatigue  loading  in  all  SEM  pictures  of  fretted  surfaces  is  towards 
the  upper  right  comer  of  the  picture.  The  fretting  fatigue  cracks  run  generally 
perpendicular  to  the  direction  of  fatigue  loading  from  upper  left  to  lower  right  in 
SEM  pictures  of  fretted  surfaces.  This  fretted  surface  is  characteristic  of  the 
surfaces  of  all  specimens  tested.  Numerous  areas  where  relatively  great  damage 
has  occurred  can  be  seen  in  the  upper  right  portion  of  the  fretting  pattern.  These 
are  the  areas  where  cracks  tended  to  nucleate.  Most  of  the  fretted  surface  pictures 
which  follow  wiU  be  of  such  areas  of  relatively  great  damage. 

Figure  4  shows  a  crack  that  has  nucleated  in  and  propagated  through  the 
fretted  surface  of  a  specimen.  Although  generally  perpendicular  to  the  direction  of 
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fatigue  loading,  locally  the  crack  propagated  along  the  contour  of  the  heavy  fretting 
so  that  at  each  point  along  its  path  it  tended  to  circle  the  dominant  local  area  of 
relatively  great  damage.  The  influence  of  local  fretting  conditions  on  crack 
nucleation  and  early  propagation  was  apparent  throughout  the  investigation. 


Figure  3.  An  example  of  the  surface  of  a  specimen  tested  in  air.  The  fretted  area 
is  shown  in  the  center  of  the  picture.  The  direction  of  fatigue  loading  is  towards 
the  upper  right  comer  of  the  picture.  (17x  magnification) 

Figures  5  and  6  show  two  8090  specimens  tested  under  identical 
conditions  and  for  the  same  number  of  cycles  except  that  the  specimen  in  Figure  5 
was  tested  in  the  SEM  vacuum  and  the  specimen  in  Figure  6  was  tested  in  air.  The 
effects  of  adhesive  wear  and  smearing  are  apparent  for  the  specimen  in  Figure  5, 
but  the  damage  is  less  severe  than  that  of  the  specimen  in  Figure  6.  The  smearing 
was  a  characteristic  of  specimens  tested  in  the  SEM  for  both  materials.  The  more 
severe  damage  was  characteristic  of  specimens  tested  in  air  for  both  materials.  A 
crack  is  propagating  through  the  damaged  area  on  the  surface  of  the  specimen  in 
Figure  6.  Figure  7  shows  a  cross  section  through  this  area  with  the  specimen 
material  on  the  bottom  of  the  picture.  Cracks  to  the  left  and  right  are  growing 
towards  each  other  and  might  be  expected  to  enlarge  the  pit  as  opposed  to 
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propagating  to  failure.  One  branch  of  a  crack  on  the  left  side  of  the  picture  appears 
to  have  started  to  grow  into  the  material. 


Figure  4.  The  fretted  surface  of  a  specimen  showing  a  meandering  crack  that  at 
any  point  is  generally  circling  the  dominant  local  area  of  relatively  great  damage. 
(150x  magnification) 


Figure  5.  The  fretted  surface  of  an  8090  specimen  which  was  tested  in  a  vacuum. 
The  effects  of  adhesive  wear  and  smearing  are  apparent,  but  the  damage  is  less 
severe  than  that  seen  on  specimens  tested  in  air  such  as  the  specimen  in  Figure  6. 
(150x  magnification) 
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Figure  6.  A  crack  propagating  from  a  pit  on  the  surface  of  an  8090  specimen 
which  was  tested  in  air.  A  cross  section  through  this  area  of  damage  is  shown  in 
Figure  7.  (470x  magnification) 


Figure  7.  A  cross  section  through  the  damaged  area  shown  in  Figure  6.  Cracks  to 
the  left  and  right  are  growing  towards  each  other.  One  branch  of  a  crack  on  the 
left  side  appears  to  have  started  to  grow  into  the  material.  (490x  magnification) 


Figure  8  shows  the  fretted  surface  of  a  7075  specimen  tested  in  air. 
Cracking  can  be  seen  in  each  of  the  more  severe  areas  of  damage  in  this  picture.  In 
each  case  the  cracking  is  oriented  towards  the  center  of  the  localized  damage. 
Figure  9  shows  a  cross  section  through  an  area  of  damage  on  this  specimen.  As 
with  the  8090,  specimen,  cracks  to  the  left  and  right  are  growing  towards  each 
other  and  might  be  expected  to  enlarge  the  pits.  However,  for  the  7075  material 
the  cracking  is  closer  to  the  surface  of  the  specimen. 


Figure  8.  The  fretted  surface  of  a  7075  specimen.  Three  regions  of  localized 
damage  can  be  seen.  Each  has  cracks  that  are  oriented  towards  the  center  of  the 
localized  damage.  (270x  magnification) 

Figures  10  and  1 1  show  cross  sections  of  the  fretted  surfaces  of  an  8090 
and  a  7075  specimen  respectively.  These  specimens  were  cut  so  that  their  long 
axes  were  across  the  thickness  of  the  rolled  plates  of  material.  This  meant  that  the 
fatigue  loading  would  tend  to  cause  crack  propagation  along  the  grain  direction 
near  the  center  axis  of  the  plate  of  material.  The  tendency  for  the  cracks  to  connect 
causing  pits  or  troughs  is  evident  for  the  8090  material  in  Figure  10  as  it  was  for 
both  materials  in  the  previous  figures  where  specimens  were  cut  with  long  axes 
parallel  to  the  surfaces  of  the  plates  of  material.  However  the  crack  in  the  7075 
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material  is  oriented  along  the  grain  direction  and  propagating  directly  into  the 
material.  That  this  would  cause  more  rapid  propagation  to  failure  in  the  7075 
material  was  reflected  in  the  test  parameters.  During  the  pretest  it  was  determined 
that  the  7075  specimen  life  at  the  planned  fatigue  loading  would  not  provide 
sufficient  cycles-to-failure  to  allow  meaningful  test  results.  Therefore  the  fatigue 
loading  on  the  through-the-thickness  7075  specimens  was  reduced. 


Figure  9.  A  cross  section  through  a  damaged  area  shown  in  Figure  8.  Cracks 
from  the  surface  tend  to  join,  but  in  regions  closer  to  the  surface  than  was 
observed  with  8090  specimens  such  as  shown  in  Figure  7.  (990x  magnification) 

The  bottom  half  of  Figure  12  shows  the  fracture  surface  of  an  8090 
specimen  tested  in  air.  A  nucleation  site  and  several  levels  of  cracks  can  be  seen  at 
the  center  of  the  picture.  These  are  further  from  the  viewer  than  the  features  at  the 
very  bottom  of  the  picture.  It  is  as  if  the  viewer  is  looking  at  the  nucleation  site 
from  the  top  of  a  cliff.  The  fracture  surface  of  the  7075  specimen  shown  in  the 
bottom  half  of  Figure  13  has  more  features  but  less  relief  and  fewer  facets  than  the 
8090  specimen  in  Figure  12.  Again,  this  would  indicate  that  different  localized 
mechanisms  of  crack  nucleation  and  early  propagation  are  involved  for  different 
materials. 
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Figure  10.  Cross  section  through  a  damaged  area  on  the  surface  of  a  through- the- 
thickness  8090  specimen.  The  tendency  of  cracks  to  connect  causing  pits  or 
troughs  is  evident.  (990x  magnification) 


Figure  11.  Cross  section  through  a  damaged  area  on  the  surface  of  a  through-the- 
thickness  7075  specimen.  The  tendency  of  cracks  to  connect  causing  pits  or 
troughs  is  not  evident.  Instead,  the  crack  is  oriented  to  the  grain  direction  and 
propagating  directly  into  the  material.  (470x  magnification) 
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Figure  12.  The  fracture  surface  of  an  8090  specimen  showing  a  nucleation  site 
and  several  levels  of  cracking.  These  are  further  from  the  viewer  than  the  features 
at  the  bottom  of  the  picture.  It  is  as  if  the  viewer  is  looking  at  the  nucleation  site 
from  the  top  of  a  cliff.  (7  lx  magnification) 


Figure  13.  A  partial  pit  can  be  seen  on  the  fracture  surface  of  a  7075  specimen 
which  has  more  features,  but  less  relief  and  fewer  facets  than  the  fracture  surface 
of  8090  specimens  such  as  that  in  Figure  12.  (78x  magnification) 
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Conclusions 


Based  on  the  fretting  surface,  metallographic,  and  fractographic 
investigations  conducted  in  this  research  program,  it  is  concluded  that: 

•  Areas  of  relatively  great  damage  on  a  fretted  surface  can  be  caused  by 
adhesive  wear. 

•  These  areas  have  more  severe  damage  if  corrosive  mechanisms  also  are 
active. 

•  Cracks  nucleate  at  or  near  surfaces  of  areas  of  relatively  great  damage  and 
propagate  within  and  from  these  areas. 

•  For  specimens  tested  in  air,  these  cracks  can  join  other  cracks  or  return  to 
the  surface,  forming  pits  or  troughs  and  debris,  or  propagate  into  the 
material. 

•  Cracks  also  propagate  into  the  material  from  the  bottoms  of  the  pits  and 
troughs. 

•  Crack  nucleation  and  early  propagation  are  influenced  by  material, 
environment,  and  local  conditions  in  the  area  of  relatively  great  damage 
where  the  crack  is  formed. 

Recommendations 

Where  the  effects  of  fretting  fatigue  are  an  issue,  the  local  conditions  that 
influence  the  fretting  fatigue  process  should  be  investigated.  Understanding  these 
conditions  for  the  materials  being  investigated  may  help  in  developing  procedures 
for  alleviating  the  effects  of  fretting  fatigue. 
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1.  INTRODUCTION 

The  development  and  evaluation  of  models  to  predict  fatigue  crack  growth  in  metallic 
structures  under  spectrum  loading  is  a  topic  which  has  received  a  great  deal  of 
attention  in  the  last  20  years.  References  1  to  9  detail  some  of  the  studies  which  have 
been  undertaken.  There  are  a  large  number  of  computer  programs  incorporating 
various  models  for  predicting  fatigue  crack  growth.  The  author  is  aware  of  at  least  12 
(Broek,  Cracks,  Cracks  IV,  Cracks  84,  MODGROW,  AFGROW,  CORPUS,  ONERA, 
PREFFAS,  FASTRAN,  STRIPY,  CG90).  The  various  programs  and  models  fall  into 
one  of  two  categories  as  follows: 

a.  Empirical.  Empirical  models  employ  various  techniques  for  modelling  the 
baseline  constant  amplitude  crack  growth  rate  data  (eg  tabular  inputs  which  are 
then  subject  to  interpolation/extrapolation,  and  equations  such  as  Forman  or 
Walker  which  fit  a  curve  to  the  measured  data).  Empirical  load  interaction  models 
with  “calibration”  constants  are  then  used  to  account  for  effects  such  as  retardation 
and  therefore  match  a  prediction  with  observed  or  measured  behaviour. 

b.  Analytical.  Some  models  utilise  the  crack  closure  theory  to  explain  the  R-ratio  and 
load  interaction  effects  and  apply  this  to  the  prediction.  A  certain  degree  of 
empiricism  still  exists,  however  the  method  by  which  the  R-ratio  and  load 
interaction  effects  are  accounted  for  is  based  on  a  more  scientific  theory  than 
simply  curve  fitting  observed  behaviour  or  incorporating  an  otherwise  meaningless 
constant  which  can  be  adjusted  to  give  a  desired  result. 

The  aim  of  the  work  presented  in  this  paper  was  to  assess  the  performance  of  the 
empirical  models  compared  with  the  analytical  closure  model  approach.  Test  data 
was  available  for  relatively  simple  centre  crack  specimens  machined  from  7075-T651 
Aluminium  Alloy  plate  material  (Reference  10).  The  specimens  were  pre-cracked  and 
then  subjected  to  a  simplified  variable  amplitude  fighter  aircraft  load  spectrum.  The 
peak  load  in  the  baseline  spectmm  (equal  to  a  eg  acceleration  of  N^  =  7.5  g)  was 
approximately  28%  of  the  yield  stress  of  the  material.  The  spectrum  was  altered  to 
simulate  placard  flight  restrictions  which  may  be  imposed  to  conserve  the  fatigue  life 
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of  an  aircr^.  The  restrictions  (at  6.5  g  and  5  g,  see  Figures  2  and  3)  and  an  artificial 
peak  load  increase  (to  8.5  g,  see  Figure  3)  resulted  in  significant  changes  to  the  crack 
growth  life  which  would  not  be  predicted  by  a  model  which  does  not  account  for  load 
interaction  effects. 


2.  SPECIMEN  TESTING 

The  specimen  test  results  used  in  this  work  were  obtained  from  Reference  10.  The 
specimens  were  machined  from  7075-T651  Aluminium  Alloy.  The  geometry  used 
was  a  fixed  end,  centre  crack  specimen  as  shown  in  Figure  1  below. 


Figure  1.  Fixed  End  Centre  Crack  7075-T651  Specimen 
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A  central  slot  was  machined  into  the  specimens  using  Electrical  Discharge  Machining 
(EDM)  and  they  were  then  subjected  to  constant  amplitude  pre-cracking  until  a  crack 
length  of  0.15  inch  (a=0.15  inch,  2a=0.30  inch)  was  reached.  The  final  maximum 
load  during  pre-cracking  was  3,000  pounds.  The  spectrum  loading  applied  later  was 
such  that  this  load  level  was  exceeded  in  the  first  10  cycles.  Residual  effects  from 
pre-cracking  were  therefore  not  anticipated. 

The  specimens  were  subjected  to  four  variations  of  a  simplified  Mirage  fighter  aircraft 
spectrum; 

(i)  the  unmodified  Mirage  spectrum, 

(ii)  clipped  at  6.5  g, 

(iii)  clipped  at  5  g,  and 

(iv)  with  the  7.5  g  peak  (which  occurs  once  in  a  100  flight,  66.6  hour  block) 
increased  to  8.5  g. 

All  loading  variables  other  than  the  spectrum  itself  were  kept  constant  throughout  the 
tests  in  order  to  isolate  the  effect  of  spectrum  modification.  A  scaling  factor  of  remote 
stress  per  g  of  2.5  ksi  per  g  was  used,  The  spectrum  is  detailed  in  Figure  2  and  3,  and 
Tables  1  and  2  below. 
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Figure  2.  Representative  Segments  From  the  Unmodified  and  6.5  g  Limit  Spectra 
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8-5g  PEAK  SPECTRUM 

Figure  3.  Representative  Segments  From  the  5  g  Limit  and  8.5g  Peak  Spectra 
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Table  1 .  Flight  Segment  Spectra 
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Table  2.  Sequence  of  Flight  Segments  in  a  Block  Representing  100  Flights,  66.6 

Hours  (1989  cycles) 


Three  specimens  were  subjected  to  the  unmodified  spectrum,  two  specimens  to  the  6.5 
g  limit  spectrum,  three  specimens  to  the  5  g  limit  spectrum  and  one  specimen  to  the 
8.5  g  limit  spectrum.  All  specimens  were  tested  to  failure  which  occurred  at  a  half 
crack  length  a  of  between  0.5  and  0.6  inch.  The  crack  growth  curves  are  plotted  in 
Figure  4  below.  There  is  very  little  scatter  in  the  results  where  several  tests  were 
performed  under  the  same  spectrum.  It  is  therefore  considered  reasonable  to  conclude 
that  the  differences  in  crack  growth  rate  can  be  attributed  solely  to  the  differences  in 
the  load  spectra. 


I  ♦  Unmodified  Spectrum: 
laSg  Limit  Spectrum 
A6.5g  LimitSpectnjm 
X  8,5  g  Limit  Spectrum 


20.000  40,000  60.000  80.000  100,000  120.000  140.000  160.000  180.000 

Load  Cycles 

Figure  4.  Experimental  Crack  Growth  Results  for  Four  Mirage  Spectra  Variations 

(Reference  10) 
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The  crack  growth  results  clearly  demonstrated  the  load  interaction  effects  expected 
under  spectrum  loading.  When  the  peak  load  was  clipped  at  6.5  g,  there  was  a  loss  of 
retardation  effect  and  the  crack  growth  life  was  shorter.  Clipping  at  5  g  led  to  an  even 
shorter  life.  Increasing  the  peak  load  from  7.5  to  8.5  g  (a  load  which  occurs  once  per 
block  of  1989  cycles)  doubled  the  life.  These  results  were  therefore  able  to  be  used  to 
assess  the  adequacy  of  the  various  models  to  predict  the  load  interaction  effects. 


3.  PREDICTIVE  MODELS 

The  predictive  models  examined  here  can  be  classified  either  as  empirical 
(CRACKS84,  Broek,  AFGROW)  or  based  on  the  analytical  closure  model 

(FASTRAN  II).  A  description  of  the  significant  features  of  each  is  presented  in  this 
section. 

3.1  CRACKS84 

A  program  called  CRACKS84  was  used  in  the  Reference  10  work.  The  results  are 
reproduced  here  for  comparison  purposes.  CRACKS84  was  based  on  a  program 
called  CRACKGRO  and  was  obtained  from  the  USAF  and  modified  to  run  on  a 
Gould  9080  computer  running  a  UNIX  operating  system.  CRACKS84  is  based  on  the 
standard  Linear  Elastic  Fracture  Mechanics  (LEFM)  approach  to  fatigue  crack 
propagahon.  The  code  includes  stress  intensity  factor  solutions  for  various 
geometries,  several  different  crack  growth  rate  models  and  several  load  interaction 
effect  models.  Reference  1 0  gives  full  details  on  the  program  and  how  it  was  applied 
in  this  case.  The  load  interaction  models  used  were  none,  Basic  Willenborg, 
Generalised  Willenborg  and  Willenborg/Chang. 


3.2  Broek 

A  suite  of  software  has  been  developed  by  Dr  D.  Broek  from  FracfuREsearch  Inc  in 
the  USA  (Reference  11).  The  software  also  utilises  the  LEFM  approach.  The 
retardation  models  examined  were  none,  Willenborg,  Calibrated  Wheeler  and 
Calibrated  Broek.  The  software  is  windows  based  and  runs  on  a  PC. 


3.3  AFGROW 

AFGROW  is  a  fatigue  crack  growth  software  package  developed  by  the  US  Air  Force 
at  Wright  Patterson  Air  Force  Base.  OHIO.  USA  (Reference  12).  AFGROW  is  a 
workstation  based,  graphically  interactive  computer  program  for  simulation  of  fatigue 
crack  growth  in  common  structural  geometries  subject  to  spectrum  loading.  It  is 
designed  to  run  either  on  a  PC  under  Windows  95,  or  under  a  UNIX  operating  system. 
For  the  current  work  the  PC  Windows  95  version  was  used.  The  program  uses  the 

LEFM  approach  and  the  retardation  models  examined  were  none,  Willenborg 
Wheeler  and  Closure. 
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3.4  FASTRAN II 


FASTRAN  II  is  a  PC  based  life  prediction  code  based  on  the  crack  closure  concept 
and  is  used  to  predict  crack  length  against  cycles  from  a  specified  initial  crack  size  to 
failure  for  many  common  crack  configurations  found  in  structmal  components.  The 
model  is  based  on  plasticity  induced  fatigue  crack  closure  and  is  used  to  calculate  the 
stress  level  at  which  the  crack  tip  becomes  fully  open  during  cyclic  loading 
(Reference  13),  leading  to  the  concept  of  an  effective  applied  stress  intensity  range 
(AKeff).  The  model  is  expected  to  provide  a  more  accurate  measure  of  the  crack 
growth  rate  for  a  spectrum  loading  situation. 


4.  INPUT  DATA 

4. 1  Common  Inputs 

The  inputs  which  were  common  to  all  programs  were  as  follows: 

a.  Specimen  Width  =  1 .5  inch 

b.  Specimen  Thickness  =  0.2  inch 

c.  Initial  Crack  Length,  Uj  =  0. 1 5  inch 

d.  Final  Crack  Length,  a,  =  0.60  inch 

e.  Specimen  Geometry:  Finite  width  centre  crack  sheet 
f  Material:  7075-T651  Aluminium  Alloy 

g.  Loading:  Uniaxial  Tension 

h.  Spectra:  Expressed  in  terms  of  “g”  with  a  scale  of  2.5  ksi  per  g 

i.  Material  Yield  Strength,  a  „  =  67  ksi 

j.  Plane  Strain  Fracture  Toughness,  K,(-  =  '})5ksi4inch 

4.2  CRACKS  84 

4.2.1  Crack  Growth  Rate  Data 

As  detailed  in  Reference  10,  the  crack  growth  rate  data  was  modelled  using  the 
Walker  equation  as  follows: 


671 


^=C,[  (  I-«)"''aA:]"  for  R>0.0 


-J^  =  CA(l*R‘rK„^r  fotR-^0.0 

Where: 

da 

—  ==  Crack  growth  rate 

AAT  =  Stress  intensity  factor  range 

^max  =  Maximum  stress  intensity  in  the  cycle 

C^,m,n,q  =  material  constants 

As  detailed  in  Reference  10,  the  available  material  data  (Reference  15)  was  modelled 
with  three  straight  line  segments  with  C„and  n  values  as  follows  (inch  and  ksi  units): 

Q  =4.9329x10''° 

n^  =  4.4533 

=2.4379x10'® 

=  2.60 

C„3=  7.3126x10'" 

«3  =5.1766 

A  value  of  w  =  0.6  was  determined. 

A  typical”  value  of  ^  =  1.0  was  used. 

A  value  of  threshold  stress  intensity  range,  AKj-h,  of  was  used. 

CRACKS  84  included  an  option  to  allow  to  change  as  a  function  of  stress  ratio, 
R.  This  option  was  not  used,  ie  was  kept  constant.  This  had  no  impact  on  the 
prediction  in-  any  case  since  the  smallest  AK  cycles  even  at  the  starting  crack  length 
exceeded  25ksiy]inch . 
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Positive  and  negative  stress  ratio  limits,  ie  values  of  R  above  and  below  which  stress 
ratio  is  presumed  not  to  have  an  effect  on  crack  growth  rate  were  required  to  be  input. 
“Typical  values”  of  +0.75  and  -0.99  were  used. 


4.2.2  Load  Interaction  Model  Input  Data 

The  load  interaction  models  used  in  CRACKS  84  (Basic  Willenborg,  Generalised 
Willenborg  and  Willenborg/Chang)  are  all  influenced  by  the  size  of  the  plastic  zone 
which  is  calculated  as  follows: 

__L 

•’'“yTC 

Where: 
r^.  =  radius  of  plastic  zone 

=  maximum  stress  intensity  due  to  the  current  load 

ma.x 

CT  j,,.  =  material  yield  stress 
Y  =  2  for  plane  stress,  6  for  plane  strain 

For  a  through  crack  case,  CRACKS  84  assumes  plane  stress  conditions. 

The  Generalised  Willenborg  model  introduced  another  parameter  known  as  the 
overload  shut-off  ratio,  R^.  A  typical  value  of  2.3  was  used. 


CT 


J 


4.3  BROEK 

4.3.1  Crack  Growth  Rate  Data 

The  Broek  software  includes  a  Walker  equation  option  for  modelling  the  crack  growth 
rate  data,  but  it  does  not  allow  the  user  to  input  the  equation  in  several  “regions”  each 
with  a  different  slope  (n)  and  intercept  (C)  value.  It  also  does  not  incorporate  a 
threshold  in  the  equation.  There  is.  however,  the  option  of  a  tabular  input  where  the 
crack  growth  rate  for  a  range  of  R  ratios  is  input  as  a  function  of  AAT .  The  Mirage 
spectra  consisted  of  sub  blocks  of  constant  amplitude  cycles  which  covered  a  limited 
range  of  R  ratios.  There  were  eight  R  ratios  in  total,  ie  R  =  -0.500,  -0.385,  -0.333,  - 
0.300,  -0.294,  0,  0.125  and  0.333.  The  Broek  software  could  only  accommodate  a 
maximum  of  six  R  ratios  in  the  table,  so  the  data  at  R  =  -0.333  and  R  =  -0.294  were 
left  out.  It  was  considered  that  the  results  would  not  be  significantly  affected  since 
there  was  data  at  fairly  close  R  ratios  in  the  table.  The  table  was  adjusted  so  that  an 
extremely  small  crack  growth  rate  corresponded  to  AK  levels  below  the  threshold 
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value  of  25ksiylinch  ,  thus  modelling  the  threshold  behaviour  in  a  similar  fashion  to 
CRACKS  84.  The  Walker  equation  representation  of  the  crack  growth  rate  relation  is 
shown  for  selected  R  ratios  in  Figure  5. 


--0.5 
-R=0  . 

-R=0.33i 


Delta  K  (ksi  Sqrt  Inch) 

Figure  5.  Walker  Equation  Model  of  Crack  Growth  Rate  Data  for  7075-T651 


4.4  AFGROW 

4.4. 1  Crack  Growth  Rate  Data 

AFGROW  was  able  to  accept  up  to  5  straight  line  segments  for  the  Walker  equation, 
so  the  three  sets  of  slope  and  intercept  (n  and  C)  as  per  CRACKS  84  were  used' 
AFGROW  does  not  request  a  value  of  the  parameter  “q”  (see  Section  4.2.1).  The 
program  deals  with  negative  R  ratio  data  by  using  the  “Harter  T  Method”  (Reference 
12).  It  does  prompt  for  positive  and  negative  stress  ratio  limits.  Although  +0.75  and  - 
0.99  were  used  with  CRACKS  84.  AFGROW  only  allows  up  to  +1.0  and  -0.5  on  the 
negative  side.  Because  the  lowest  R  ratio  in  the  Mirage  spectra  was  -0.5,  it  was 
considered  reasonable  to  set  the  negative  limit  at  -0.5  and  this  should  give  a  consistent 
result  with  the  CRACKS  84  predictions.  The  positive  limit  was  set  at  0.75  as  per 


4.5  FASTRAN  II 

4.5.1  Crack  Growth  Rate  Data 

The  FASTRAN  II  approach  requires  that  you  obtain  the  effective  stress  intensity 
range,  ,  as  a  function  of  crack  growth  rate  for  the  material  of  interest.  Newman 
(Reference  13)  has  developed  a  computer  program  to  process  raw  data  using  a  given 
constraint  factor,  a  ,  to  produce  the  data.  Data  was  available  from  Reference 
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14  for  7075-T6  material.  This  data  is  plotted  along  with  the  Walker  equation  data  at 
R=0  and  R=0.75  in  Figure  6  below.  The  Newman  data  was  considered  to  be 
consistent  with  the  other  data  and  was  used  to  perform  the  predictive  runs  with 
FASTRAN  II.  Some  modification  was  made  to  the  constraint  factor,  a  ,  and  this  is 
discussed  in  Section  5.6. 
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;  Data  from  Reference  1 4 
■  —^Walker  Equation  Data  at  R  =  0,75 
—A— Walker  Equation  Data  at  R=0 


Figured.  Comparison  of  Newman’s  7075-T6  data  with  Walker  Equation  Data 

5.  RESULTS 
5.1  No  Retardation 

The  first  case  to  be  considered  was  a  prediction  without  any  retardation  effects  being 
considered.  This  does  not  apply  to  FASTRAN  II  because  it  can  only  be  run  with  load 
interaction  effects  being  considered.  The  results  are  summarised  in  Table  3  below: 


675 


Spectrum 

Experimental 

Result 

Predictive  Model  j 

CRACKS  84 

BROEK 

AFGROW 

Unmodified 

85,630 

48,552 

49,422 

48,746 

6.5  g 

71,714 

48,945 

49,473 

48,887 

5.0  g 

44,006 

52,402 

52,426 

*  52,431 

8.5  g 

178,965 

46,567 

49,365 

48,592 

Table  3.  Comparison  of  Experimental  and  Predicted  Crack  Growth  Lives  (in  cycles) 
from  ai=  0.15  inch  to  Failure  (approximately  a  =  0.6  inch)  with  No  Retardation 

Modelling. 


These  results  are  considered  to  be  reasonably  consistent.  They  have  been  obtained  by 
essentially  using  the  same  material  data,  albeit  converted  into  a  different  format  to  suit 
each  particular  program.  In  each  model,  the  predicted  life  is  increased  slightly  for  the 
6.5  g  Md  5.0  g  spectra,  and  reduced  for  the  8.5  g  spectrum  (compared  to  the 

unmodified  spectrum).  This  is  expected  because  the  retardation  effects  are  not 
predicted. 


5.2  Willenborg  Retardation  Model 


The  results  when  the  Willenborg  retardation  model  was  selected  are  shown  in  Table  4 
below.  In  the  case  of  CRACKS  84  there  are  three  variations  of  the  Willenborg  model. 


Spectrum 

Experimental 

Result 

CRACKS  84 

Broek 

AFGROW 

Basic 

Willenborg 

Generalised 

Willenborg 

Willenborg  / 
Chang 

Unmodified 

85,630 

90,323 

66,457 

64,334 

99,969 

75,329 

6.5  g 

71,714 

80,377 

63,876 

64,335 

100,242 

75,655 

5.0  g 

44,006 

64,715 

56,482 

57,855 

114,721 

58,832 

8.5  g 

178,965 

122,144 

76,400 

66,450 

99,885 

76,402 

Table  4.  Comparison  of  Experimental  and  Predicted  Crack  Growth  Lives  (Cycles) 
from  a,  =  0.15  inch  to  Failure  (approximately  a  =  0.6  inch)  Using  the  Willenborg 

Retardation  Model 
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5.2.1  CRACKS  84,  Basic  Willenborg 


For  the  unmodified  spectrum,  the  prediction  with  this  retardation  model  (90,323 
cycles)  has  increased  significantly  from  the  no  retardation  case  (48,552  cycles)  to  be 
much  closer  to  the  experimental  result  (85,630  cycles).  Relative  to  that  result,  the 
model  predicts  a  shorter  life  for  the  6.5  g  spectrum  and  a  shorter  life  again  for  the  5  g 
spectrum.  A  longer  life  is  predicted  for  the  8.5  g  spectrum.  The  absolute  values  of 
the  predictions  are  not  accurate,  but  the  trend  follows  the  same  pattern  as  the 
experimental  results. 

5.2.2  CRACKS  84,  Generalised  Willenborg 

The  trends  here  are  similar  to  the  Basic  Willenborg  model. 


5.2.3  CRACKS  84,  Willenborg  /  Chang 

Once  again,  the  trends  are  in  the  expected  direction  but  the  absolute  values  of  the 
predictions  vary.  The  Basic  Willenborg  model  appears  to  have  given  a  better  result 
in  this  case  than  either  the  Generalised  Willenborg  or  the  Willenborg/Chang  model. 


5.2.4  Broek,  Willenborg 

As  per  the  CRACKS  84  predictions,  the  unmodified  spectrum  result  has  increased 
significantly  (49,422  to  99,969  cycles).  Relative  to  this,  the  predicted  lives  at  6.5  g 
and  5  g  have  increased  which  is  opposite  to  the  experimental  results  and  contradicts 
the  trends  observed  with  CRACKS  84.  The  predicted  life  at  8.5  g  has  decreased 
which  is  also  the  opposite  of  what  would  be  expected.  The  author  has  raised  concerns 
about  these  results  with  Dr  Broek  at  FractuREsearch  Inc  in  the  USA. 


5.2.5  AFGROW,  Willenborg 

As  per  the  previous  observations,  the  unmodified  spectrum  result  has  increased 
significantly  (from  48,746  to  75,329  cycles).  Relative  to  this,  the  6.5  g  spectrum 
result  is  a  slightly  longer  life  which  is  not_in  the  expected  direction.  The  5  g  result  is 
shorter  than  the  unmodified  spectrum  and  the  8.5  g  result  is  longer.  Both  of  these 
have  moved  in  the  expected  direction. 


5.3  Wheeler  Retardation  Model 

The  results  for  the  Wheeler  retardation  model  are  shown  in  Table  5  below: 
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Spectrum 

Experimental 

Result 

Wheeler  Model  | 

Broek 
m  =  0.835 

AFGROW 
m  =  0.86 

Unmodified 

85.630 

85,546 

85,274 

6.5  g 

71,714 

80,789 

80,313 

5g 

44.006 

61,043 

61,477 

8.5  g 

178,965 

104,193 

105,436 

Table  5.  Comparison  of  Experimental  and  Predicted  Crack  Growth  Lives  (Cycles) 
from  aj  =  0. 1 5  inch  to  Failure  (approximately  a  =  0.6  inch)  Using  the  Wheeler 

Retardation  Model 


In  this  case  the  trends  are  in  the  expected  direction  for  both  sets  of  software  and  the 
results  are  consistent  with  each  other.  The  calibration  parameters  of  m  =  0.835  and 
0.86  for  Broek  and  AFGROW  respectively  were  obtained  by  getting  a  “match”  for 
the  unmodified  spectrum. 


5.4  Broek  Retardation  Model 

The  Broek  retardation  model  was  calibrated  to  the  Unmodified  Spectrum  results  and 
the  predictions  are  shown  in  Table  6  below: 


Spectrum 

Experimental 

Result 

Broek  Model 

Factor  =  1.3 

Unmodified 

85,630 

87,860 

6.5  g 

71,714 

79,257 

5g 

44,006 

50,432 

8.5  g 

178,965 

206,477 

Table  6.  Comparison  of  Experimental  and  Predicted  Crack  Growth  Lives  (Cycles) 
from  a^  =  0. 1 5  inch  to  Failure  (approximately  a  =  0.6  inch)  Using  the  Broek 

Retardation  Model 


A  full  plot  of  these  crack  growth  predictions  compared  with  the  experimental  data  are 
shown  in  Figures  7  to  10  below. 
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Load  Cycles 

Figure  7.  Unmodified  Spectrum,  Broek  Retardation  Model  (Calibration  Factor  =  1.3), 

Experimental  Data  and  Prediction 


—Prediction 
A  ExperimentatOata 


Figure  8.  6.5  g  Limit  Spectrum,  Broek  Retardation  Model  (Calibration  Factor  =  1.3), 

Experimental  Data  and  Prediction 


Predicton 

■  ExpenmentalData 


Figure  9.  5  g  Limit  Spectrum,  Broek  Retardation  Model  (Calibration  Factor  =  1 .3), 

Experimental  Data  and  Prediction 


—  Prediction 
X  Experimental  Data 


Figure  10.  8.5  g  Limit  Spectrum,  Broek  Retardation  Model  (Calibration  Factor  = 

1 .3),  Experimental  Data  and  Prediction 

5.5  Closure  Model  (AFGROW) 

AFGROW  included  a  “Closure  Retardation  Model”  option,  but  it  appears  to  be 
significantly  different  to  the  analytical  model  used  in  FASTRAN  11.  The  model  relies 
on  the  user  inputting  the  Crack  Opening  Load  Ratio  (OLR)  at  R  =  0.  For  the  results 
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presented  here  the  default  value  of  OLR  =  0.15  was  used.  The  results  are  shown  in 
Table  7  below: 


Spectrum 

Experimental 

Result 

AFGROW  Closure  Model 
OLR  (R=0)  =  0.15 

Unmodified 

85,630 

45,792 

6.5  g 

71,714 

45,342 

5g 

44,006 

46,508 

8.5  g 

178,965 

47,716 

Table  7.  Comparison  of  Experimental  and  Predicted  Crack  Growth  Lives  (Cycles) 
from  ai=  0.15  inch  to  Failure  (approximately  a  =  0.6  inch)  Using  the  AFGROW 

Closure  Retardation  Model 


5.6  FASTRAN  II 

The  crack  growth  rate  data  used  with  the  FASTRAN  II  program  (see  Sect  4.5.1)  was 
obtained  from  centre  crack  specimens  of  2.3  mm  thickness.  The  specimens  tested 
under  the  Mirage  spectra  were  5  mm  thick  and  this  would  have  a  significant  effect  on 
the  constraint  factor,  a  .  The  constraint  factor  in  the  “constraint  loss”  regime 
(Reference  14)  was  adjusted  until  FASTRAN  II  gave  the  correct  prediction  for  the 
unmodified  spectrum.  This  was  essentially  a  “calibration”  of  the  da/dN  vs  data. 

The  “calibrated”  data  was  then  used  to  predict  the  crack  growth  under  the  other 
spectra.  The  a  conditions  necessary  were  as  follows: 

da/dN  <  2.76  x  1 0'^  inches  per  cycle,  a  =  1 .8 

da/dN  >  2.76  x  10"^  inches  per  cycle,  a  =1.57 

For  intermediate  rates,  a  was  varied  linearly  with  the  logarithm  of  crack  growth  rate 
(Reference  13). 
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The  crack  growth  predictions  obtained  are  shown  in  Table  8  below: 


Spectrum 

Experimental 

Result 

FASTRAN  II 
Prediction 

Unmodified 

85,630 

85,080 

6.5  g 

71,714 

73,182 

5g 

44,006 

46,276 

8.5  g 

178,965 

201,092 

Table  8.  Comparison  of  Experimental  and  Predicted  Crack  Growth  Lives  (Cycles) 
from  a,  =  0.15  inch  to  Failure  (approximately  a  =  0.6  inch)  Using  FASTRAN  II  with 

a  =  1 .57  in  the  constraint  loss  regime 


The  full  plot  of  the  crack  growth  predictions  for  FASTRAN  II  compared  with  the 
experimental  data  are  shown  in  Figures  1 1  to  14  below. 


♦  Experimental  Data 
—Prediction 


Figure  11.  Unmodified  Spectrum,  FASTRAN  II,  Experimental  Data  and  Prediction 
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Crack  Length  (inches)  uraui^  Lcngui 


A  Expen  mental  Data 
——Prediction 


igure  12.  6.5  g  Limit  Spectrum,  FASTRAN II,  Experimental  Data  and  Prediction 


■  ExpenmentalData 
—Prediction 


Figure  13.  5  g  Limit  Spectrum.  FASTRAN  11,  Experimental  Data  and  Prediction 


X  ExpenmentalOata' 
—— Predtction 


Figure  14.  8.5  g  Limit  Spectrum,  FASTRAN II,  Experimental  Data  and  Prediction 


6.  DISCUSSION  AND  CONCLUSIONS 

A  range  of  fatigue  crack  growth  prediction  models  and  a  range  of  software  packages 
have  been  evaluated  for  their  ability  to  deal  with  load  sequence  effects  for  long  fatigue 
crack  test  data.  The  test  data  utilised  clearly  demonstrated  a  load  sequence  interaction 
effect.  Observations/  comments  are  as  follows: 

a.  The  empirical  models  all  gave  similar  results  for  the  no  retardation  case.  As 
expected,  the  spectrum  changes  did  not  alter  the  predictions  greatly. 

b.  The  Willenborg  model  produced  inconsistent  and  illogical  results  in  the  Broek  and 
AFGROW  packages,  and  more  sensible  results  with  CRACKS  84.  There  appear  to 
be  many  variations  with  the  Willenborg  model,  and  the  reason  for  the 
inconsistencies  could  only  be  ascertained  by  examining  the  source  code  for  the 
various  programs.  The  result  of  this  study  will  be  conveyed  to  the  authors  of  the 
Broek  software  and  AFGROW  because  these  programs  are  currently  in  use. 

CRACKS  84  is  obsolete,  and  it  is  not  (to  the  author’s  knowledge)  being  actively 
used. 

c.  AFGROW  and  Broek  gave  reasonably  consistent  results  for  the  Wheeler 
retardation  model.  The  model  correctly  predicted  the  direction  in  which  the  crack 

growth  would  change  with  a  spectrum  change,  but  the  ratio  of  predicted  lives  was 
not  well  predicted. 

d.  The  calibrated  Broek  retardation  model  gave  good  results  when  you  examine  the 
final  life  tables  (Table  6).  However  the  plots  (Figures  7  to  10)  show  that  the  total 
curve  is  not  always  well  modelled.  The  prediction  for  the  8.5  g  peak  load  case 
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(Figure  10)  in  particular  exhibits  some  unusual  bumps.  These  bumps  appear  to  be 
more  pronounced  when  the  spectrum  contains  higher  peak  loads.  The  prediction 
for  the  5  g  case  is  a  smooth  curve  with  little  indication  of  bumps.  The  bumps  are 
not  a  phenomenon  which  occur  in  the  experimental  results. 

e.  The  AFGROW  closure  model  did  not  give  good  results.  The  5  g  result  goes 
against  the  expected  trend,  and  the  magnitude  of  the  variation  in  growth  life 
compared  with  the  unmodified  spectrum  is  significantly  less  than  observed  in  all 
cases.  Further  investigation  is  required. 

f  FASTRAN  II  produced  the  best  predictions  (see  Figures  11  to  14).  An 
encouraging  aspect  to  this  is  that  there  is  a  physical  explanation  to  the  “calibration” 
process  which  took  place.  The  baseline  crack  growth  data  had  been  obtained  from 
thinner  (2.3  mm)  specimens.  The  thicker  specimens  used  in  the  Mirage  spectrum 
testing  would  be  expected  to  have  a  higher  constraint  level.  Having  calibrated  the 
model  to  the  unmodified  spectrum,  very  good  predictions  were  obtained  for  the 
other  spectra. 

The  results  of  this  exercise  demonstrate  the  superior  performance  of  Newman’s 
analytical  crack  closure  model  for  predicting  fatigue  crack  growth  under  spectrum 
loading.  This  is  consistent  with  the  findings  in  References  5  and  6.  The  empirical 
nature  of  the  other  models  leads  to  inconsistent  and  illogical  results.  These  results  are 
for  a  relatively  simple  2-D  crack  configuration  and  the  models  would  not  be  expected 
to  perform  any  better  for  a  more  complex  geometry  or  more  complicated  loading. 
There  also  appears  to  be  inconsistency  in  how  the  models  have  been  coded  into  the 
software,  particulairly  for  the  Willenborg  retardation  model. 

The  disadvantage  with  the  FASTRAN  II  approach  is  the  issue  of  the  baseline  crack 
growth  rate  data.  However,  this  exercise  has  shown  that  as  observed  by  Newman 
(References  13  and  14),  data  from  high  R-ratios  can  be  used  to  approximate  the  da/dN 
vs  relation  (see  Figure  6).  The  constraint  factor,  a  ,  can  be  used  as  a  calibration 

factor.  It  is  not,  however,  simply  an  empirical  or  arbitrary  calibration.  The  constraint 
factor  needs  to  be  adjusted  to  reflect  the  level  of  constraint  present  in  the  particular 
case  being  analysed  when  compared  to  the  constant  amplitude  test  data.  There  will 
always  be  some  degree  of  compromise  inherent  in  this  because  the  constraint  varies 
along  any  crack  front  (ie  whether  it  is  a  through  crack  or  a  part  through  crack).  The 
analytical  crack  closure  model  is  considered  to  be  the  best  option  for  predicting 
fatigue  crack  growth  under  spectrum  loading. 
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Crack  Propagation  Load  Measurement  Method  (CPLM) 
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Atlantic  City,  NJ  08405 
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DETERMINING  FLIGHT  LOADS 
AND  CRACK  GROWTH  RATES 
FROM  FAILED  AIRCRAFT 
STRUCTURAL  COMPONENTS 

Donald  A.  Shockey,  Takao  Kobayashi, 
Charles  G.  Schmidt,  and  Richard  W.  Klopp 
SRI  International,  Menlo  Park,  CA  94025 

and 

Thomas  H.  Flournoy 
Federal  Aviation  Administration 
Technical  Center,  Atlantic  City,  NJ  08405 

SUMMARY — The  possibility  of  deducing  load  spectrum 
parameters  from  fatigue  failure  surfaces  is  explored  by 
applying  innovative,  three-dimensional  topographic  char¬ 
acterization  and  analysis  techniques  to  failure  surfaces  in 
aluminum  sheet.  Precise,  high-resolution  elevation  maps  of 
fracture  surfaces  were  obtained  using  confocal  optics 
scanning  laser  microscopy.  Elevation  power  spectral  den¬ 
sity  curves  resulting  from  a  fast  Fourier  transform  of  the 
elevation  data  appear  sensitive  to  stress  intensity  range 
and  environment.  A  conjugate  fracture  surface  matching 
procedure,  FRASTA,  can  detect  and  may  provide  a  way  to 
quantify  overloads. 

""  Poulter  Laboratory 


Gray-scale  images  and  Perspective  Views 
of  Topography  of  Conjugate  Fracture  Surfaces 
Produced  Under  Low  and  High  AK  Fatigue  Conditions 
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(b)  HIGHAK  SURFACES 
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Fast  Fourier  Transform  Power  Spectra 
for  Conjugate  Fatigue  Failure  Surfaces 
Produced  Under  High  and  Low  AK 
Loading  Conditions 


WAVELENGTH 


1mm  lOOum  lO^im 
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International 


EPSD  Curves  for  Two  Specimens  Tested 
Under  Identical  Conditions 


WAVELENGTH 


1mm 
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1mm 
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4,4  MPa  Vm 
1 1  MPa  Vm 
0.6 
5  Hz 


SPATIAL  FREQUENCY  (m'') 
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Peak  Stress  Intensity  Versus 
Crack  Length  for  Fatigue  Experiment 

Showing  Overloads 


CRACK  EXTENSION  AS  MEASURED 
FROM  BASE  OF  NOTCH  (in.) 

CAM-5048-83 


Poulter  Laboratory 


PEAK  STRESS  INTENSITY  (MPaVm) 


Lffilmlj 


A  Series  of  Plan  Views 
of  a  Crack  Front  at  Increasing 
Topograph  Displacements 

(White  areas  indicate  separated  materiai; 
biack  areas  indicate  intact  material) 
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CONJUGATE  SURFACE  SPACING  (1  =  936.9  urn) 

CAM-5048-87 
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CRACK  EXTENSION  (mm) 


CONJUGATE  SURFACE  SPACING  (1  *  936.9  ^m) 


CAM-604M8 


Poulter  Laboratory 


CRACK  EXTENSION  (mm) 


AREA  (%) 


Fracture  Progression  Curve 
Relation  to  Overioads  Shown 
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CRACK  EXTENSION  (mm) 


♦  Important  but  previously  unobtainable 
information  regarding  the  conditions  that 
produced  fatigue  faiiures  can  iikeiy  be 
obtained  from  fracture  surfaces. 

♦  Keys  to  extracting  this  information  are  a  fast 
and  accurate  method  for  characterizing 
fracture  surface  topography  and  a  three- 
dimensional  analysis  of  the  data. 

♦  Fourier  analyses  to  fracture  surfaces  may 
provide  the  basis  for  a  library  of  reference 
curves  useful  in  determining  the  load 
spectrum  and  environment  conditions  that 
caused  a  service  failure. 


Poulter  Laboratory 
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♦  Comparison  of  conjugate  fracture  surface 
topographs  may  indicate  and  quantify 
details  of  the  overload  spectrum 
experienced  by  a  component  that  failed  in 
fatigue. 

♦  Shapes  of  striations  may  be  quantifiable  by 
stereoscopy  or  COSLM  and  provide  a  means 
to  determine  maximum  and  minimum  values 
of  the  cyclic  load  spectrum. 

♦  These  expectations  should  be  explored  to 
seek  advances  in  failure  analysis 
technology. 


Poulter  Laboratory 
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Full-Scale  Testing  of  Fuselage  Panels  Obtained  from  Retired  Aircraft 

D.Y.  Jeong/*  SJ.  Kokkinsr*  T.H.  Flournoy,^  J.V.  Canha/ 

G.W.  Neat/  and  D.E.  Niese/ 

^  U.S.  Department  of  Transportation 
Research  and  Special  Programs  Administration 
John  A.  Volpe  National  Transportation  Systems  Center 
Cambridge,  MA  02142 

^  Foster-Miller,  Inc. 

Waltham,  MA  02154 

^  U.S.  Department  of  Transportation 
Federal  Aviation  Administration 
William  J.  Hughes  Technical  Center 
Atlantic  City,  NJ  08405 

“U.S.  Air  Force 
OC-ALC/LACRA 
Tinker  AFB,  OK  73145 


ABSTRACT.  Under  the  joint  sponsorship  of  the  U.S.  Department  of  Transpor¬ 
tation  and  the  U.S.  Air  Force,  Foster-Miller,  Inc.  is  carrying  out  a  research 
program  to  conduct  fatigue  and  residual  strength  testing  of  fuU-scale,  curved,  stiff¬ 
ened,  fuselage  sections  from  retired  Boeing  707  and  C/KC-135  aircraft.  This 
program  also  includes  non-destructive  inspection  to  detect  the  presence  of  any 
cracking  or  corrosion  damage  prior  to  panel  testing.  The  specific  objectives  of  this 
program  are:  (1)  to  demonstrate  the  feasibility  of  performing  realistic  cyclic  test¬ 
ing  of  fuselage  structures  removed  from  actual  aircraft;  and  (2)  to  perform  full- 
scale  tests  to  determine  fatigue  and  residual  strength  of  fuselage  panels  with  lap 
splices.  This  paper  will  report  on  the  status  of  the  work-in-progress. 


1.0  BACKGROUND 

In  support  of  the  Federal  Aviation  Administration  William  J.  Hughes  Technical 
Center’s  (FAAAVJHTC)  National  Aging  Aircraft  Program,  a  facility  to  test  full- 
scale  fuselage  panels  was  designed  and  built  by  Foster-MiUer,  Inc.  (FMI)  in  the  late 
1980s  [1]  under  contract  to  the  John  A.  Volpe  National  Transportation  Systems 


Presenters  at  the  1997  USAF  Aircraft  Structural  Integrity  Program  (ASIP)  Conference. 
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Center  (Volpe  Center).  This  facility,  referred  to  as  the  FAA  Aging  Aircraft  Test 
Fixture  (Figure  1),  has  been  used  to  examine  the  fatigue  and  residual  strength  of 
fuselage  panels  containing  widespread  fatigue  damage  (WFD)  in  lap  splices  [2].  A 
unique  feature  of  this  fixture  is  that  it  employs  water  rather  than  air  as  a  pressuri¬ 
zation  medium  to  create  biaxial  loading.  Whiffle- tree  systems  are  used  to 
distribute  force  uniformly  across  the  ends  of  the  test  panel.  Hydraulic  cylinders  at 
one  end  of  the  panel  provide  this  pulling  force.  A  series  of  tumbuckles  anchor  the 
panel  on  the  sides  which  incorporate  specially-designed  elastic  restraints  to  provide 
proper  loading  of  both  skins  and  circumferential  frames.  An  inflatable  seal  closes 
the  gap  between  the  panel  and  the  tub  to  prevent  leakage. 


Figure  1.  FAA  Aging  Aircraft  Test  Fixture. 

Prior  to  1997,  only  laboratory  panels  with  structural  details  resembling  a  Boeing 
737  (B-737)  airplane  had  been  used  at  the  Aging  Aircraft  Test  Facility.  Validation 
of  the  laboratory  panels  was  achieved  when  strain  data  collected  from  the  panels 
agreed  reasonably  well  with  strain  data  collected  during  a  ground-pressurization 
test  of  an  actual  B-737  airplane  [3]. 

The  U.S.  Department  of  Transportation  (U.S.  DOT)  and  the  U.S.  Air  Force 
(USAF)  are  now  collaborating  on  a  joint  research  program  to  examine  the  fatigue 
and  static  strength  of  fuselage  panels  obtained  from  disassembly  of  retired  aircraft. 
Specifically,  panels  have  been  obtained  from  retired  Boeing  707  (B-707)  and 
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C/KC-135  aircraft.  In  other  words,  panels  from  actual  aircraft  are  being  used  in 
this  on-going  program  which  comprises  two  parts.  The  first  part  involves  non¬ 
destructive  inspection  (NDI)  characterization  of  the  panels  to  detect  any  type  of 
damage  (i.e.,  fatigue  cracking  and/or  corrosion)  prior  to  testing.  The  second  part 
entails  a  testing  program  to  evaluate  the  fatigue  and  residual  strength  of  the  panels. 

This  paper  will  describe  the  work  in-progress  on  this  research  program.  To  date, 
NDI  evaluations  have  been  completed  on  four  aircraft  panels  (two  B-707  panels 
and  two  C/KC-135  panels),  and  testing  has  been  completed  on  two  B-707  panels. 


2.0  DESCMPTION  AND  USAGE  HISTORY  OF  AIRCRAFT  PANELS 

The  FAA  Aging  Aircraft  Test  Facility  accommodates  panels  that  are  120  inches  in 
length  and  68  inches  along  the  circumference  with  a  constant  radius  of  curvature 
of  75  inches.  This  radius  may  vary  over  a  limited  range,  but  due  to  the  modular 
design  of  the  fixture,  substitution  of  certain  components  can  be  made  to  accom¬ 
modate  a  wide  range  of  radii  from  commuter  to  wide-body  types  of  aircraft.  In  the 
current  test  program,  panels  with  limited  structural  interaction  with  windows, 
doors,  and  other  design  details  were  preferred.  Because  of  these  restrictions,  fu¬ 
selage  panels  were  taken  from  the  crown  area  of  various  retired  airplanes.  Table  1 
gives  a  brief  description  of  the  fuselage  panels  collected  for  the  present  test  pro¬ 
gram  to  date.  Structural  variations  among  the  panels  are  noted  in  the  table  as  well. 


Table  1.  Summary  of  Panel  Descriptions. 


Panel 

No. 

USAP 

ID  No. 

Type 

A/C 

Stringer 

Pitch 

Body 

Stations 

Stringer-Skin 

Connection 

Skin 

Thickness 

Remarks 

1 

C2186 

B-707 

9.5  in. 

BS540  to 
BS600 
Port  Side 

Spot  welds 
and  rivets  at 
frames 

0.046  in. 
to 

0.047  in. 

Tear  strap  width  is  3 
in.  Stringer-frame  ties 
at  all  stringers. 

2 

CZ186 

B-707 

9.5  in. 

BS1080  to 
BS1240 
Port  Side 

Continuous 
rivets  only 

0.062  in. 
to 

0.064  in. 

Tear  strap  width  is  3 
in.  Less  than  half  of 
stringer-frames 
crossings  had  ties, 
but  were  removed 
before  testing. 

3 

CA029 

C-135 

9.0  in. 

BS980  to 
BS1140 
Port  Side 

Spot  welds 
and  rivets  at 
frames 

0.068  in. 
to 

0.070  in. 

No  tear  straps.  No 
stringer-frame  ties. 

4 

CA029 

C-135 

9.0  in. 

BS440  to 
BS680 
Port  Side 

Spot  welds 
and  rivets  at 
frames 

0.048  in. 
to 

0.049  in. 

No  tear  straps.  No 
stringer-frame  ties. 

NOTE:  The  fuselage  radius  is  74  inches  for  both  the  B-707  and  C/KC-135  aircraft. 
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Two  panels  were  obtained  from  a  B-707  airplane  that  had  accumulated  22,071 
flight  cycles  during  77,742  flight  hours  of  in-service  usage  between  January  1967 
and  October  1990.  Although  both  panels  were  taken  from  the  same  airplane,  the 
skin  thicknesses  were  nominally  different.  A  panel  taken  from  the  crown  area  for¬ 
ward  of  the  wing  between  Body  Stations  BS540  and  BS600  had  a  skin  thickness 
of  0.046  to  0.047  inch,  including  coatings.  A  panel  taken  from  the  crown  area  aft 
of  the  wing  between  Body  Stations  BS1080  and  BS1240  had  a  skin  thickness  be¬ 
tween  0.062  and  0.064  inch.  Other  structural  differences  were  observed  in  these 
two  fuselage  sections  in  addition  to  the  difference  in  skin  thickness.  For  example, 
in  the  forward  section  of  the  fuselage,  the  stringers  were  attached  to  the  skins  via 
spot  welds.  In  the  aft  section,  the  stringers  were  attached  to  the  skin  by  riveting. 
Crossings  of  the  frames  and  the  stringers  were  reinforced  with  tie  clips  in  the  for¬ 
ward-section  panel.  In  the  aft  section,  ties  were  located  at  less  than  half  of  the 
stringer-frames  crossings  when  the  panel  was  delivered  to  the  test  facility.  Before 
testing,  however,  the  stringer-frame  ties  were  removed  from  the  aft-section  panel 
to  ensure  uniform  load  transfer.  Also,  variations  in  stringer  cross-section  are  evi¬ 
dent,  even  within  the  same  runs  in  a  given  panel. 

Two  panels  were  obtained  from  a  C/KC-135  airplane  that  accumulated  2,792 
flight  cycles  during  14,267  flight  hours  between  February  1958  and  July  1992.  As 
in  the  case  of  the  B-707  airplane,  the  two  panels  from  the  same  C/KC-135  airplane 
had  different  skin  thicknesses.  A  panel  t^en  forward  of  the  wing  between  Body 
Stations  BS440  and  BS680  had  a  nominal  skin  thickness  of  0.048  and  0.049  inch. 
A  panel  taken  aft  of  the  wing  between  Body  Stations  BS940  and  BS1140  had  a 
nominal  skin  thickness  of  0.068  and  0.070  inch.  Riveting  of  the  lap  joints  in  the 
C/KC-135  fuselage  panels  did  not  employ  the  commercial- type  machined  counter¬ 
sink  on  the  skin  (as  used  in  the  B-707  lap  joints).  Rather,  after  positioning  and 
drilling  of  the  rivet  holes,  the  two  skins  were  dimpled  together  in  a  countersink  die, 
producing  a  cone-shaped  engagement  at  every  rivet  hole,  after  which  a  flush-head 
rivet  was  installed. 


3.0  PANEL  PREPARATION 

The  aircraft  sections  removed  from  the  retired  airplanes  required  some  modifica¬ 
tions  to  accommodate  the  test  fixture.  These  modifications  included: 

•  sizing  the  test  area  of  the  panel  to  ensure  a  longitudinal  lap  joint  near  the  center 

•  cleaning  and  degreasing  the  panels 

•  stripping  paint  at  the  lap  joints  and  edges  for  instrumentation 

•  removing  stringers  at  locations  where  the  panel  is  attached  to  the  fixture 
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•  cutting  skin  edges  and  frame  ends  to  final  size 

•  installing  doubler  strips  on  each  side  of  the  panel  (three  layers  per  side) 

•  removing  of  stringer-frame  ties  (in  the  second  B-707  panel  only) 

•  reinforcing  the  frame  ends 

•  drilling  loading  holes 

During  the  panel  preparations  it  was  discovered  that  the  longitudinal  lap  splices 
were  not  bonded  with  an  adhesive.  However,  it  was  found  that  a  zinc  chromate 
primer  was  applied  to  both  faying  surfaces,  and  was  probably  wet  when  the  panel 
skins  were  joined.  This  primer  dried  in  place,  and  in  many,  but  not  all,  places  con¬ 
nected  the  upper  and  lower  skins  in  the  overlap. 


4.0  NON-DESTRUCTIVE  INSPECTION  (NDI)  OF  FUSELAGE  PANELS 

The  B-707  and  C/KC-135  panels  were  inspected  with  five  different  NDI  tech¬ 
niques  to  detect  for  prior  damage:^ 

•  Thermal  Wave  Imaging  (W ayne  State  University) 

•  Magneto-Optic  Imaging  (PRI  Instrumentation) 

•  D-Sight  (National  Research  Council  Canada) 

•  Pulsed  Eddy  Current  (Iowa  State  University) 

•  Ultrasound  via  “dripless  bubbler”  technique^  (Iowa  State  University) 

A  description  of  these  various  NDI  techniques  is  beyond  the  scope  of  this  paper, 
but  specific  details  of  these  methods  can  be  found  in  the  open  literature  [4-8]. 

No  evidence  of  fatigue  cracking  was  found  in  any  of  the  fuselage  panels.  In 
most  areas,  the  NDI  surveys  found  little,  if  any  corrosion  in  these  panels  (which,  as 
noted  previously,  were  taken  from  aircraft  over  30  years  old).  However,  in  Panel 
#2  (the  thicker  of  the  two  B-707  panels  which  was  taken  from  an  aft  crown  sec¬ 
tion),  a  few  areas  of  hidden  lap  joint  corrosion  in  the  range  of  10  to  20%  thickness 
loss  were  suggested  by  two  NDI  methods:  pulsed  eddy  current  and  ultrasound. 
The  effects  of  this  will  be  determined  during  the  testing  of  the  panel.  Also,  in  all 
cases,  post-test  dismantling  of  the  joints  will  provide  direct  viewing  of  these  inter¬ 
nal  joint  surfaces. 


'■  Affiliations  of  the  technical  staff  perfonning  the  various  NDI  techniques  are  listed  in  paren¬ 
theses. 

■  “Dripless  bubbler”  essentially  refers  to  a  captured  water  column  that  incorporates  a  water  vac¬ 
uum  return  system. 
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5.0  CORRELATION  OF  STRAIN  GAGE  DATA  WITH  ANALYSIS 


Strain  gage  measurements  were  taken  at  various  locations  on  a  B-707  panel  to  en¬ 
sure  proper  loading.  Figure  2  shows  a  schematic  of  the  strain  gage  layout  for  the 
instrumented  B-707  panel.  The  panel  was  instrumented  with  ten  gages,  located 
mostly  along  the  vertical  centerline  (midway)  between  frames.  All  but  two  gages 
were  oriented  to  measure  strains  in  the  circumferential  (hoop)  direction.  The  other 
two  gages  measured  strains  in  the  longitudinal  (axial)  direction. 

Finite  element  (FE)  analyses  based  on  large  deformation  theory  were  conducted 
by  the  Volpe  Center  as  part  of  the  present  research  program.  The  development  of 
the  finite  element  model  for  the  B-707  fuselage  panel  was  facilitated  based  upon 
previous  experience  in  analyzing  strain  fields  in  a  B-737  fuselage  [3].  Moreover, 
the  FE  model  developed  for  the  B-737  fuselage  panel  was  modified  to  examine  the 
B-707  panels  in  the  present  study. 
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Figure  2.  Schematic  of  strain  gage  layout  for  B-707  test  panel. 

The  actual  strains  measured  from  the  B-707  panel  at  an  internal  pressure  of  9.5 
psi  are  compared  with  the  strains  calculated  from  the  finite  element  model  at  the 
same  pressure  level  in  Table  2.  The  strain  level  in  the  lap  joint  area  could  not  be 
calculated  accurately  with  the  present  model  without  mesh  refinement.  Otherwise, 
the  calculated  strains  are  within  20%  of  the  actual  strains. 

Table  3  lists  hoop  and  longitudinal  stresses  calculated  at  the  gage  locations  by 
the  FE  model.  From  elementary  strength  of  materials  considerations  of  a  pressur¬ 
ized  thin-waUed  cylinder,  the  theoretical  ratio  between  the  longitudinal  stress  and 
the  hoop  stress  is  equal  to  one-half.  In  Table  3,  the  ratio  is  shghtly  less  than  the 


736 


theoretical  value,  but  the  FE  results  are  reasonable  because  the  fuselage  structure 
is  stiffened.  Also,  in  the  previous  study  for  the  B-737  aircraft  [3],  the  FE  analysis 
calculated  longitudinal-to-hoop  stress  ratios  between  0.475  and  0.503.  The  theo¬ 
retical  hoop  stress  for  a  thin-waUed  cylinder  with  radius  (74  inches)  and  skin 
thickness  (0.040  inch)  equal  to  that  of  a  B-707  pressurized  to  9.5  psi  is  equal  to 
17.6  ksi.  The  FE  results  for  hoop  and  longitudinal  stresses  at  the  midbay  locations 
in  the  B-707  panel  are  about  77%  and  72%  of  the  thin-walled  cylinder  approxima¬ 
tion,  respectively. 


Table  2.  Measured  and  Calculated  Strains  in  B-707  Panel  at  9.5 psi. 


Gage 

I  Strains  1 

No. 

Actual 

FE 

%  diff. 

Kb) 

157 

224 

-42.7% 

15 

2 

1404 

1119 

20.3% 

15.5 

3 

956 

1102 

15.3% 

2.5 

4 

1533 

(0 

- 

0 

5 

1161 

1082 

6.8% 

2 

6 

1281 

1113 

13.1% 

4 

7 

1297 

1114 

14.1% 

6 

8 

1176 

1118 

4.8% 

13.5 

9(b) 

195 

224 

-14.9% 

14 

10(d) 

893 

940 

-5.3% 

13.5 

NOTES: 

(a)  All  gage  locations  are  in  inches  relative  to  gage  number  4  which  was  placed  1 .5  inches  from  the  bottom  edge  of 
the  lap  splice.  All  but  one  gage  is  mounted  midway  between  frames. 

(b)  This  gage  measures  strain  in  the  longitudinal  (or  axial)  direction.  All  but  two  gages  are  oriented  to  measure 
hoop  (or  circumferential)  strains. 

(c)  The  present  finite  element  model  requires  mesh  refinement  to  calculate  accurate  strains  in  the  lap  joint  area. 

(d)  This  gage  is  mounted  above  a  frame.  All  other  gages  are  mounted  midway  between  frames. 


Table  3.  Calculated  Stresses  in  B-707  Panel  at  9.5  psi. 


Gage 

1  FE  Calculations  for  Stresses  I 

No.  (a) 

Gz  (ksi) 

06  (ksi) 

0/06 

1 

6.46 

13.7 

0.472 

2 

6.47 

13.7 

0.473 

3 

6.50 

13.6 

0.481 

4 

(b) 

- 

- 

5 

6.24 

13.2 

0.472 

6 

6.40 

13.6 

0.470 

7 

6.43 

13.6 

0.472 

8 

6.49 

13.7 

0.472 

9 

6.45 

13.7 

0.472 

10 

4.97 

11.4 

0.437 

NOTES: 

(a)  Refer  to  Table  2  for  specific  location  of  strain  gages. 

(b)  The  present  finite  element  model  requires  mesh  refinement  to  calculate  accurate  strains  in  the  lap  joint  area. 
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6.0  FATIGUE  AND  RESIDUAL  STRENGTH  TESTING  OF  8-707  PANELS 


To  date,  testing  has  been  completed  on  two  B-707  panels.  During  the  fatigue  tests 
on  both  panels,  the  pressure  was  cycled  between  1  and  9.5  psi^  at  a  rate  of  720  cy¬ 
cles  per  hour.  Fatigue  crack  growth  was  monitored  with  a  lOx  microscope  at 
roughly  1,000-cycle  intervals  after  initial  cracking  was  observed.  This  visual 
method  had  been  performed  in  previous  fatigue  tests  conducted  by  FMI.  Moreo¬ 
ver,  MSD-type  cracks  emanating  0.050  inch  from  the  rivet  head  have  been 
successfully  detected  using  this  method  in  laboratory  panels  which  were  made  from 
bare,  smooth  aluminum.  In  the  actual  B-707  aircraft,  multiple  layers  of  paint  made 
visual  detection  of  initial  cracking  difficult'^.  Visual  detection  of  cracking  was  en¬ 
hanced  when  the  panel  was  slightly  pressurized  which  opened  the  cracks. 
Furthermore,  the  panel  pressurization  was  considered  safe  since  water  was  used 
for  the  internal  loading.  Initially,  the  panel  paint  surface  had  been  left  intact  on  the 
panel,  but  after  the  first  visible  sign  of  cracking,  the  surface  of  the  panel  was 
stripped  and  cleaned  to  increase  the  likelihood  of  finding  additional  cracks. 

6.1  Panel  #1 

The  first  B-707  panel  was  taken  from  the  crown  area,  forward  of  the  wing  be¬ 
tween  Body  Stations  BS540  and  BS600E. 

6.1.1  Fatigue  Test 

The  significant  events  in  the  fatigue  test  conducted  on  the  first  B-707  panel  are 
summarized  as  follows: 

•  The  first  visible  sign  of  fatigue  cracking  was  observed  after  36,000  cycles  m 
the  fixture^ 

•  First  linkup  of  adjacent  MSD-type®  cracks  occurred  after  47,500  cycles  in  the 
fixmre. 


^  The  use  of  water  as  a  pressurization  medium  requires  tension-to-tension  cycling.  The  corre¬ 
sponding  pressure  range  in  the  fatigue  test  was  8.5  psi  which  is  approximately  equal  to  the 
nominal  operating  pressure  of  8.6  psi  for  the  B707  aircraft. 

Eddy  current-based  NDI  methods  can  be  used  for  crack  detection  if  the  exact  point  of  crack 
formation  is  desired  during  such  testing. 

^  The  cycle  counts  reported  in  this  section  are  the  number  of  cycles  accumulated  in  the  test  fix¬ 
ture.  A  representative  value  for  the  fatigue  life  of  the  panel  can  be  obtained  by  adding  the 
number  of  cycles  accumulated  during  in-service  usage  to  the  number  of  cycles  accumulated  in 
the  fixture.  The  B-707  panels  accumulated  22,071  cycles  during  in-service  usage. 
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•  All  cracking  occurred  in  the  upper  row  of  rivets  in  a  bay  characterized  as  hav¬ 
ing  “light”  corrosion  by  D-sight. 

•  Fatigue  cycling  was  stopped  after  48,616  cycles  in  the  fixture,  at  which  time 
successive  hnkup  of  MSD  cracks^  had  created  a  single  isolated  crack  with  an 
overall  length  of  14.5  inches. 

Figure  3  shows  the  crack  growth  data  obtained  from  the  first  B-707  panel  test. 
The  visible  crack  length  is  the  distance  measured  from  the  edge  of  the  rivet  hole  to 
the  crack  tip.  The  number  of  cycles  are  those  accumulated  while  the  panel  was  in 
the  test  fixture.  Rivets  were  numbered  consecutively  across  the  panel.  The  letters 
“L”  and  “R”  refer  to  the  left  and  right  sides  of  the  rivets.  The  figure  indicates  that 
initial  cracking  occurred  at  Rivets  14  and  15  which  are  roughly  midbay  between 
frames  (rivet  spacing  in  this  row  was  1  inch  with  frames  at  20-inch  intervals). 


Visible  Crack  Length  (inch) 


Figure  3.  Crack  growth  in  upper  rivet  row  of  B-707 panel. 


®  MSD  is  an  acronym  for  multiple-site  damage,  and  is  a  source  of  widespread  fatigue  damage 
(WFD)  characterized  by  the  simultaneous  presence  of  fatigue  cracks  in  the  same  structural 
element. 

’  It  was  possible  for  cycling  to  take  place  even  with  relative  long  cracks  because  the  panel  was 
sealed  internally  with  a  soft  rubber  dam  to  minimize  leakage  without  artificially  restraining 
the  edges  of  the  skin. 
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6.1.2  Residual  Strength  Test 


A  14.5-inch  lead  crack  with  MSD-type  cracks  was  created  as  a  result  of  the  fatigue 
testing  (Figure  4).  When  the  panel  was  pressurized  to  0.5  psi,  the  crack  surfaces 
were  separated  by  about  0.25  inches  due  to  bulging  in  the  lap  skins.  With  a  starting 
pressure  of  0.5  psi,  the  pressure  was  slowly  increased  at  a  rate  of  0.2  psi  per  sec¬ 
ond  until  panel  failure*  occurred  at  13.0  psi.  At  the  failure  pressure,  the  crack  tip 
closest  to  the  aft  end  turned  along  the  teai‘  strap.  The  other  crack  tip  apparently 
began  to  curve  along  the  tear  strap.  Moreover,  the  failure  resulting  from  this  test 
would  have  been  considered  a  “controlled”  depressurization  in  an  actual  aircraft. 


Figure  4.  One-bay  crack  created  by  successive  linkups  ofMSD  cracks. 

6.2  Panel  #2 

The  second  B-707  panel  (Panel  #2)  was  taken  from  the  same  airplane  as  Panel  #1. 
Panel  #2  was  cut  out  of  a  section  aft  of  the  wing,  and  had  a  greater  skin  thickness 
than  Panel  #1.  The  structural  differences  between  the  two  panels  were  described 
previously  in  Section  2.0. 

6.2.1  Fatigue  Test 

Fatigue  testing  is  now  continuing  past  100,000  cycles  with  no  visible  signs  of  fa¬ 
tigue  damage. 


®  In  this  context,  “failure”  means  that  the  damage  (i.e.,  aacking)  in  the  panel  was  such  that 
hydraulic  pressure  could  not  be  maintained  to  continue  testing. 
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6  2.2  Residual  Strength  Test 


Depending  on  the  results  of  the  fatigue  testing,  a  residual  strength  test  may  be 
conducted  on  this  thicker  panel. 


7.0  ESTIMATION  OF  WFD  THRESHOLD  FOR  B-707  PANELS 

The  onset  of  widespread  fatigue  damage  (WFD)  was  estimated  for  the  B-707  air¬ 
craft  by  applying  a  methodology  developed  during  previous  research  in  the  Aging 
Aircraft  Research  Program  [9],  This  methodology  requires  an  analysis  tool  to  de¬ 
termine  the  residual  strength  of  the  fuselage  structure  with  WFD  and  fatigue  crack 
growth  data  from  tests. 

In  the  present  study,  the  residual  strength  of  B-707  panels  was  examined  using  a 
compatibility  displacement  (CD)  analysis.  Plastic  coUapse  of  Hgaments  between 
adjacent  cracks  was  assumed  as  a  failure  criterion  for  linkup  of  MSD-type  cracks. 
The  CD  analysis,  however,  relies  on  the  knowledge  of  a  bulging  factor  to  predict 
the  residual  strength  in  a  curved  stiffened  panel.  A  bulging  factor  for  a  one-bay 
crack  in  a  B-707  panel  was  back-calculated  from  the  results  of  the  residual 
strength  test.  The  bulging  factor  for  the  one-bay  crack  was  then  modified  to  ac¬ 
count  for  a  two-bay  crack  in  the  WFD  threshold  estimates.  The  CD  analysis  also 
relies  on  the  knowledge  of  rivet  flexibility.  Two  values  for  the  rivet  flexibUity  were 
assumed  in  the  CD  analysis  to  estimate  upper  and  lower  bounds  of  residual 
strength.  The  lower  bound  estimate  represents  a  panel  without  adhesive  bonding. 
The  upper  bound  estimate  represents  a  panel  with  infinitely  stiff  rivets.^  Based  on 
the  CD  analysis,  the  corresponding  critical  MSD-crack  lengths  were  found  to  be 
0.09  and  0.15  inch. 

The  most  rapid  crack  growth  rate  observed  in  the  fatigue  test  was  used  to  esti¬ 
mate  the  number  of  cycles  for  an  MSD-type  crack  to  grow  to  the  critical  length. 
Based  on  the  data  shown  in  Figure  3,  the  MSD-crack  lengths  of  0.09  and  0.15  inch 
correspond  to  31,400  and  37,600  cycles  in  the  test  fixture,  respectively.  Adjusting 
these  values  to  account  for  the  difference  in  the  stress  ratio  ,  the  number  of  cycles 
to  reach  0.09  and  0.15  inch  are  35,100  cycles  and  42,000  cycles.  Adding  the  num¬ 
ber  of  cycles  from  in-service  usage,  the  WFD  threshold  is  estimated  to  be  between 

57,800  and  64,100  cycles. 


^  The  test  fixture  applies  stress  cycles  with  non-zero  minimum  stress.  The  ratio  of  minimum 
suess  to  maximum  stress  is  roughly  equal  to  0.1.  In  an  actual  airaaft,  the  minimum  stress, 
and  therefore  the  stress  ratio,  is  zero. 
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h  R  737  the  WFD  thresh- 
When  this  methodology  was  cycles'[9].’  These  estimates 

old  was  estimated  to  be  between  3^00  “  debonded  lap-sphce 

however,  were  ba^d  »P“"  Xc  estimates  of  the  WFD  ^ 

panel  which  may  have  '““'Xwee^te  B-707  and  the  B-737  ancralt  is  the  skin 
The  main  structural  difference  o  036  versus  0.040  inch), 

thickness  (thinner  thickness  m  the  B-737,  o.u. 


8.0  CONCLCSIONS  AND  FUTURE  WORK 


8  41  (  ^  _ 

Although  the  ioint  FAAmSAF  research  program  is  on-going,  the  foUowmg  co  - 

elusions  can  be  stated.  c,ircraft 

•  lAf  rwn  fuselage  panels  from  actual  aircratt 

.  Based  on  testing  and  mlys^t  the  ons^  ^hiTkS^estimLd  to  be  between 
in  a  B-707  panel  with  0.04U  men 
58  000  and  64,000  cycles. 


Fntureworkinthisprogramwillincludethefollowing: 


future  worK  m  uno  o — 

Patigueandres.dualstmngthte.w..^^^^^ 

pIS'b  S  lap  joint  “X^tte 

■«  in  terms  of  fatigue  and  static  ^ «  „2t  bonding  because  the  nv- 
stronger  than  the  conventionahy-nv  J  without  any  relative  mouon 

ets  and  the  engaged  skins  should  ^  ^  the  dimpled  joint  would  be 

between  them.  It  is  expected  *a' s  ?  ^d  weU-bonded  joint.  How- 
parable  to  that  of  a  oonven^Xot  B  expected  to  be  different  due 
ever,  the  ‘XPe’Prom  me  mLhined  countersink  and  the  presence 

PePll^tPeholefromdredimpie. 

.  The  lap  splices  in  the  test  panels  will  be  disassembled  after  residual  stren.t 

PPVverlfytheNDlmeasuremen.^^^^^^^^^^^^^^^^ 

JfadSt  SPl  SPned^o  the  research 
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CODSTRAN  SIMULATION  CYCLE 
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PLY  MICROSTRESSES  THROUGH  COMPOSITE  STRESS  PROGRESSIVE  DECOMPOSITION 
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Ply  (Lamina)  Fracture  Modes  Tracked  by 
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Intra  and  Inter  Ply 


Interply  Layer  (Matrix)  Fracture  Modes 
Tracked  by  CODSTRAN 
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Interply  Relative  Rotation  (A0  ) 
Scissoring  Effect 
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STRESS-STRAIN  RELATIONS  AND  DAMAGE 
PROGRESSION  FOR  WOVEN  AND  NON- WOVEN 
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composites. 

Increasing  the  shear  component  of  loading  increased  the  dif¬ 
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sitive  to  the  presence  and  magnitude  of  in-plane  shear 
stresses. 
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Abstract:  An  efficient  and  accurate  computational  methodology  and  associated  software  (PC-Rep)  for  the 
analysis  and  design  of  bonded  composite  repairs  is  described.  PC-Rep  runs  on  personal  computers 
including  laptops  and  features  a  convenient  graphical  user  interface.  The  software  can  be  used  to  compute 
a  variety  of  parameters  important  in  repair  design  such  as  patched  stress  intensity  factors  and  crack  growth 
life.  These  parameters  can  be  used  to  determine  residual  strength,  inspection  intervals  and  remaining 
structural  life.  Comparisons  with  standard  reference  solutions  are  carried  out  to  validate  the  software. 
Parametric  studies  typical  of  repair  design  are  presented. 

1.  INTRODUCTION 

Infrastructure  the  world  over  is  aging.  This  infrastructure  includes  buildings,  bridges, 
roads,  pipelines,  and  aircraft.  As  it  ages  there  is  growing  concern  about  the  structural 
integrity  of  much  of  this  infrastructure.  Economic  and  market  forces  have  forced 
industry  and  government  agencies  to  utilize  many  of  the  structures  which  make  up  the 
infrastructure  well  beyond  their  original  design  lives.  For  example,  in  1988,  22%  of 
Boeing  727's  were  operating  beyond  their  original  design  life.  This  percentage  has 
continued  to  increase  since  that  time. 

There  are  essentially  two  ways  of  dealing  with  the  aging  infrastructure  problem.  The 
first  is  to  retire  existing  infrastructure  and  replace  it.  The  second  is  to  repair  the  existing 
infrastructure.  While  the  first  solution  is  perhaps  more  desirable,  it  is  not  always  viable 
from  an  economic  perspective.  Thus,  it  is  imperative  that  strategies  to  repair  aging 
infrastructure  be  developed.  By  doing  so,  the  service  lives  of  the  infrastructure  can  be 
extended  beyond  that  for  which  it  was  originally  designed. 

A  promising  infrastructure  repair  technique  is  the  use  of  adhesively  bonded  composite 
repairs.  Bonded  composite  repairs  offer  an  economical  means  of  repairing  cracked  or 
damaged  structure  which  does  not  drastically  increase  the  weight  of  the  structure. 
Adhesively  bonded  composite  repairs  have  many  advantages  over  mechanically  fastened 
repairs  such  as:  (i)  new  stress  concentrations  are  not  introduced  into  the  repaired 
structure  due  to  rivet  holes;  (ii)  the  repair  patches  are  readily  formed  into  complex 
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shapes;  (iii)  high  stiffness  to  weight  and  strength  to  weight  ratios  of  the  patch;  and  (iv) 
high  corrosion  and  fatigue  resistance  of  the  composite.  This  repair  technique  has  been 
primarily  used  in  the  area  of  military  aviation.  However,  it  is  equally  well  suited  for  use 
in  commercial  aviation  as  well  as  the  repair  of  civil  infrastructure  such  as  buildings  and 
bridges. 

In  order  to  properly  design  an  adhesively  bonded  composite  repair  of  a  structure,  many 
factors  must  be  considered.  These  include: 

•  global  stiffening  of  the  structure  due  to  the  high  stiffness  of  the  composite  patch; 

•  the  effect  of  size,  shape,  thickness  (including  taper)  and  material  properties  of  the 
composite  patch  on  the  crack-tip  stress  intensity  factors,  skin  and  patch  stresses, 
adhesive  shear  strain,  and  peel  stresses; 

•  the  effect  of  the  material  properties  of  the  adhesive  on  the  crack-tip  stress 
intensity  factors,  skin  and  patch  stresses,  adhesive  shear  strain,  and  peel  stresses; 

•  the  effect  of  thermal  cycling  on  the  composite  repair;  and 

•  the  effect  of  disbonds  on  the  performance  of  the  composite  repair. 

This  paper  describes  an  efficient  and  accurate  computational  methodology  and  associated 
software  (PC-Rep)  for  the  analysis  and  design  of  bonded  composite  repairs.  It  builds  on 
the  work  of  Park,  Ogiso  and  Atluri  (1992);  and  Nagaswamy,  Pipkins  and  Atluri  (1996). 
PC-Rep  runs  on  personal  computers  including  laptops  and  features  a  convenient 
graphical  user  interface.  The  software  can  be  used  to  compute  a  variety  of  parameters 
important  in  repair  design  such  as  patched  stress  intensity  factors  and  crack  growth  life. 
These  parameters  can  be  used  to  determine  residual  strength,  inspection  intervals  and 
remaining  structural  life. 

2.  ANALYSIS  OF  REPAIRS 

There  are  several  methods  for  the  analysis  of  patched  cracks.  They  can  be  broadly 
classified  as  either  analytical  or  numerical.  A  review  of  these  methods  can  be  found  in 
Nagaswamy,  Pipkins  and  Atluri  (1996).  In  this  paper  the  focus  will  be  on  the  usability 
of  these  methods  in  a  maintenance  environment.  The  repair  design  requirements  in  this 
environment  are  conflicting;  fast  turnaround  (minutes)  yet  accurate.  From  the  structural 
design  point  of  view,  accuracy  usually  means  utilizing  a  numerical  method  such  as  the 
finite  element  method.  However,  finite  element  analyses  (FEA)  may  not  always  be  done 
quickly  in  a  maintenance  environment  for  several  reasons.  These  include:  lack  of 
personnel  knowledgeable  in  FEA;  the  unavailability  of  FEA  software  (pre  and  post 
processors  and  solver)  due  to  cost;  and  the  unavailability  of  computer  hardware  needed 
for  FEA.  As  a  result  of  one  or  more  of  these  reasons,  traditional  FEA  software  is  not 
typically  of  great  interest  to  maintenance  organizations. 

In  recognition  of  this  fact,  several  pieces  of  software  featuring  easy  to  use  graphical  user 
interfaces  (GUI),  fast  solution  times,  and  minimal  hardware  requirements  have  come  into 
use  in  recent  years.  These  include  CalcuRep,  RAPID  and  PC-Rep.  CalcuRep  and  PC- 
Rep  are  applicable  to  bonded  composite  repairs  while  RAPID  is  for  mechanical  doubler 
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type  repairs.  CalcuRep  [Fredell  (1994)]  is  based  on  the  analytical  approach  of  Rose 
(1981).  RAPID,  developed  by  the  Federal  Aviation  Administration  (FAA),  is  also  based 
on  an  analytical  approach  wherein  displacement  compatibility  is  enforced  between  the 
repair,  structure  and  mechanical  fasteners.  PC-Rep,  developed  by  Knowledge  Systems, 
Inc.,  has  an  outstanding  feature  which  distinguishes  it  from  these  two  codes.  The 
analysis  performed  by  PC-Rep  is  based  upon  the  finite  element  method.  However,  PC- 
Rep  relieves  the  user  of  the  tedium  often  associated  with  finite  element  analysis.  This  is 
so  because  of  the  automated  state-of-the-art  meshing  and  computational  fracture 
mechanics  algorithms  that  have  been  implemented  in  PC-Rep.  Thus,  the  user  can  enjoy 
the  benefits  of  finite  elements  (i.e.  arbitrary  geometries  and  unsurpassed  accuracy,  etc.) 
without  having  to  worry  about  any  details  of  finite  element  analysis. 

The  analysis  carried  out  by  PC-Rep  is  based  on  the  Finite  Element  Alternating  Method 
(FEAM)  [Atluri(1997)].  Because  the  repaired  structure  is  not  homogeneous,  the  FEAM 
must  be  used  in  a  two  stage  analysis  [Nagaswamy,  Pipkins  and  Atluri  (1996)].  The 
stages  are  as  follows. 

1 .  Evaluation  of  the  stresses  exerted  by  the  adhesive  and  patch  on  the  base  sheet, 

using  a  coarse  mesh. 

2.  The  stresses  obtained  from  stage  1  are  applied  as  body  forces  on  the  base  sheet, 

and  the  Finite  Element  Alternating  Method  is  used  to  find  the  stress  intensity  factor. 

In  stage  1,  a  traditional  finite  element  methodology  is  used  to  deduce  the  stresses  exerted 
by  the  patch.  As  the  crack  tip  is  not  meshed  for  the  singularity,  a  coarse  mesh  is 
sufficient  for  this  purpose.  The  sheet  and  the  patch  are  modeled  with  eight  noded  2D 
elements.  The  adhesive  is  modeled  by  shear  elements  as  in  Jones  and  Callinan  (1981). 

Stage  2  consists  of  an  analysis  using  the  2D  Finite  Element  Alternating  Method  [Atluri 
(1997)].  The  FEAM  is  based  upon  superposition  as  illustrated  in  Figure  2.1.  In  the 
FEAM,  the  stresses  in  the  uncracked  body  are  first  analyzed,  by  a  traditional  finite 
element  method,  for  the  given  system  of  external  loading.  To  model  the  crack,  the 
tractions  at  the  locations  of  the  crack  in  an  otherwise  uncracked  body  must  be  erased.  To 
erase  the  crack  face  stresses,  the  tractions  found  by  the  finite  element  solution  are 
reversed.  As  is  often  the  case,  analytical  solutions  exist  for  cracks  subjected  to  arbitrary 
crack  face  stresses  but  for  infinite  bodies.  In  the  present  2D  problem  a  solution  is  used  in 
which  the  crack  face  tractions  are  given  in  terms  of  Chebyschev  polynomials.  Therefore, 
if  these  solutions  are  to  be  used  for  erasing  crack  face  stresses,  the  residual  stresses  at  the 
finite  body  extent  in  the  infinite  body  have  to  be  erased.  This  is  done  by  reversing  the 
residual  stresses  applied  to  the  finite  uncracked  body.  This  results  in  residual  tractions 
on  the  crack  face.  To  erase  these,  the  analytical  solution  is  used  again.  This  is  repeated 
until  convergence  is  achieved.  Convergence  is  achieved  when  the  stress  intensity  factors 
for  each  analytical  iteration  become  small.  The  final  stress  intensity  factor  is  the  sum  of 
all  the  stress  intensity  factors. 
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Summarizing,  the  steps  in  the  2D  Finite  Element  Alternating  Method  are: 

1.  Solve  for  the  uncracked  finite  body  under  given  loads  by  traditional  finite 
element  method. 

2.  Compute  the  stresses  at  the  crack  location  from  step  1. 

3.  To  reverse  the  above  found  stresses,  use  a  least-squares  approximation  to  fit  the 
stresses  to  Chebyschev  polynomials. 

4.  Use  the  above  solution  and  find  the  stresses  at  the  boundaries  of  the  finite  body. 

5.  Sum  the  Chebyschev  coefficients  to  find  the  stress  intensity  factors  for  this 
iteration. 

6.  If  the  stress  intensity  factor  for  this  iteration  is  small,  then  convergence  is 
obtained. 

7.  If  convergence  is  not  obtained,  apply  the  reversed  residual  stresses  on  the  finite 
body  boundaries  and  use  the  traditional  finite  element  method  to  solve  for  the 
displacements  and  go  to  step  2. 

The  solution  of  the  finite  size  body  with  a  crack  is  then  obtained  by  summing  the  results 
(i.e.  displacements,  stress  intensity  factors,  etc.)  from  all  iterations. 

3.  PC-REP 

The  finite  element  based  analysis  technology  described  in  the  previous  section  has  been 
combined  with  a  graphical,  menu  driven  front  end  and  an  automatic  mesh  generator  and 
ported  onto  a  personal  computer.  This  computer  package,  referred  to  as  PC-Rep,  is 
extremely  robust  and  is  capable  of  fast  and  accurate  calculation  of  various  parameters 
relating  to  bonded  patch  repair.  A  schematic  of  a  typical  patched  crack  configuration  is 
illustrated  in  Figure  3.1 

3.1.  Features 

At  present  this  robust  package  has  many  attractive  features  including: 

•  Convenient  Graphical  User  Interface  (GUI) 

•  Physically  Based  Inputs 

•  Automated  Analysis  and  Mesh  Generation 

•  Finite  Element  Kernel 

•  Rectangular  Metal  Skins 

•  Adhesively  Bonded  Metal  or  Composite  Patches;  Centrally  Located  on  Skin 

•  One  Straight  Crack  Centrally  Located  under  the  Patch 

•  Remote  Uniform  Loading  Normal  to  Crack 

•  Partial  Disbond/No  Disbond 

•  Stress  Intensity  Factor  Calculation 

•  Fatigue  Lifetime  Determinations 

•  Adhesive  Shear  Stress  Computation 

•  Material  Property  Database 

•  Patch  Shape;  Rectangular,  Octagonal,  Elliptical 

•  Patch  Taper 
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•  Crack  Size;  Can  Extend  Beyond  Patch 

•  Stiffeners;  Broken,  Unbroken 

•  Residual  Stresses  from  Curing 

•  Thermal  Cycling  due  to  Temperamre  Differential  Between  Ground  and  Cruise 

Perhaps  the  most  significant  feature  of  PC-Rep  is  that  all  inputs  to  the  program  are  actual 
physical  quantities  as  distinct  from  finite  element  related  data.  Based  on  these  physical 
data,  such  as  patch  size,  stiffener  dimensions,  adhesive  properties,  a  pre-processor  will 
automatically  generate  the  optimum  finite  element  mesh  and  the  relevant  boundary 
conditions.  This  greatly  simplifies  the  usage  of  PC-Rep. 

Another  major  advantage  of  PC-Rep  is  that  it  allows  the  user  a  considerable  level  of 
complexity  with  respect  to  the  design  and  analysis  of  bonded  patches.  Despite  this 
complexity,  PC-Rep  executes  extremely  quickly  on  a  personal  computer;  a  matter  of 
minutes.  This  efficiency  is  achieved  through  the  use  of  the  finite  element  alternating 
method.  Some  of  these  features  are  very  important  when  it  comes  to  the  design  of 
patches  for  aircraft  components.  Broken/unbroken  stiffeners  can  interact  significantly 
with  the  crack  growth  process.  Residual  stresses  and  thermal  stresses  play  an  important 
role  in  the  growth  of  a  crack.  Patch  taper,  along  with  disbond  must  always  be 
considered.  All  of  these  features  are  captured  within  a  standard  PC-Rep  analysis.  Many 
of  these  will  be  illustrated  in  Section  4. 

3.2.  Graphical  User  Interface 

The  PC-Rep  package  consists  of  three  parts;  (a)  a  GUI,  (b)  a  finite  element  kernel  and  (c) 
an  automatic  mesh  generator.  The  GUI  is  written  in  tcl/tk  while  the  mesh  generator  and 
the  kernel  are  written  in  FORTRAN.  All  user  interaction  takes  place  through  the  GUI. 
The  user  does  not  have  to  be  concerned  with  the  mesh  generation  or  other  finite  element 
related  details. 

The  GUI  uses  a  convenient  point  and  click  system  that  allows  the  user  to  select  various 
options  such  as  patch  configuration,  stiffener  type,  etc.  and  to  specify  the  appropriate 
inputs.  Figure  3.2  illustrates  a  typical  screen  that  the  user  can  generate.  The  right  hand 
side  contains  buttons  that  can  be  used  to  activate  the  different  options.  The  left  hand  side 
of  Figure  3.2  contains  a  schematic  diagram  of  the  patch  configuration  that  is  being 
prepared  for  analysis.  This  allows  the  user  to  view  each  update  that  is  made  to  the 
analysis  model. 

When  the  user  selects  a  particular  option  from  the  buttons  on  the  right  hand  side,  a  new 
screen  appears  that  allows  the  user  to  input,  in  a  tabular  manner,  various  quantities 
relevant  to  that  option.  For  example.  Figure  3.3  illustrates  the  table  that  appears  when 
the  Patch  option  is  selected.  The  user  can  select  from  a  range  of  different  patch 
materials,  including  composites  and  metals,  or  else  can  input  the  material  properties 
directly.  The  user  can  also  select  the  patch  shape  (rectangular,  tapered,  octagonal  or 
elliptical)  along  with  the  patch  dimensions  using  this  screen. 
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Using  the  Analysis  Control  option,  the  user  can  select  the  type  of  analysis  that  is  desired. 
This  can  be  the  direct  calculation  of  the  stress  intensity  factors  for  the  patched  crack  or 
else  it  can  be  the  calculation  of  the  fatigue  life  of  the  component,  based  on  the  Forman 
fatigue  relationship.  Once  the  analysis  is  complete,  the  user  can  utilize  the  post¬ 
processing  options  to  view  the  most  significant  results.  These  include  plots  of  the  stress 
intensity  factor  or  the  fatigue  life  (number  of  load  cycles)  versus  the  crack  length.  The 
user  can  also  examine  important  quantities  such  as  the  maximum  shear  stress  in  the 
adhesive  for  different  patch  configurations. 

4.  ANALYSIS  AND  RESULTS 

The  PC-Rep  code  is  used  here  to  analyze  a  wide  variety  of  patched  aircraft 
configurations. 

4.1.  Validation  Problems 

To  validate  the  PC-Rep  analyses,  comparisons  were  made  with  three  results  from  the 
literature. 

4.1.1.  Rose's  Solution 

A  number  of  comparisons  were  made  with  the  solutions  developed  by  Rose  (1982). 
Here,  a  rectangular  patch  was  placed  over  a  rectangular  plate  and  a  crack  is  centrally 
located  underneath  the  patch.  The  relevant  geometrical  and  material  property  data  are 
given  in  Table  4.1.  In  this  table,  tj  is  the  skin  thickness,  Ej  is  the  elastic  modulus  of  the 
skin.  Vs  is  the  Poisson  ratio  of  the  skin.  Ex  is  the  stiffness  of  the  patch  parallel  to  the  crack 
(x  direction),  Ey  is  the  stiffness  of  the  patch  perpendicular  to  the  crack  (y  direction),  Gxy 
is  the  shear  modulus  of  the  patch,  Vyx  is  the  Poisson  ratio  of  the  patch,  Ga  is  the  shear 
modulus  of  the  adhesive  and  ta  is  the  thickness  of  the  adhesive. 

For  the  Rose  solution,  the  plate  and  patch  are  infinitely  wide  but  the  width  of  the  plate 
and  the  patch  are  taken  as  200  mm  in  the  PC-Rep  calculation.  A  uniform  unit  normal 
traction  (1  MPa)  is  applied  to  the  base  plate  in  the  direction  normal  to  the  crack.  There 
are  no  stiffeners  or  disbonds  included  in  this  calculation.  The  computed  stress  intensity 
factors  are  presented  in  Table  4.2  for  a  range  of  different  values  of  the  patch  thickness.  It 
can  be  seen  that  there  is  excellent  agreement  between  the  two  sets  of  results. 

While  the  results  in  Table  4.2  are  in  very  good  agreement,  it  is  important  to  point  out  a 
number  of  differences  that  exist  between  the  two  approaches. 

(a)  The  Rose  (1982)  solution  is  applicable  when  the  patched  plate  is  of  infinite  width 
but  PC-Rep  considers  the  case  of  a  finite  plate.  When  the  patch  width  is  less  than  the 
plate  width,  the  PC-Rep  results  will,  not  surprisingly,  be  higher  than  the  Rose  solution. 
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To  minimize  this  effect,  the  width  of  the  plate  is  selected  to  be  much  larger  than  the 
length  of  the  crack 


(b)  The  Rose  solution  assumes  the  flexibility  of  the  adhesive  layer  to  be 

F  =  ^  (6.1) 

However,  PC-Rep  uses  the  Jones  and  Callinan  (1981)  approach  where  the  adhesive  layer 
flexibility  is  modified  by  considering  the  effects  of  shear  deformation  of  the  patch  and 
the  base  plate.  In.this  case. 


G„  8  G,  8  G„ 


(6.2) 


This  effect  is  more  pronounced  as  the  thickness  of  the  patch  increases  and  is  also 
significant  when  the  width  of  the  patch  is  less  than  the  width  of  the  plate. 

4.1.2.  Experimental  Data 

Comparisons  were  also  made  with  the  experimental  data  of  Denney  (1995).  Table  4.3 
contains  the  appropriate  material  property  and  geometrical  data  and  in  this  instance,  the 
width  of  the  patch  is  considerably  less  than  the  width  of  the  plate  (tp  is  the  thickness  of 
the  patch).  The  maximum  and  minimum  remote  uniform  stresses  are  Cmax  and  Gmin 
respectively.  The  shear  modulus  of  the  adhesive  was  taken  as  21.5  MPa  for  the  purpose 
of  these  calculations  [Chow  and  Atluri  (1997)].  This  shear  modulus  is  consistent  with 
the  operating  stress  in  the  adhesive  rather  than  the  initial  elastic  value  of  405MPa.  The 
results,  showing  crack  length  as  a  function  of  the  number  of  load  cycles,  are  presented  in 
Figure  4. 1  and  the  excellent  agreement  between  the  experimental  and  computed  results  is 
obvious. 


4.1.3.  Other  Comparisons 

Comparisons  have  also  been  made  with  the  results  obtained  by  Park,  Ogiso  and  Atluri 
(1992).  The  configuration  analyzed  in  this  reference  is  given  in  Table  4.4  and  the 
comparison  of  computed  stress  intensity  factor  are  given  in  Table  4.5.  The  lower  results 
in  [Park,  Ogiso  and  Atluri  (1992)]  are  likely  due  to  the  fact  that  the  patch  is  of  infinite 
width. 
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4.2.  Stress  Intensity  Factor  Variation 


For  the  parameters  in  Table  4.4,  a  typical  plot  of  the  stress  intensity  factor  as  a  function 
of  crack  length  is  given  in  Figure  4.2.  The  patch  thickness  is  taken  as  0.6  mm  and  the 
applied  load  is  100  MPa.  As  can  be  seen,  the  stress  intensity  tends  to  level  off  as  the 
crack  length  approaches  the  edge  of  the  patch.  After  the  crack  extends  beyond  the  patch, 
the  stress  intensity  factor  increases  dramatically.  These  trends  are  typical  of  the  stress 
intensity  factor  behavior  that  is  observed  for  patched  cracks. 

4.3.  Effect  of  Stiffeners  on  Crack  Growth 

PC-Rep  can  also  be  used  to  consider  the  influence  of  stiffeners  on  crack  growth.  In  the 
first  example,  two  horizontal  and  two  vertical  stiffeners  are  included.  This  configuration 
is  illustrated  schematically  in  Figure  4.3  and  the  stiffeners  are  symmetrical  with  respect 
to  the  crack.  The  skin,  patch  and  adhesive  geometric  and  material  properties  are  the 
same  as  in  Table  4.3.  All  stiffeners  are  identical  and  have  a  cross  section  area  of  80 
mm  .  The  horizontal  stiffeners  are  75  mm  from  the  crack  and  the  vertical  stiffeners  are 
50  mm  from  the  crack  center.  Figure  4.4  illustrates  the  computed  fatigue  lifetimes  in  the 
presence  of  stiffeners.  It  can  be  seen  that  this  lifetime  is  much  greater  than  the  case 
where  there  are  no  stiffeners.  The  effect  of  the  vertical  stiffeners  is  to  retard  the  stress 
intensity  factor  increase  as  the  crack  approaches  the  stiffener. 

The  second  example  illustrates  the  effect  of  a  broken  stiffener.  Here  a  single  vertical 
stiffener  bisects  the  crack  and  in  one  instance  the  stiffener  is  broken  and  in  the  second  it 
is  unbroken.  These  results  are  shown  in  Figure  4.5.  As  expected,  the  computed  lifetime 
is  much  greater  when  the  stiffener  is  unbroken.  These  examples  illustrate  how  easy  it  is 
to  include  the  important  effects  of  stiffeners  in  a  bonded  patch  analysis. 

4.4.  Effects  of  Patch  Size  and  Shape 

Given  the  flexibility  of  PC-Rep,  it  is  relatively  easy  to  conduct  several  parametric  studies 
relating  to  the  size,  shape  and  properties  of  bonded  composite  patches.  A  selection  of 
these  results  now  follows.  In  each  case,  the  base  problem  is  that  described  in  Table  4.3. 
There  are  no  stiffeners  or  disbond  and  the  effect  of  thermal  stresses  has  been  ignored. 

4.4.1.  Patch  Shape 

Three  different  patch  shapes  are  considered;  rectangular,  octagonal  and  elliptical.  The 
overall  width  and  height  of  each  patch  is  the  same.  For  the  octagonal  patch,  the  comer 
cut-outs  are  all  10  mm  in  height  and  10  mm  in  width.  The  fatigue  lifetimes  are  presented 
in  Figure  4.6  and  it  can  be  seen  that  there  is  not  a  great  deal  of  difference  in  the  results. 
The  maximum  adhesive  stresses  were  also  relatively  similar  in  each  case. 
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4.4.2.  Patch  Thickness 


The  patch  thickness  is  varied  to  investigate  the  effect  on  both  the  stress  intensity  factor 
and  also  the  maximum  adhesive  shear  stress.  As  expected,  the  stress  intensity  factor 
shows  a  steady  decrease  as  the  thickness  is  increased  and  this  is  illustrated  in  Figure  4.7. 
Interestingly,  the  adhesive  stress  exhibits  a  minimum  with  the  stresses  rising  for  very 
thick  patches  and  also  very  thin  patches.  This  result  is  shown  in  Figure  4.8. 

4.4.3.  Taper 

The  effect  of  tapering  is  very  important  as  it  can  help  to  reduce  the  magnitude  of  the 
maximum  adhesive  shear  stress.  In  this  example,  the  taper  width,  W,,  and  height,  Ht,  are 
taken  as  5  mm  and  8  mm  respectively  and  the  thickness  of  the  patch  at  the  edge  of  the 
taper  is  0.15  mm.  Figure  4.9  compares  the  computed  maximum  adhesive  shear  stress 
with  and  without  taper.  The  taper  has  the  effect  of  reducing  this  maximum  stress  by  up 
to  15%.  However,  for  very  thin  patches,  the  maximum  stress  does  not  reduce 
dramatically.  This  suggests  that  the  maximum  stress  in  this  instance  is  at  the  edge  of  the 
crack.  It  was  also  observed  that  the  stress  intensity  factors  were  also  not  dramatically 
affected  by  the  taper. 

4.4.4.  Vertical  Height  of  Patch 

Figure  4.10  shows  the  variation  in  the  stress  intensity  factor  as  the  patch  height  is 
increased.  The  general  trend  is  for  the  stress  intensity  factor  to  decrease  moderately  with 
height.  The  maximum  adhesive  shear  stress  also  exhibits  a  moderate  increase  with 
increasing  patch  height. 

4.4.5.  Horizontal  Width  of  Patch 

As  expected,  the  stress  intensity  factor  increases  as  the  patch  width  decreases.  This 
result  is  illustrated  in  Figure  4.11.  The  maximum  adhesive  shear  stress  demonstrates  a 
similar  trend. 

4.5.  Residual  Stresses  and  Temperature 

As  discussed  earlier,  it  is  very  important  to  consider  the  effects  of  temperature  when 
designing  a  bonded  composite  patch  due  to  the  residual  stresses  set  up  during  the  curing 
process  and  due  the  thermal  stresses  set  up  during  the  ground-air-ground  cycle  because 
of  the  thermal  property  mismatch.  In  this  example,  the  curing  temperature  is  taken  as 
120°C,  the  room  temperature  as  20°C,  the  maximum  ground  temperature  as  50°C  and  the 
minimum  ground  temperature  as  -50°C.  All  other  parameters  are  as  in  Table  4.3  with  the 
addition  of  horizontal  and  vertical  stiffeners  that  are  75  mm  and  50  mm  from  the  crack 
center  respectively  (cross-section  area  80  mm^)  and  with  the  adhesive  stiffness  at  21.5 
MPa. 
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The  effect  of  these  thermal  stresses  is  to  increase  the  stress  intensity  factor.  This  is 
reflected  in  reduced  fatigue  lifetimes.  This  is  illustrated  in  Figure  4.12  where  four 
separate  results  are  presented; 

•  no  residual  or  thermal  stresses 

•  residual  stresses  only 

•  thermal  stresses  only 

•  residual  and  thermal  stresses 

This  result  shows  the  relative  importance  of  including  each  of  these  effects. 

5.  DISCUSSION 

This  paper  clearly  demonstrates  the  importance  of  new  generation  software  tools  that 
utilize  enhanced  automation  to  obtain  accurate  solutions  to  highly  complex  problems. 
PC-Rep  is  distinguished  from  PC  based  repair  codes  as  the  analysis  is  performed  by 
software  that  is  based  upon  the  finite  element  method.  However,  PC-Rep  relieves  the 
user  of  the  tedium  often  associated  with  finite  element  analysis.  This  is  so  because  of  the 
automated  state-of-the-art  meshing  and  computational  fracture  mechanics  algorithms  that 
have  been  implemented  in  PC-Rep.  Thus,  it  is  possible  to  enjoy  the  benefits  of  finite 
elements  (i.e.  arbitrary  geometries  and  unsurpassed  accuracy,  etc.)  without  having  to 
worry  about  any  details  of  finite  element  analysis. 

The  examples  presented  in  this  paper  illustrate  the  feasibility  of  a  personal  computer 
based  design  procedure  for  bonded  composite  repair  of  aircraft  structures.  It  was  easy  to 
consider  many  significant  features  such  as  residual  curing  stresses,  disbond,  taper  and 
stiffeners.  PC-Rep  can  be  viewed  as  the  new  generation  design/analysis  tool  that  merges 
an  easy  to  use  GUI  with  a  powerful,  yet  efficient  computational  kernel. 
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Tables 


2Ws  =  200  mm;  2Hs  =  635  mm;  ts  =  2.9  mm _ 

2Wp  =  200  nnom;  2Hp  ~  152  mm _ _ _ 

Es  =  73,000  MPa;  v,  =  0,33 _ 

Ex  =  25,400  MPa;  Ey  =  208,000  MPa;  Gxv  =  7,240  MPa;  Vy,  =  0. 168 

Ga  -  965  MPa;  tj  =0.102  mm  _ 

2a  =  38.1  mm 


Table  4.1  Geometrical  and  Material  Property  Data  used  in  Comparison  with  Rose  (1982) 


Patch 

Thickness 

(mm) 

PC-Rep  Stress  Intensity 

Factor 

(MPa  mm''^) 

Rose  Stress  Intensity 

Factor 

(MPamm''^) 

0.1 

3.04 

3.16 

0.2 

2.40 

2.46 

0.3 

2,03 

2.07 

0.4 

1.78 

1.80 

0.5 

1.59 

1.61 

0.6 

1.45 

1.45 

0.7 

1.33 

1.33 

0.8 

1.23 

1.22 

0.9 

1.14 

1.13 

1.0 

1.07 

1.06 

Table  4.2  Computed  Stress  Intensity  Factor  Variation  with  Patch  Thickness 


2Ws  =  152  mm;  2Hs  =  508  mm;  t;  =  1.0  mm _ 

2Wp  =  50  mm;  2Hp  =  56  mm;  tp  -0.381  mm _ 

Es  =  72,400  MPa;  v,  =  0.33 _ 

Ex  =  25,000  MPa;  Ey  =  210,000  MPa;  G^y  =  20,500  MPa;  v^x  =  0.168 

Ga  =  405  MPa;  ta  =  0.127  mm _ 

an,ax  =  120  MPa;  =  1 2  MPa;  R=  0.1 


Table  4.3  Geometrical  and  Material  Property  Data  used  in  Comparison  with  Denney  (1995) 

2Ws  =  152  mm;  2Hs  =  508  mm;  t;  =  2.0  mm 
2Wp  =  50  mm;  2Hp  =  80  mm 
Es  =  72,400  MPa;  v,  =  0.32 

Ex  =  26,250  MPa;  Ey  =  210,000  MPa;  Gxv  =  75,000  MPa;  Vy,  =  0.16 

Ga  =  965  MPa;  ta  =  0.102  mm _ 

g  =  1  MPa 

2a  =  20  mm  ~  ~~ 


Table  4.4  Geometrical  and  Material  Property  Data  used  in  Comparison  with  Park,  Oeiso  and  Atluri 

(1992) 
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Patch 

Thickness 

(mm) 

PC‘Rep  Stress  Intensity 

Factor 

(MPamm'^) 

Park  et  al.  Stress 

Intensity  Factor 
(MPamm'^-) 

0.1 

2.88 

2.80 

0.2 

2.39 

2.19 

0.4 

1.89 

1.68 

0.6 

1.60 

1.35 

0.8 

1.40 

1.18 

1.0 

1.25 

1.01 

Table  4.5  Computed  Stress  Intensity  Factor  Variation  with  Patch  Thickness 
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Stoiiesifer,  Brust,  and  Leis  (1993  -  Wang,  Briist,  and  Atliiri  (1997  - 

Interacting  Mixed  Mode  Series  of  Three  Papers) 
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Further  consideration  needed  to: 

-  verify  trends 

-  evaluate  the  impact  to  structural  Integrity 
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ABSTRACT 

Proper  models  of  fastened  structural  connections  must  include  all  the  significant  effects  that 
influence  their  performance.  Inherent  nonlinear  effects  make  the  numerical  treatment  of  these 
models  very  cumbersome  and  far  from  the  reach  of  non-specialists  in  finite  element  analysis.  To 
overcome  this  difficulty,  we  have  developed  a  unique  capability  within  our  p-version  FEA 
software  Stress  Check^  which  provides  for  the  computation  of  the  structural  and  the  strength 
responses  of  fastened  structural  connections.  The  implementation  accounts  for  the  following: 

•  Partial  contact  between  fastener  and  plate 

•  Fastener  shear  stiffness 

•  Material  and  geometric  nonlinearities 

•  Interference  fitting 

•  Initial  clearance  between  fastener  and  plates 

It  is  assumed  that  the  effects  of  bending  are  negligible  and  that  there  is  no  friction  between  the 
contact  surfaces.  In  other  words,  all  dominant  effects  in  shear  connections  are  accounted  for  in  a 
single  analysis  tool. 

A  unique  library  feature  of  Stress  Check  allows  for  the  creation  of  handbook-like  problem  defini¬ 
tions  for  fastened  structural  connections  parameterized  by  topological  description,  material  prop¬ 
erties  and  loadings.  The  output  data,  such  as  forces  acting  on  fasteners,  stress  distributions, 
estimates  of  the  size  of  plastic  zones,  etc.  can  be  produced  automatically  in  tabular  or  graphical 
formats.  Combined  with  Stress  Check’s  easy  to  use  handbook  framework,  this  capability  provides 
a  unique  tool  for  analyzing  structural  joints  in  a  reliable  and  effective  way. 

Solutions  of  high  quality  and  reliability  can  be  produced  in  a  reasonably  short  time  by  persons 
who  need  not  have  been  trained  in  the  use  of  finite  elements.  This  is  because  advanced  nonlinear 
finite  element  analysis  procedures  are  available  from  the  simple  handbook-like  interface.  This  fea¬ 
ture  delivers  important  benefits:  The  ability  to  standardize  the  analysis  of  fastened  structural  con¬ 
nections  on  the  basis  of  the  most  advanced  FEA  procedures  available  today,  while  substantially 
reducing  time  and  costs  of  design  and  analysis  computations.  This  feature  also  allows  the  compu¬ 
tation  of  stress  intensity  factors  for  cracks  following  repair  by  fastened  doubler  plates. 


1.  Stress  Check  is  a  trademark  of  Engineering  Software  Research  and  Development,  Inc. 
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intermediate  and  local  meshes  sufficient  to 
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Boundary  Conditions  (displacement,  load) 
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LUNCH  PRESENTATION 


/.  Morgan 

Oklahoma  City  Air  Logistics  Center 


B-IB  BOMBER  HORIZONTAL  STABILIZER 
SUBSTRUCTURE  FAILURES 


Prepared  By:  John  Morgan,  Aerospace  Engineer 
Structural  Engineering  Section 
B-IB  System  Support 
USAF 


INTRODUCTION 

In  May  of  1993,  a  B-IB  aircraft  sustained  a  lightning  strike  to  its  left  horizontal  stabilizer. 
In  an  attempt  to  identify  any  possible  damage,  the  stabilizer  was  X-rayed  at  various 
locations.  Some  of  these  X-rays  identified  cracks  in  the  substructure  that  could  not  be  the 
result  of  a  lightning  strike.  These  cracks  later  proved  to  be  a  part  of  a  larger  fleetwide 
problem  with  the  B-IB  aircraft  that  has  only  now  been  fully  understood. 

In  an  attempt  to  understand  and  control  the  problem,  the  US  Air  Force  formed  a  team  of 
engineering  experts  from  its  own  ranks  and  that  of  the  prime  contractor,  Boeing  North 
American  (BNA).  This  team  was  responsible  for  monitoring  the  status  of  the  fleet’s 
structural  integrity,  identifying  the  cause  of  the  problem  and  developing  a  structural 
enhancement  to  correct  the  problem.  The  USAF  is  currently  in  the  process  of 
implementing  the  structural  enhancement. 

The  purpose  of  this  paper  is  to  provide  an  overview  of  the  process  that  this  team  used  to 
deal  with  the  horizontal  stabilizer  problem. 

PHYSICAL  DESCRIPTION 

To  fully  understand  the  problem  with  the  B- IB  horizontal  stabilizers,  knowledge  of  the 
stabilizer  structure  is  required. 
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The  stabilizers  are  a  symmetric  airfoil  approximately  20  feet  long  (see  figure  1).  The 
width  tapers  from  16  feet  at  the  root  to  5  feet  at  the  tip.  Each  stabilizer  is  mounted  on  the 
aircraft  via  spindles  that  protrudes  from  the  side  of  the  aircraft.  The  stabilizers  can  be 
rotated  symmetrically  or  asymmetrically  about  these  spindles  to  control  both  the  pitch 
and  roll  of  the  aircraft. 

The  stabilizer  itself  consists  of  a  main  structural  box,  which  is  the  part  that  is  failing  and 
composite  fairings  which  are  attached  to  the  leading  edge,  trailing  edge  and  tip  of  this 
box. 

As  showu  in  figure  2,  the  main  structural  box  is  composed  of  relatively  thick(0.125  in.  to 
0.470  in)  one  piece  aluminum  upper  and  lower  skins  attached  to  a  thin  egg  crate  type 
substructure.  The  substructure,  as  shown  in  figure  3,  is  comprised  of  spars  (running  in 
the  spanwise  direction)  and  riblines  (running  in  the  chordwise  direction). 

The  front  and  rear  spars  are  machined  aluminum  I-beams.  Integral  flanges  are  machined 
into  the  spars  for  attachment  of  the  riblines.  The  intermediate  spars  are  titanium  I-beam 
weldments,  however,  the  webs  are  not  flat,  but  have  a  sinusoidal  shape.  These  spars  are 
commonly  refered  to  as  Sine  Wave  Beams(SWBs).  This  configuration  provides  increased 
strength  at  reduced  weight  when  compared  to  a  conventional  beam.  At  three  foot 
intervals  the  spar  web  flattens  out  to  provide  an  attach  point  for  the  riblines.  The  SWB 
spars  are  continuous  from  inboard  to  outboard  and  range  up  to  20  feet  in  length.  Typical 
thicknesses  are  0.012  in.  for  the  web  and  0.020  in.  for  the  cap 

The  riblines  are  composed  of  formed  aluminum  rib  segments  that  are  fitted  in  between 
the  spars.  These  rib  segments  are  of  a  two  piece  design:  (1)  A  rib  web  with  a  top,  bottom 
and  forward  attachment  flange,  and  (2)  an  aft  clip.  This  design  provides  only  fore  and  aft 
adjustment,  and  would  therefore  require  shimming  on  installation.  Thickness  for  the  ribs 
varies  from  0.040  to  0.125  in. 
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FIGURE  1 :  Relative  size  and  shape  of  the  horizontal  stabilizers  to  the  B-IB  aircraft. 
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Figure2.  B-!3  Ikii  izorital  Stabilizer 


FIGURE  3:  Typical  substructure  construction  of  the  outboard  section  of  the  B-IB 
Horizontal  Stabilizer. 


In  all,  the  design  of  the  stabilizer  provides  both  a  strong  and  lightweight  structure.  As 
stated  before,  the  stabilizers  have  the  strength  to  pitch  or  roll  a  477,000  lb  aircraft, 
however,  they  only  weighs  1655  pounds  each. 

FAILURE  BACKGROUND 

After  the  initial  discovery  of  the  damage  in  the  lightning  strike  stabilizer,  the  stabilizer 
was  subjected  to  a  detailed  inspection  The  first  finding  was  that  65  percent  of  the 
fasteners  attaching  the  top  skin  to  the  substructure  were  loose  or  missing  in  the  area 
where  the  initial  damage  was  discovered.  The  stabilizer  was  partially  disassembled  by 
removing  the  upper  skin  to  further  investigate  the  problem.  The  subsequent  visual 
inspection  identified  several  interesting  anomalies  (see  figure  4).  First,  the  damage, 
although  widespread  was  limited  to  the  outboard  two  thirds  (2/3)  of  the  main  structural 
box.  The  damage  itself  took  several  different  forms: 

•  Cracks  in  rib  segments 

•  Cracks  in  the  SWB  spars 

•  Cracks  in  the  integral  flanges  on  the  front  and  rear  spars 

•  Loose  or  failed  substructure-to-substructure  fasteners 

•  Degraded  corrosion  coating  on  aluminum  structure 

Lastly,  the  inspection  revealed  gaps  between  the  skins  and  the  substructure  that 
substantially  exceeded  the  original  production  drawing  tolerances  of  0.010  in..  These 
excessive  gaps  were  caused  by  three  different  phenomena 

•  Lack  of  shimming  to  level  adjacent  parts 

•  Premature  hardening  of  sealant  applied  to  the  weld  bead  that  resulted  in 
standing  off  the  skin 

•  Severe  weld  distortion  of  the  SWB  caps 

As  a  result  of  the  inspection  of  this  stabilizer,  the  engineering  team  turned  its  attention  to 
determining  the  extent  of  the  problem  in  theB-lB  fleet. 
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FIGURE  4  (Sheet  1  of  7):  Inspection  of  the 'i  Jghtning  Strike'  Horizontal  Stabilizer. 
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INSPECTION 


The  engineering  team  had  three  initial  goals  to  achieve  with  the  inspection  program.  The 
first  was  to  take  care  of  any  fleetwide  problem  with  loose  or  missing  skin  fasteners.  This 
was  accomplished  with  a  one-time  inspection  through  the  fleet  for  the  loose  or  missing 
fasteners.  Any  problems  encountered  were  corrected  in  conjunction  with  this  inspection. 

The  second  goal  was  to  determine  the  extent  of  substructure  damage  in  the  fleet.  To 
accomplish  this,  the  team  developed  a  borescope  and  X-ray  inspection  procedure  to 
identify  the  damage  in  the  substructure.  Inspection  of  a  sampling  of  the  aircraft  soon 
proved  that  the  problem  was  fleetwide,  therefore,  a  fleetwide  inspection  was  issued. 

This  inspection  showed  that  practically  all  stabilizers  had  some  form  of  substructure 
damage,  and,  in  a  large  majority  of  the  cases  the  damage  was  substantial  (see  figure  5). 

The  third  goal  was  to  assess  the  structural  integrity  of  any  stabilizer  exhibiting  damage 
and  make  a  recommendation  for  removal  from  service  or  reinspection  intervals.  This  was 
done  by  modifying  the  NASTRAN  finite  element  model  to  reflect  the  actual  damage, 
processing  for  internal  loads  and  routing  this  data  through  automated  stress  analysis 
programs  for  the  stabilizer  structure  to  determine  if  acceptable  residual  strength  remained. 

The  reliability  of  the  substructure  inspection  has  proven  to  be  less  than  desirable.  Visual 
inspections  of  some  stabilizers  disassembled  for  repair,  have  shown  that  less  than  half  of 
the  damage  noted  during  the  borescope/X-ray  inspections  was  actually  there.  In  addition, 
significant  damage  was  found  that  had  not  been  detected  in  the  borescope/X-ray 
inspections.  Contributing  to  the  low  reliability  is  the  size  of  the  inspection:  1 1 5  square 
feet  of  X-ray  film  and  200+  borescope  locations  for  each  stabilizer.  However,  the 
borescope/X-ray  inspection  remains  the  only  viable  technique  for  inspecting  the  stabilizer 
substructure. 
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FIGURE  5  (Sheet  2  of  2):  Quantifying  the  substructure  damage  in  the  fleet.  The 
quantities  reflect  the  number  of  stabilizers  in  the  fleet  that  contain  damage  in  a  particular 
SWB  Spar.  Quantities  shown  on  the  outside  periphery  of  the  diagram  and  centered  over 
the  ribline  locations  reflect  damage  in  integeral  flanges.  Based  on  195  stabilizers. 


SAFETY  OF  FLIGHT  ISSUE 


An  assessment  of  safety  of  flight  issues  regarding  an  inflight  failure  of  a  horizontal 
stabilizer  was  aided  by  a  1992  midair  collision  between  a  B-IB  and  a  KC-135  tanker 
during  a  midair  refueling.  The  B-IB  lost  the  outboard  2/3  of  the  left  hand  horizontal.  This 
damage  is  consistent  with  that  which  would  be  expected  due  to  failure  precipitated  by  the 
substructure  damage.  The  B-IB  pilot  reported  no  loss  in  flight  control  or  increase  in 
workload  and  the  aircraft  landed  safely  without  further  mishap.  The  B-IB  flight  simulator 
was  used  to  evaluate  the  aircraft/pilot  response  to  an  inflight  stabilizer  failure  during 
critical  high  workload  tasks.  These  findings  confirmed  that  the  aircraft  was  controllable 
and  the  loss  of  the  aircraft  was  not  probable. 


FAILURE  INVESTIGATION 

Metallurgical 

The  first  step  the  team  took  in  investigating  the  failures  was  to  evaluate  failed  samples 
from  the  lightning  strike  stabilizer.  The  metallurgical  results  showed  the  failure  mode  to 
be  high  cycle  fatigue  with  the  presence  of  a  high  mean  stress  -  cycle  counts  consistent 
with  a  high  frequency  event. 

Review  of  Original  Development  and  Testing 

With  the  metallurgical  results  in  mind,  the  team  began  reviewing  the  process  used  in  the 
original  design  development  and  testing  of  the  stabilizer. 

The  stabilizer  design  strength  had  been  verified  by  completing  a  full  scale  ultimate  load 
static  test  and  a  three  life,  flight  by  flight  fatigue  test  with  no  significant  failures.  A 
ground  and  flight  test  program  investigated  the  actual  loads  environment  of  the  aircraft 
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including  the  engine  acoustic  loading  on  the  empennage,  loads  envelope  and  actual 
maneuver  and  gust  spectra  loads.  No  surprises  were  found  during  this  test  program. 

The  review  obviously  keyed  on  environmental  aspects  that  could  affect  the  fatigue  life. 
There  were  no  apparent  problems  found  with  the  process,  however,  several  interesting 
facts  were  noted.  First ,  during  wind  tunnel  testing  of  the  aircraft,  a  vortex  was  found 
that  originated  high  on  the  forward  fuselage  and  passed  under  the  stabilizer.  This  vortex 
was  determined  to  be  quite  powerful  and  under  certain  conditions  could  significantly 
impact  the  loading  on  the  stabilizer.  Secondly,  due  to  the  size  of  the  stabilizers  and  the 
high  levels  of  the  engine  acoustic  noise,  no  full  scale  test  could  be  performed  to  verify 
the  acoustic  fatigue  life  of  the  stabilizer.  Since  the  design  process  assumed  that  the 
primary  type  of  response  would  be  a  high  frequency  panel  type  response,  test  boxes  were 
used  that  attempted  to  replicate  the  various  “cells”  of  the  stabilizer  structure..  Thirdly, 
some  early  acoustic  test  boxes  failed  prematurely  during  life  cycle  testing.  The  cause  was 
determined  to  out-of-tolerance  assembly  gaps  similar  to  those  found  in  the  lightning 
strike  stabilizer. 

Interim  Repair 

The  on  going  inspection  of  stabilizers  resulted  in  a  significant  number  being  withdrawn 
from  service  due  to  unacceptable  levels  of  damage.  With  no  spares  available,  this  was 
having  an  impact  on  the  readiness  of  the  B-IB  fleet.  With  the  information  concerning  the 
premature  failures  in  the  original  acoustic  test  boxes  and  the  results  from  the  inspections, 
a  preliminary  conclusion  was  reached  that  at  least  one  cause  of  the  substructure  failures 
were  assembly  anomalies.  Therefore,  an  interim  repair  was  developed  and  implemented 
that  would  correct  the  assembly  anomalies  and  replace  any  damaged  substructure.  It  was 
recognized  that  the  repair  might  not  be  the  final  solution  to  the  problem,  but  that  it  would 
maintain  the  serviceability  of  the  fleet,  until  further  investigation  and  testing  could 
validate  this  cause. 
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Verification  of  the  Cause  of  Failure 


With  the  pressure  off  to  take  care  of  the  immediate  stabilizer  problem,  efforts  turned  to 
investigating  and  verifying  the  causes  of  failure.  The  first  steps  in  this  process  were:  (1) 
verify  the  response  of  the  stabilizer  to  input  loads,  and  (2)  verify  the  source  and 
magnitude  of  cyclic  loads  on  the  stabilizer. 

Dynamic  Response:  To  verify  the  response  of  the  stabilizer  to  dynamic  input  loads, 
software  tools  and  hardware  capabilities  unavailable  during  the  original  development  of 
the  stabilizers,  were  used.  These  tools  allowed  detailed  Finite  Element  Models  (FEMs)  to 
be  developed  that  could  determine  the  dynamic  response  of  the  structure  to  specified 
inputs. 

The  first  use  of  these  models  were  used  to  determine  the  natural  frequencies  and 
associated  mode  shapes  of  the  stabilizer.  The  results  were  surprising  because  they 
showed  the  response  of  the  stabilizer  in  the  neighborhood  of  200  Hz(  where  the  engine 
acoustic  spectrum  has  a  peak)  consisted  of  complex  mode  shapes  with  the  shapes 
dependent  on  the  overall  geometry  of  the  stabilizer  acting  like  a  trapazoidal  plate  with  the 
skins  vibrating  in  unison.  As  previously  stated,  the  acoustic  test  box  for  the  stabilizer 
assumed  that  individual  panel  modes  (bounded  by  ribs  and  spars)  would  be  the  primary 
driver  for  the  stabilizer  response. 

Ground  and  Flight  Test:  The  purpose  of  the  ground  and  flight  tests  were  to  investigate 
any  unaccounted  for  load  conditions  for  the  stabilizers  and  to  verify  the  source, 
magnitude  and  response  characteristics  of  the  critical  cyclic  loads  on  the  stabilizer.  To 
accomplish  this  data  was  collected  from  an  instrumented  stabilizer  in  all  ground  and 
flight  regimes  that  might  contribute  to  the  stabilizer  loading.  No  significant  unaccounted 
for  load  conditions  were  found  during  this  test  program.  Although  the  test  verified  that 
the  acoustic  noise  of  the  engines  during  full  afterburner  use  during  the  take-off  roll  was 
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the  most  critical  condition,  the  most  interesting  result  was  the  verification  that  the 
stabilizer  responded  to  the  acoustic  input  as  predicted  by  the  dynamic  model. 

Dynamic  Modelin2 

With  the  discovery  of  this  new  acoustic  response  mode,  emphasis  was  put  into  improving 
the  accuracy  of  the  dynamic  model.  It  was  felt  that  this  model  could  be  used  to  quickly 
and  efficiently  model  both  the  current  configuration  and  new  configurations  of  the 
stabilizer.  Otherwise,  the  team  would  have  to  resort  to  slow  and  expensive  laboratory 
tests. 

Improvement  of  the  dynamic  model  consisted  of  refinements  of  the  mesh  and  changes  to 
modem,  higher  order  elements  that  would  more  accurately  model  the  complex  motion  of 
the  stabilizer  and  the  interaction  between  the  various  components.  Simulation  of  the 
acoustic  environment  and  coupling  of  those  loads  to  the  refined  NASTRAN  model  was 
performed  by  SDRC  who  supplies  an  I-DEAS  module  for  performing  vibro-acoustic 
analysis  (see  figure  6).  A  standand  NASTRAN  frequency  analysis  was  performed  with 
the  results  being  compared  with  data  from  ground  ranup  tests.  In  comparison  with  test 
results,  the  model  was  able  to  accurately  model  the  acoustic  input  from  the  engines  and 
the  response  of  the  stabilizer  (see  figure  7). 

Fatigue  Analysis 

The  dynamic  model  allowed  the  team  to  determine  the  stress  within  the  substmcture 
during  acoustic  excitement.  With  this  information,  the  team  could  finally  perform  a 
fatigue  analysis  of  the  substructure.  This  analysis  provided  disheartening  news  to  the 
team.  Acoustics  alone  could  not  have  caused  the  failures  seen  in  the  stabilizer. 

Therefore,  the  team  was  forced  to  look  for  other  contributing  factors. 


996 


oinparlson  of  Acceleration  Heoponsos 

ctober  1995  Ground  Runup  vs.  NASTRAN  S/N  50  Repair  Configuration 
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The  team  did  not  have  to  look  far.  A  high  enough  residual  stress  in  the  substructure,  due 
to  attempting  to  pull  out  gaps  during  assembly,  could  cause  fatigue  failure  when  loaded 
by  the  acoustic  noise  (see  figure  8).  A  few  analytical  models  and  verification  tests, 
proved  that  this  was  the  case.  In  fact,  at  certain  locations  pulling  out  gaps  well  within  the 
tolerances  provided  on  the  engineering  drawings  would  provide  sufficient  residual 
stressed  to  cause  a  fatigue  failure  It  should  be  noted  that  once  the  contribution  of  the 
residual  stresses  were  appreciated,  the  stabilizer  was  analyzed  for  the  flight  load  spectrum 
and  found  to  have  more  than  adequate  life. 

Failure  Investigation  Summary 

The  final  determination  of  the  investigating  team  was  that  the  vast  majority  of  the 
damage  noted  in  the  substructure  was  due  to  a  combination  of  high  residual  stress  and 
acoustics.  The  high  residual  stresses  were  due  to  pulling  out  assembly  gaps  during 
fastener  installation. 

Cure  for  the  Problem 


The  team  had  two  options  for  curing  the  problem  of  the  substructure  failures:  (1)  Lower 
the  residual  stress,  (2)  Lower  the  cyclic  stress.  Improvements  in  design,  fabrication  of  the 
SWBs  and  assembly  techniques  have  been  shown  to  significantly  reduce  substructure 
gaps,  but  sufficient  residual  stresses  can  be  reached  by  pulling  out  allowable  assembly 
gaps.  The  team  quickly  agreed  that  reducing  this  stress  was  not  the  sole  solution  to  the 
problem  Therefore,  the  only  option  available  was  to  lower  the  cyclic  stress  in 
conjunction  with  controlling  the  residual  stresses.  This  could  be  achieved  by  stiffening 
the  substructure.  The  interim  repair  arbitrarily  stiffened  the  riblines,  however,  stopped 
short  of  changing  the  spars.  The  final  repair  that  was  developed,  included  stiffening  the 
sparlines.  This  repair  has  been  proven  analytically  to  solve  the  problem.  Currently,  the 
team  is  testing  a  prototype  of  the  repair  to  verify  the  analytical  findings.  Plans  have  been 
enacted  to  implement  the  final  repair  beginning  in  the  year  2000. 
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FIGURE  8  (Sheet  1  of  2):  Predicted  life  of  one  rib  segment  in  the  original  production 
stabilizer.  The  improved  repair  lowers  the  cyclic  stress  in  the  rib,  therefore,  provides 
significant  improvement  in  life. 


SUMMARY 


The  substructure  failure  of  the  B-IB  horizontal  stabilizers  provided  a  challenging  and 
interesting  problem  for  the  investigating  team.  In  addition,  it  provided  some  important 
lessons  learned.  For  one,  the  failures  showed  the  possible  effect  of  using  such  contrasting 
structure  (ie.  thin  substructure  beneath  relatively  thick  skins).  Also,  the  failures  show 
how  residual  stresses  occurring  from  common  assembly  processes  can  cause  problems. 
Therefore,  these  stresses  need  to  be  dealt  with  during  design  and  the  development  of  the 
manufacturing/assembly  plan  Lastly,  the  use  of  an  important  new  engineering  tool  was 
demonstrated  during  the  failure  investigation.  This  tool  is  the  capability  of  new  software 
to  accurately  model  acoustic  sources  and  the  dynamic  responses  of  complex  structure. 
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Background/Introduction  Cont’d 

Horizontal  designed  for  high  cycle  vibratory  stresses 


O) 

c 


o 

E 
a 
o 

0) 

0) 

o  ■o 

^  O) 

s  -g  « 


0) 
0 

75 
£ 
O 
3 


0 

To 

E 


■o 

c 

0 

7? 

0 

> 


1  = 

■o 

Q)  0 
0  0 

i  ? 

0  0 

•S  -52 

O)  0 

c  >» 

®  75 
£  £ 
O)  0 

S  0 

o  -s 
0 

■O  0 

-S 

o  0 

X  o 
LU  O 


■D 
■D 
0 

0 
■D 
3 

O 
£ 

S' 

o> 
o> 

0  T- 


0 

3 

O) 


■o 

0 

■o 

3 

o 


0 

E 

a. 

o 

0 

> 

0 

0  o 
o  — 
0  o 


> 

3 


o> 

o> 


o 

0 


O) 

£ 


■o 

0 

O) 

£ 

0 

o 


S  0 

£  O 
£  :£ 

2  ■£ 
>  .E> 
£  0 
0  0 

«  s 

0  £= 


0 

3 

0 

0 

+-  >- 
C  0 


0 


I  I  .? 


0 
HH* 

U.  0 


£  0 

£  g 

<D  0 

O  "O 

0  c 
OC  0 


0 

> 

0 

£ 

0 


T5  i5 


O 
o 

®  ■£» 
♦2  O) 
0  ™ 

£  'V 
■o  c 
c  T 

o  Z 
o 

D)m- 


0 

*-« 

0 

0 


o 

o  0 
o  g 

0 

S  i 

£  JQ 
O  O 

§■  « 
0  0 


.=  O  ^ 
0  ®  0  ? 

H  S  H  < 

I  I  I 


1007 


Predict  responses  close  to  those  measured  during  test 
Identify  damage  mechanism(s) 

Develop  analytic  acoustic  environment  independent  of 
structural  models 


B-1B  Horizontal  Stabilizer 
Structural  Breakdown 
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Outb’d  Rib  (Typical) 


Horizontal  Rib  Configuration  and 
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Horizontal  Stabilizer 
esentative  Soar  Dama 
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External  and  internal  acoustic  measurements 

External  accelerometers 

Limited  number  of  internal  strain  gages 


Significant  Resuits  From  Test  Data 
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Low  Frequency  Analysis 
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Displacements  are  a  function  of  actual  structure 
Method  results  in  high  stresses  at  enforced  points 
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Internal  Acoustic  Analyses  Cont’d 
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High  Frequency  Models  and 
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Spars  were  shear  elements  with  shear  modulus  reduced  20% 
to  simulate  reduced  shear  stiffness  of  sine  wave  beams 
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Coupled  Acoustic/Structural 

Model  Analysis 

Boeing  North  American  using  i-Deas™  software 
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Deveioped  in  France  by  STRACO 

Acoustic  field  is  modeied  by  boundary  eiement  model 
(BEM) 

Determines  acoustic  and  eiasto-acoustic  modes 


Predicted  Vs  Measured  Acoustic 
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Measured  top  surface  levels  10  dB  below  bottom 
Therefor  accurate  top  surface  modeling  not  required 
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Enhanced  NASTRAN  Model 
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Enhanced  NASTRAN  Model  Cont’d 
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96  Mesh 


BNA  Response  Analysis 
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Predicted  Vs  Measured  Stress 
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repaired  horizontal  stabilizers  to  meet  their  remaining  life 
requirements 


structural  Integrity  Computer  Program 
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Structural  Integrity  Computer  Program 
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Structural  Integrity  Computer  Program 
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-  Ratio  of  Steady  to  Applied  Stress  Highly  Scattered 

•  COMPLEX  ROTATING  COMPONENTS  GEOMETRIES 
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structural  Integrity  Computer  Program 
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Structural  Integrity  Computer  Program 
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structural  Integrity  Computer  Program 
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Structural  Integrity  Computer  Program 
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;  all  Above  Contract  Sponsored  by  WR-ALC 
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Structural  Integrity  Computer  Program 
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Structural  Integrity  Computer  Program 
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Common  Error  Trapping  (Force  Total  Usage  =  100%  ...) 
Will  Allow  User  to  Control  Where  Adjustments  are  Made 
Eliminate  Time  Allocated  to  Unachievable  Regimes 
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Structural  Integrity  Computer  Program 
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•  STUDY  POTENTIAL  of  SICP  USE  with  HUMS  DATA 


SESSION  XII 

FORCE  MANAGEMENT 

Chairman  -  J,  Turner 
San  Antonio  Air  Logistics  Center 


Detecting  High  Cycle  Fatigue  with  User  Defined  Regime  Recognition 


John  A.  Cicero,  Ph.D. 

Senior  Software  Engineer 
Systems  &  Electronics,  Inc. 

190  Gordon  Street 
Elk  Grove  Village,  IL  60007-1 120 
Tel:  (630)  829-6556 
FAX:  (630)  829-6551 
E-mail:  jcicero@ben.edu 


Introduction 

Systems  &  Electronics,  Inc.  (SEI)  has  developed  a  system  that  detects  high 
cycle  fatigue  through  user  defined  regime  recognition.  The  system  includes  an 
airborne  recorder  and  a  ground-based  PC.  The  airborne  recorder  records  aircraft 
parameters  such  as  airspeed,  altitude,  Nz,  angular  acceleration,  etc.  This  data  is 
compressed  and  then  stored  in  an  electronic  (FLASH)  memory  module.  At  the 
end  of  a  flight  (or  series  of  flights)  the  memory  module  is  downloaded  into  a  PC. 
A  decompression  algorithm  is  used  to  expand  the  compressed  data  into  time- 
contiguous  data.  This  time-contiguous  data  is  then  processed  by  a  flight 
condition  code  processor  program.  This  program  is  used  to  identify  flight 
maneuvers  that  cause  fatigue. 

Classifying  the  Maneuvers 

What  makes  this  system  unique  is  that  the  user  can  classify  a  specific 
maneuver  by  its  individual  parameter  ranges.  For  example,  assume  that  the  user 
defines  an  Angle  Of  Bank  50  Degree  Left  Turn  @  0.4Vh  as  one  in  which: 

0.35%  <  airspeed  (%Vh)  <  0.45% 

95%  <  rotor  RPM  <  1 1 5% 
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20%  < 
0.5g< 
-60  degrees  < 
-45  degrees  < 
20  %  < 
20  %  < 
20  %  < 
-500  feet/min  < 


engine  torques 
Nz 

Roll  Attitude 
Pitch  Attitude 


<110% 

<3g 

< -6  degrees 
<  45  degrees 


Latitudinal  Stick  Position  <  100% 
Longitudinal  Stick  Position  <  100% 
Rudder  Stick  Position  <  1 00% 

Rate  of  Climb  <  500  feet/min 


Using  the  flight  condition  code  processor  program  shown  in  Figure  1,  the  user 
can  create  a  new  flight  condition  code  or  modify  an  existing  flight  condition  code. 
For  each  input  parameter  the  user  selects  whether  or  not  the  input  is  used  to 
classify  the  current  flight  condition  code.  If  it  is  selected,  the  user  must  determine 
the  lower  and  upper  limits  of  this  input  parameter.  For  any  given  flight  condition 
code  the  user  can  also  specify  that  the  following  additional  flight  information  be 
logged: 


1)  the  parameter  value  at  the  start  of  the  maneuver, 

2)  the  parameter  value  at  the  end  of  the  maneuver, 

3)  the  maximum  parameter  value  during  the  maneuver,  and 

4)  the  minimum  parameter  value  during  the  maneuver. 

The  user  can  also  specify  discrete  event  information  such  as  a  takeoff  and  a 
landing.  The  flight  condition  code  processor  can  be  configured  to  first  look  for  a 
takeoff  before  it  will  recognize  a  landing.  After  a  takeoff  occurs,  the  flight 
condition  code  processor  can  be  configured  to  look  for  a  landing  before  it  will 
recognize  another  takeoff 

The  Flight  Condition  Code  Processor  Output 

Once  all  of  the  flight  condition  codes  have  been  entered,  the  user  can  select  the 
ran  option  from  the  flight  condition  code  processor  menu.  There  are  essentially 
two  sets  of  input  data  to  the  flight  condition  recognition  program.  One  set  is  the 
decompressed  time-contiguous  flight  data  from  the  flight  recorder.  The  other  set  is 
the  user  specified  flight  condition  codes  described  above.  After  the  data  is  run 
through  the  flight  condition  code  processor  the  following  output  sets  will  be 
generated: 
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Figure  1 .  The  flight  condition  code  processor. 


The  first  output  is  a  flight  profile.  When  the  flight  processor  encounters 
recorded  data  within  the  parameter  limits  (that  were  determined  by  the 
user)  it  will  report  the  following  information: 

1)  a  unique  flight  condition  name  (i.e.,  Angle  Of  Bank  50  Degree  Left  Turn 
@  0.4Vh), 

2)  the  time  duration  of  the  flight  condition  (i.e.,  10  seconds), 

3)  any  coincident  peak/valley  parameter  data  that  occurred  during  the 
condition  (i.e..  If  Nz  was  chosen  as  a  coincident  peak/valley  parameter,  the 
report  might  contain  the  following  results:  Nz  valley  =  0.8  g’s  and  Nz 
peak  =  2.6  g’s), 

4)  any  coincident  parameter  data  recorded  at  the  start  of  the  flight  maneuver, 
(i.e..  If  Rate  of  Climb  is  chosen  as  a  coincident  parameter  at  the  start  of  the 
flight  maneuver,  the  report  might  contain  the  following  results:  Rate  of 
Climb  Entering  =  -300  feet/'minute),  and 
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Figure  2.  A  histogram  illustrating  the  amount  of  time  spent  in  various  maneuvers. 


5)  any  coincident  parameter  data  recorded  at  the  end  of  the  flight  maneuver, 
(i.e..  If  Rate  of  Climb  is  chosen  as  a  coincident  parameter  at  the  end  of  the 
flight  maneuver,  the  report  might  contain  the  following  results:  Rate  of 
Climb  Leaving  =  -250  feet/minute). 

The  second  output  from  the  flight  condition  code  processor  describes  the  amount 
of  time  spent  in  each  maneuver.  These  results  can  be  placed  into  a  spreadsheet  to 
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generate  a  histogram  of  the  amount  of  time  spent  in  the  most  frequently  occurring 
flight  maneuvers.  For  example,  the  histogram  in  Figure  2  shows  that  the  aircraft 
spent  a  total  of  300  seconds  in  Straight  &  Level  Flight  @  O.SVh.If  the  user  is  not 
satisfied  with  the  results  from  the  system,  the  user  can  modify  the  classification  of 
flight  conditions,  and  then  rerun  the  existing  data  through  the  flight  condition 
code  processor  to  obtain  a  new  set  of  results.  It  is  important  to  note  that  the 
aircraft  does  not  have  to  be  flown  again  to  rerun  the  flight  data  through  the  fine- 
tuned  flight  condition  code  processor. 

Neural  Network  Generated  Loads 

In  another  version  of  the  system  the  recorded  data  is  run  through  neural 
network  algorithms  developed  by  the  Navy  to  obtain  predicted  loads  at  various 
points  on  the  aircraft.  This  data  is  also  compressed,  stored  in  an  electronic 
memory  module,  downloaded  into  a  PC,  and  then  expanded  for  processing.  This 
data  is  then  run  through  the  flight  condition  recognition  algorithms  described 
above.  These  algorithms  can  identify  load  conditions  (in  addition  to  flight 
conditions)  that  cause  fatigue. 

The  output  of  the  load  information  is  displayed  by  a  viewer.  The  viewer  (Shown 
in  Figure  3.)  allows  the  user  to  graphically  compare  flight  loads  with  input 
parameters  at  various  times  during  the  flight.  Each  of  these  graphs  are  time 
stamped  and  flight  condition  code  stamped  so  that  the  user  can  see  what  type  of 
loads  occur  during  a  particular  flight. 

Conclusion 

SEI  has  developed  a  set  of  tools  that  can  detect  high  cycle  fatigue  through  the 
following  steps: 

1 . )  Display  load  profiles  at  various  points  on  the  aircraft  during  an  entire 

maneuver. 

2. )  Provide  Peak/Valley  parameter  and  load  values  for  each  maneuver  during  a 

flight. 

3. )  Display  the  amount  of  time  spent  in  each  maneuver  during  the  entire  flight. 
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4.)  Examine  the  load  conditions  at  various  points  on  the  aircraft  during  the 
entire  flight. 


Figure  3.  The  load  viewer  tool. 
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Where  the  Requirements  Come  From 
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Principle  Structural  Element  (PSE) 
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NDE  Inspection  Development  Process 


1063 


Interaction  of  Signal/Noise 
Discrimination 


1065 


Xcw  121.688 


^  U 


ft  ^ 


1066 


Inspect  Forward  Lower  Wing  Skin,  Aft  Lower  Wing  Skin,  And 
Stringer  17  At  Shaded  Fastener  Locations 


Simulated 
Stringer  C/L 

(Upper  Edge  Of  -3  Detail,  Surface  Plate) 


1067 


Longitudinal  Splice  at  Longeron  5 
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Eddy  Current  C-Scan  Imaging 
Longitudinal  Splice  at  Longeron  5 
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DC-10  Crown  Skin  Splice 
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NDE  Inspection  Development  Process 
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C-17A  INDIVIDUAL  AIRCRAFT 
TRACKING  PROGRAM 

Rick  Selder 
C-17  Loads  and  Dynamics 
The  Boeing  Company 
USAF  Airlift  and  Tanker  Programs 

Ko-Wei  Liu 

C-17  Durability  and  Damage  Tolerance 
The  Boeing  Company 
USAF  Airlift  and  Tanker  Programs 


1997  USAF  Structural  Integrity  Program  Conference,  San  Antonio,  TX,  Dec.  2-4 


The  C-17A  Individual  Aircraft  Tracking  Program  (lATP)  is  being  developed  as 
part  of  the  C-17  full  scale  engineering  development  program.  lATP  computer 
system  programming  is  nearly  complete  with  full  operational  capability 
expected  by  mid  1998.  The  lATP  system  is  being  developed  by  the  USAF 
Airlift  and  Tanker  Programs  business  unit  of  the  Boeing  Company  in  Long 
Beach,  California. 
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Outline 


•  Benefit  of  iATP 

•  IATP  approach 

•  Improvements  over  past  lATPs 

•  Data  processing  flow 

•  SFDR  flight  recorded  parameters 

•  Flight  data  editing 

•  Stress  spectra  generation 

•  Control  point  selection 

•  Gap  filling  missing  data 

•  Future  Damage  Projection 

•  Content  of  output  lAT  data  report 

•  Computer  system  characteristics 
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Benefit  of  lATP 


•  Minimize  maintenance;  Inspection  of  critical  areas  are  defined  for 

each  aircraft  based  on  its  individual 
damage  assessment. 

•  Efficient  fleet  management;  Modifications  are  scheduled  and  aircraft 

base  rotations  and  retirements  are  based 
on  actual  aircraft  usage. 

•  Extended  service  life;  Appropriate  inspections  and  efficient  fleet 

management  will  maximize  the  service  life 
of  each  aircraft. 

•  Better  field  support;  lATP  data  used  to  help  identify  parameters 

that  attribute  to  problems  found  in  the  field. 
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Today,  the  cost  of  buying  and  maintaining  military  aircraft  is  very  high. 
Therefore,  there  is  a  great  need  to  extend  the  service  life  of  aircraft  as  long  as 
possible  but  without  jeopardizing  safety  due  to  long  term  structural  degradation 
of  the  airframe  as  a  result  of  application  of  repeated  loads.  An  I  ATP  which 
gives  the  ability  to  track  the  structural  status  of  individual  aircraft  is  the 
primary  tool  to  assure  a  long  service  life  without  risk.  The  main  benefits  are 
described  as  follows; 

Minimize  maintenance:  The  structural  status  of  each  aircraft  is  used  to  define 
critical  areas  and  inspection  intervals  for  each  aircraft  based  on  its  individual 
damage  assessment.  This  will  minimize  required  maintenance  without 
jeopardizing  aircraft  safety. 

Efficient  flept  manapernent:  Modifications  for  each  aircraft  can  be  scheduled 
based  on  its  unique  structural  assessment.  Also,  decisions  on  rotating  aircraft 
between  bases  or  retiring  aircraft  can  be  based  on  actual  usage  statistics. 

Extended  service  life:  Performing  inspection  at  appropriate  intervals  and 
efficient  management  of  the  fleet  will  maximize  the  service  life  of  each 
aircraft. 

Retter  field  support:  The  lATP  system  gives  valuable  flight  data  that  can  be 
used  to  help  solve  problems  found  in  the  field.  For  example,  if  a  problem  was 
found  to  be  associated  with  an  incident  that  occurred  during  a  specific  flight, 

I  ATP  flight  data  can  be  evaluated  to  find  the  circumstances  surrounding  the 
incident. 
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lATP  Approach 


•  Record  aircraft  CG  Nz  cycles  and  associated  flight  parameters  for 
each  flight  of  every  aircraft. 

•  Apply  analytic  methods  to  recorded  data  to  generate  flight-by-flight, 
cycle-by-cycle  stress  spectra  at  selected  airframe  control  points. 

•  Calculate  damage  based  on  cycle-by-cycle  crack  growth  analysis 
at  each  control  point. 

•  Compile  usage  statistics  by  aircraft,  duty  base,  and  fleet. 


The  purpose  of  an  lATP  system  is  to  calculate  damage  at  selected  airframe 
locations  in  order  to  determine  differences  in  aircraft  usage  and  to  adjust 
inspection  intervals  of  critical  areas  of  the  airframe.  To  accomplish  this,  the 
approach  used  for  the  C-17  starts  by  recording  aircraft  vertical  acceleration 
(Nz)  cycles  and  associated  flight  parameters  for  each  flight  of  every  aircraft. 
Analytical  methods  are  then  applied  to  calculate  flight-by-flight,  cycle-by¬ 
cycle  stress  spectra  at  selected  airframe  locations  (control  points).  Damage  at 
each  control  point  is  then  calculated  based  on  cycle-by-cycle  crack  growth 
analysis.  In  addition,  usage  statistics  are  compiled  for  each  aircraft,  by  duty 
base,  and  for  the  C-17  fleet  as  a  whole. 
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Improvements  Over  Past  lATPs 


•  Eliminate  dependency  on  the  use  of  strain  gauges. 

°  Past  experience  shows  that  strain  gauges  are  prone  to  faults. 

o  Production  strain  gauges  expensive  to  install  and  maintain. 

°  Gives  flexibility  in  choosing  or  changing  control  points. 

•  Damage  determined  directly  from  recorded  flight  profiles  rather 
than  assumed  to  be  equivalent  to  the  damage  from  the  closest 
matching  design  mission  type. 

•  Reduces  flight  crew  work  load:  Uses  a  Standard  Flight  Data 
Recorder  (SFDR)  eliminating  the  need  for  flight  logs. 

•  Spectra  generation  methods  are  as  rigorous  and  detailed  as  that 
used  for  design  analysis  of  the  C-17  airframe. 


The  C-17  I  ATP  has  many  improvements  over  lATPs  used  in  the  past.  These 
improvements  are  described  as  follows: 

•Control  point  stress  spectra  generation  does  not  depend  on  the  use  of  strain 
gauges.  Production  strain  gauges  are  expensive  to  install  and  maintain  and  past 
experience  shows  that  they  are  prone  to  faults.  Also,  eliminating  the  need  to 
locate  a  gauge  at  each  control  point  location  gives  greater  flexibility  in 
choosing  or  changing  control  points. 

•Some  past  lATPs  obtain  the  damage  due  to  a  given  flight  by  assuming  this 
damage  is  equivalent  to  that  given  by  the  closest  matching  design  mission  type. 
The  C-17 1  ATP  determines  damage  due  to  a  given  flight  directly  from  its 
recorded  flight  profile  parameters. 

•All  flight  data  needed  by  the  C-17 1  ATP  is  recorded  by  a  Standard  Flight  Data 
Recorder  (SFDR).  This  eliminates  the  need  to  collect  flight  log  data  thus 
reducing  the  work  load  of  the  flight  crew. 

•  The  SFDR  records  28  separate  flight  parameters  allowing  for  the  generation 
of  flight  profiles  as  detailed  as  the  design  mission  profiles  used  for  design 
analysis  of  the  C-17  airframe.  This  then  allows  for  a  rigorous  generation  of 
stress  spectra. 
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Data  Processing  Flow 


The  goal  of  I  ATP  is  to  collect  and  process  100%  of  operational  flight  data  for  every 
aircraft  in  the  C-17  fleet.  Flight  data  is  recorded  by  use  of  a  digital  Standard  Flight 
Data  Recorder  (SFDR),  produced  by  Smiths  Industries,  installed  on  each  aircraft. 

Data  collection  starts  by  downloading  compressed  flight  data  from  the  SFDR  onto 
floppy  diskettes  by  the  AMC  maintenance  crews.  The  data  is  then  sent  to  the 
Oklahoma  City  Air  Logistics  Center  (OC-ALC)  where  it  is  uploaded  and  archived  onto 
an  IBM  mainframe  computer.  Then,  before  data  can  be  processed  by  the  lATP  system, 
it  is  decompressed  by  use  of  decompression  software  supplied  by  Smiths  Industry. 

The  first  step  of  the  lATP  process  is  to  edit  the  flight  data  which  includes  reformatting 
data  into  separate  flights,  parameter  editing,  and  segment  identification.  Flights 
containing  erroneous  data  that  cannot  be  corrected  are  discarded.  Once  complete  the 
resulting  flight  profiles  and  Nz  sequences  are  stored  into  a  database.  In  addition  to 
SFDR  data,  maintenance  records  and  horizontal  stabilizer  and  duty  assignments  are 
entered  into  the  database  by  the  user.  All  valid  flights  are  then  processed  through  a 
spectra  generation  program  to  produce  control  point  stress  spectra. 

Once  all  stress  spectra  have  been  generated  for  all  fights  in  a  given  report  period,  they 
are  processed  through  a  damage  program  to  calculate  damage  at  each  control  point. 
Damage  for  missing  flights  is  calculated  using  a  gap  filling  process.  In  addition,  usage 
statistics  are  compiled  for  the  report  period  and  accumulated  with  all  previous  periods. 
These  usage  statistics  and  damage  results  are  then  output  as  the  lAT  data  report. 

The  lATP  system  also  provides  a  backup  utility  to  archive  edited  flight  data  and  to 
delete  no  longer  needed  data  following  report  generation. 
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SFDR  Parameters 

H  Mmr 

Data  Taq 

Event  Parameters 

Aircraft  Serial  Number 

Auto-oilot  On 

Seauential  Fliaht  Number 

Thrust  Reverser  Deoloved 

Airframe  Time 

Slats  Out 

Mission  Tvoe 

UARRSI  Ready 

Date  of  Fliaht 

UARRSI  Latched 

Deoarture  Base 

Carao  Door  Open 

Landina  Gear  Ud 

Time  Tracked 

Weiaht  On  Wheels 

Pressure  Altitude 

Parkina  Brake  Enaaaed 

Radar  Altitude 

Flao  Deoloved 

Indicated  Airsoeed 

Soeed  Brake  Deoloved 

Gross  Weiaht 

TrooD  Door  Open 

Fuel  Weiaht 

Aircraft  CG 

Peak/Vallev  Time  History 

Enaine  Pressure  Ratio 

ICG  Vertical  Acceleration  (Nz) 

Flao  Position 

Wheel  Speed 

$IJDE7 

Four  types  of  data  are  recorded  by  the  SFDR.  These  are  data  tag,  event,  time 
tracked,  and  peak/valley  time  history  parameters.  Data  tag  parameters  give 
information  on  the  flight  as  a  whole  and  are  recorded  once  at  the  beginning  of 
each  flight.  Event  parameters  identify  a  change  in  aircraft  configuration  and 
are  recorded  when  a  change  in  configuration  occurs.  Time  tracked  parameters 
change  steadily  with  time  and  are  recorded  whenever  their  value  changes  by  a 
predetermined  increment  (gate  value).  Peak/valley  tinie  history  parameters  are 
cyclic  in  nature  and  are  processed  by  an  onboard  algorithm  which  finds  and 
records  peak/valley  cycles  outside  a  predetermined  threshold. 
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Flight  Data  Editing 


•  Reformatting; 

°  Divide  data  into  individual  flights. 

°  Check  continuity  of  aircraft  serial  and  flight  numbers,  airframe 
time,  and  flight  date  against  previously  processed  flights. 

•  Individual  Parameter  Editing; 

°  Check  if  parameter  value  is  within  an  acceptable  range. 

“  Check  if  rate  of  change  of  value  is  within  an  acceptable  slope. 

°  Check  if  parameter  was  recorded  at  predetermined  increments 
(gate  value) 

°  Check  if  event  parameter  (flaps  deploy,  gear  extend,  etc.)  cycles 
within  an  acceptable  time  span. 
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Flight  data  editing  is  accomplished  in  four  steps  described  as  follows; 
Reformatting: 

One  download  of  SFDR  data  contains  several  flights.  This  step  divides  the 
data  into  the  separate  flights  then  checks  the  continuity  of  aircraft  serial  and 
flight  numbers,  airframe  time,  and  date  against  the  same  data  for  previously 
processed  flights. 

Individual  parameter  editing: 

Each  parameter  is  checked  to  assure  that  its  value  is  within  an  acceptable  range 
and  its  change  in  value  is  within  an  acceptable  slope.  Also,  the  parameter 
value  is  checked  if  it  was  recorded  at  the  predetermined  increment  (gate  value), 
(example:  gross  weight  is  checked  to  assure  it  is  recorded  in  640  lb, 
increments.)  For  event  parameters,  a  check  is  made  to  assure  that  the  event 
cycle  time  is  within  an  acceptable  time  span,  (example:  flaps  deploy  within  30 
sec.) 
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Flight  Data  Editing _ 

•  Parameter  Comparison  Checks; 

°  Compare  parameters  against  each  other  for  conflicts. 

•  Flight  Profile  Generation: 

°  Divide  flight  into  identifiable  segments  (ascent,  cruise,  air-drop, 
etc.) 

°  Check  segment  sequence  for  logical  flow. 

°  Find  special  event  segments,  (example:  Touch  &  Go) 

°  Categorize  flight  as  one  of  35  design  mission  types.  (Provided  to 
give  usage  information,  not  used  in  damage  calculation) 

•  Note:  Methods  incorporate  experience  gained  as  a  result  of 
processing  over  4,000  flights  of  actual  SFDR  data. 


Parameters  comparison  checks: 

Parameters  are  compared  against  each  other  to  find  any  conflicts,  (example: 
Landing  gear  extended  above  a  placard  speed  of  250  knots) 

Flight  profile  generation: 

The  time  history  of  recorded  flight  parameters  is  divided  into  identifiable 
segments  and  then  the  value  of  each  parameter  within  the  time  span  of  the 
segment  is  averaged.  The  exception  is  peak/valley  time  history  parameters  for 
wWch  all  values  recorded  within  a  segment  are  retained.  The  resulting  profile 
of  flight  segments  is  then  checked  to  assure  that  the  segments  are  sequenced  in 
a  logical  flow.  Also,  the  profile  segments  are  surveyed  to  find  combinations 
thamake  up  special  event  segments.  For  example,  the  combination  of  a 
landing  impact,  followed  by  a  landing  roll,  take-off  run,  then  take-off  rotation^ 
is  identified  as  a  touch  and  go  segment.  Once  the  flight  profile  is  generated  it  is 
categorized  as  one  of  35  design  mission  types  by  matching  its  flight 
characteristics  to  that  of  the  design  mission  types. 

It  should  be  noted  that  the  above  described  data  editing  methods  were  not 
developed  solely  on  a  preconceived  idea  of  the  nature  of  flight  recorded  data, 
but  rather  incorporate  the  experience  gained  and  are  proven  out  as  a  result  of 
processing  over  4000  flights  of  actual  SFDR  recorded  data. 
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Flight  Data  Editing 


Flight  Profile 


ft/C 

FLT  AIFRftME 

TIME  DATE 

BASE 

GMT 

9306D3 

137  375.250 

379.125  950814 

KADU 

1601  1947 

DUR 

SEGMENT 

PALT  RflLT  IAS 

GU 

FU 

CG 

EPR 

USPD 

FLAP 

594.250 

FULL  STOP 

240.  0.  0. 

503249. 

174481. 

36.2 

1.0177 

0.00 

O.QQ 

2.GQ0 

TAXI 

240.  0.  0. 

502080. 

173600. 

36.3 

1 . 4336 

18.00 

18.25 

29.0Q0 

TAKEOFF  RUN 

240.  0.  0. 

502080. 

173600. 

36.3 

1 . 4336 

77.79 

18.25 

G.OOO 

TAKE  OFF  ROTATION 

240.  0.  148. 

5Q2Q80. 

173600. 

36.3 

1 . 4336 

120.00 

18,25 

79.000 

ASCENT 

1031.  826.  187. 

502080. 

173600. 

36.3 

1.4336 

0.00 

18.25 

18.000 

ASCENT 

2028.  1956.  246. 

502080. 

173600. 

36.3 

1 . 4336 

0.00 

Q.QO 

1632.000 

ASCENT 

18549.  18531.  299. 

493516. 

165059. 

36.3 

1 . 2363 

0.00 

G.OO 

2134.000 

CRUISE 

27902.  28860.  300. 

480627. 

152157. 

36.3 

1.2187 

0.00 

0.00 

358.000 

ASCENT 

29468.  30557.  297. 

472888. 

144402. 

36.3 

1.3029 

0.00 

0.00 

7209.000 

CRUISE 

30980.  30627.  290. 

452072. 

123593. 

36.3 

1.2582 

0.00 

Q.QO 

1090.000 

CRUISE 

30992.  26932.  292. 

429336. 

100828. 

36.3 

1 . 2556 

0.00 

Q.QO 

460.000 

DESCENT 

21469.  13274.  306. 

425857. 

97518. 

36.3 

0.8463 

0.00 

0.00 

301.000 

CRUISE 

13811.  5545.  270. 

424868. 

96394. 

36.3 

1.0199 

0.00 

Q.QO 

68.000 

DESCENT 

12602.  5085.  253. 

424000. 

95400. 

36.3 

0.9697 

0.00 

0.00 

45.000 

DESCENT 

10433.  3851.  242. 

424000. 

95400. 

36.3 

0.9688 

0.00 

0.00 

17.000 

DESCENT 

9569.  3269.  212. 

424000. 

95400. 

36.3 

0.9688 

0.00 

18.50 

30.000 

LANDING  APPROACH 

8950.  2617.  187. 

424000. 

95400. 

36.3 

0.9688 

G.OO 

18.50 

140.000 

LANDING  APPROACH 

7547.  1254.  138. 

423771 . 

95177. 

36.3 

1.1308 

0.00 

38.25 

Q.QOO 

LANDING  IMPACT 

6296.  0.  128. 

423360. 

94800. 

36.5 

1.1797 

116.00 

38.25 

21.000 

LANDING  ROLL  OUT 

6296.  0.  6. 

423360. 

94800. 

36.5 

1.1015 

64.18 

38.25 

4.000 

TAXI 

6296.  0.  0. 

423360. 

94800. 

36.5 

1.0781 

14.00 

38.25 

185.750 

FULL  STOP 

6296.  0.  0. 

423360. 

94800. 

36.5 

1.0163 

0.00 

Q.QO 

The  above  is  a  generated  flight  profile  for  one  actual  flight.  The  first  record 
gives  information  on  the  flight  as  a  whole  with  each  of  the  following  records 
describing  each  segment.  Given  are  the  segment  duration  and  type  followed  by 
the  averaged  segment  parameters.  The  peak/valley  time  history  Nz  cycles 
associated  with  each  segment  are  stored  in  a  separate  unformatted  file. 
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Stress  Spectra  Generation 

•  Find  gust  encounters:  Separate  out  gust  Nz  cycles  and  time  in  gust 
(T J  by  use  of  the  following  criteria. 

°  Time  between  successive  Nz  peaks  and  valleys  less  than  1.0  sec. 

°  Time  between  recorded  threshold  crossings  (±0.1  g)  less  than  5.0 
sec.  Average  Nz  peak/valleys  time  difference  for  gust  encounter 
must  remain  less  than  1 .0  sec. 

o  Remaining  Nz’s  assumed  to  be  due  to  symmetric  maneuvers. 

•  Symmetric  maneuvers 

°  Obtain  flight  parameter  to  control  point  stress  transfer  coefficients 
by  performance  of  a  regression  analysis. 

°  Apply  transfer  coefficients  to  flight  segment  parameters  to  obtain 
stress  versus  Nz  in  terms  of  slope  (F)  and  intercept  (P). 

°  Convert  flight  segment  maneuver  Nz  cycles  to  stress  cycles: 
a  =  P  +  Nz  X  F 


For  each  flight  profile  segment  control  point  stress  cycles  due  to  loading  events 
occurring  within  the  segment  are  determined  using  that  segment’s  parameters 
and  associated  Nz  cycles.  The  loading  events  include  symmetric  maneuvers, 
atmospheric  turbulence,  dynamic  taxi,  landing  impact,  and  take-off  rotation. 

During  flight  segments  both  symmetric  maneuver  and  atmospheric  turbulence 
(gust  encoimter)  occur.  The  gust  encounters  are  found  by  applying  the 
following  criteria  to  the  recorded  Nz  cycles  to  separate  out  gust  Nz  cycles  and 
time  in  gust  (Tg): 

1 .  The  gust  encounter  begins  when  the  time  between  successive  Nz  peaks  and 
valleys  is  less  than  1.0  sec. 

2.  The  encounter  continues  if  the  above  remains  true  or  if  the  time  between 
recorded  threshold  crossing  is  less  than  5.0  sec.  and  the  average  Nz 
peak/valley  time  difference  remains  less  than  1.0  sec. 

Note:  the  5.0  sec  criteria  is  used  to  allow  for  reasonable  continuation  of  a  gust 
encounter  despite  dropped  data  values  due  to  the  recording  threshold  of  ±0.1  g. 

All  Nz’s  that  do  not  meet  the  above  criteria  are  assumed  to  be  due  to  symmetric 
maneuvers. 

Symmetric  maneuvers: 

For  each  maneuver  Nz  value,  control  point  stress  is  calculated  using  the 
following  linear  relationship: 

a  =  P  +  Nz  x  F 
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Stress  Spectra  Generation 


•  Atmospheric  turbulence 

°  Select  one  of  60  PSD  gust  cases  by  matching  on  flight  segment 
parameters. 

°  Use  gust  case  analytic  stress  response  (A,  NJ  and  segment 
time  in  gust  (T^)  to  generate  a  random  sequence  of  delta  stress 
(Ao)  cycles  by  application  of  Monte  Carlo  methods. 

°  Add  f  .Og  stress,  calculated  using  maneuver  stress  equation,  to 
Aa  cycles. 

°  Future  enhancements: 

•  Evaluate  C.G.  Nz  to  distinguish  between  non-storm  and  storm 
turbulence,  (currently  an  analytic  distribution  is  assumed) 

•  Evaluate  L/ESS  strain  gage  data  to  develop  a  relationship 
between  C.G.  Nz/Nx  and  stress  response. 
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Symmetric  maneuvers  rcont.~): 

The  intercept  (P)  and  slope  (F)  of  this  equation  were  determined  by  applying  stress 
transfer  coefficients,  determined  by  performance  of  a  regression  analysis,  to  the 
flight  segment  parameters. 

Atmospheric  turbulence: 

The  first  step  in  determining  stress  cycle  for  atmospheric  turbulence  is  to  select  one 
of  60  PSD  gust  cases  by  performing  a  best  match  of  flight  segment  to  gu^  case 
parameters.  Then  from  the  chosen  gust  cases  analytic  stress  response  (  A,  NJ  and 
the  flight  segments  time  in  gust  (TJ  generate  a  random  sequence  of  delta  stress  (Aa) 
cycles  by  application  of  Monte  Carlo  methods.  To  these  Aa  cycles  add  a  l.Og  stress 
calculated  using  the  maneuver  stress  equation. 

In  order  to  utilize  recorded  data  to  determine  the  gust  intensity  in  addition  to  the 
time  in  gust  encounters,  enhanced  methods  will  be  developed.  First,  an  effort  will  be 
made  to  develop  a  method  to  evaluate  C.G.  Nz  to  distinguish  between  storm  and 
non-storm  turbulence  (currently  an  analytic  distribution  is  assumed).  Following 
this,  when  a  sufficient  database  of  L/ESS  strain  gauge  data  is  available,  it  will  be 
evaluated  to  develop  a  C.G.  Nz  and  additional  Nx  (lateral  load  factor)  to  stress 
relationship. 
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Stress  Spectra  Generation 


•  Dynamic  taxi 

=  Select  one  of  25  dynamic  taxi  cases  by  matching  on  taxi  segment 
parameters. 

=  Use  taxi  case  analytic  stress  response  (a,  NJ  and  segment 
duration  (T)  to  generate  a  random  sequence  of  delta  stress  (Aa) 
cycles  by  application  of  Monte  Carlo  methods. 

=  Add  1  .Og  stress,  determined  by  application  of  regression  transfer 
coefficients  to  taxi  segment  parameters,  to  Aa  cycles. 

=  Future  enhancements: 

.  Evaluate  C.G.  Nz  to  distinguish  between  prepared  and  semi 
prepared  airfields,  (currently  all  airfields  assumed  to  be 
prepared) 

.  Evaluate  liESS  strain  gage  data  to  development  a  relationship 
between  C.G.  Nz  and  control  point  stress  response. 


Dynamic  taxi: 

The  first  step  in  determining  stress  cycles  for  dynamic  taxi  is  to  select  one  of  25 
dynamic  taxi  cases  by  performing  a  best  match  of  taxi  segment  to  taxi  case 
parameters.  Then  from  the  chosen  taxi  cases  analytic  stress  response  (  A,  NJ 
and  the  taxi  segments  duration  (T)  generate  delta  stress  (Aa)  cycles  by 
application  of  Monte  Carlo  methods.  To  these  Aa  cycles  add  a  l  .Og  stress 
determined  by  application  of  stress  transfer  coefficients,  determined  by 
regression  analysis,  to  taxi  segment  parameters. 

In  order  to  utilize  recorded  data  to  determine  the  runway  roughness  enhanced 
methods  will  be  developed.  First,  an  effort  will  be  made  to  develop  a  method  to 
evaluate  C.G.  Nz  to  distinguish  between  prepared  and  semi-prepared  airfields 
(currently  all  airfields  are  assumed  to  be  prepared).  Following  this,  when  a 
sufficient  database  of  L/ESS  strain  gauge  data  is  available,  it  will  be  evaluated 
to  develop  a  C.G.  Nz  to  stress  relationship. 
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Stress  Spectra  Generation 


•  Landing  impact 

°  Use  the  highest  Nz  occurring  during  landing  impact  and  along  with 
landing  GW  analytically  convert  to  aircraft  sink  speed. 

°  Using  sink  speed  and  weight  (GW,  FW,  or  PW)  match  to  one  of  15 
landing  impact  conditions. 

°  Ratio  landing  condition  analytic  stress  sequence  to  match  segment 
sink  speed  and  weight. 

•  Take-off  rotation 

°  Assume  the  elevator  deflection  required  to  rotate  the  aircraft  is  the 
same  as  used  for  analysis. 

°  Apply  regression  transfer  coefficients  to  assumed  elevator  deflection 
and  flight  segment  parameters  to  obtain  stress  during  rotation. 


Landing  impact: 

Control  point  stress  sequences  due  to  landing  impact  are  generated  as  a 
function  of  weight  (gross  weight,  fuel  weight,  or  payload  weight)  and  aircraft 
sink  speed  using  the  following  procedure: 

1 .  The  highest  recorded  Nz  occurring  during  the  landing  impact  along  with 
the  landing  gross  weight  is  converted  to  aircraft  sink  speed. 

2.  One  of  1 5  landing  impact  conditions  is  chosen  by  performing  a  best  match 
of  sink  speed  and  weight.  The  weight  used  is  dependent  on  control  point 
location  (example:  for  wing  control  points,  fuel  weight  is  used). 

3.  The  chosen  landing  conditions  stress  sequence  is  ratioed  to  match  the 
segment  sink  speed  and  weight. 

Take-off  rotation: 

Control  point  stresses  that  occur  during  take-off  rotation  are  determined  as  a 
function  of  elevator  defection  and  associated  flight  segment  parameters.  The 
elevator  deflection  used  was  assumed  to  be  the  same  as  that  used  for  analysis. 
Regression  stress  transfer  coefficients  were  applied  to  the  assumed  elevator 
deflection  and  flight  segment  parameters  to  obtain  stress  during  rotation. 
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Control  Point  Selection 


Damage  Tracking 

•  Nine  control  points  throughout  the  airframe  were  selected  to  monitor  the 
structural  status  of  each  individual  aircraft. 

•  The  control  points  were  selected  on  each  major  component  to  track  the 
crack  growth  life  based  upon  individual  aircraft  usage. 

•  Stress  intensity  solutions  for  each  control  point  are  in  tabular  form  and 
damage  calculation  is  based  on  one  dimensional  (crack  length)  cycle- 
by-cycle  crack  growth  analysis  to  minimize  computing  time. 

•  Life  at  any  location  on  the  airplane  can  be  correlated  to  one  of  the 
control  points  based  upon  the  ratio  of  the  analytical  life  provided  the  two 
locations  have  the  same  loading  characteristics. 


The  continual  assessment  of  the  in-service  structural  integrity  of  the  individual 
airplanes  will  be  based  on  damage  calculated  for  each  of  the  control  points 
selected  throughout  the  airframe.  Control  points  were  selected  on  each  major 
component  to  monitor  the  change  of  each  damage  driving  parameter  (gust, 
maneuver,  cabin  pressure,  and  etc)  and  their  impact  on  the  life  of  the 
component  as  usage  varies.  Location  of  control  points  selected  on  each 
airframe  component  was  based  on  knowledge  gained  from  the  life  assessment 
effort  and  during  full-scale  engineering  development.  Damage  accumulation  is 
based  on  cycle-by-cycle  crack  growth  at  each  control  points.  The  stress 
intensity  solution  for  each  control  point  is  stored  in  tabular  form  and  crack 
growth  analyses  are  performed  in  the  length  direction  only  to  minimize  the 
computing  time.  Structural  life  at  any  location  on  the  airplane  can  be  derived 
from  one  of  the  existing  control  points  based  upon  the  ratio  of  the  analytical 
life  of  the  two  locations.  This  is  valid  provided  the  spectra  at  the  two  locations 
have  the  same  loading  characteristics  for  the  same  mission  event. 
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Control  Point  Selection 


•  AF1  -  Aft  Fus  Canted  Bhd  Cap 

•  AF2  -  Aft  Fus  Pressure  Bhd 

•  CF1  -  Long  #1  at  Wing  Rear  Spar 

•  HS1  -  H.S.  Front  Spar 

•  PI  -  Inb’d  Pylon  Stub 

•  VS1  -  V.S.  Front  Spar 

•  WF1  -  Wing/Fus  Trapezoidal  Pnl 

•  W1  -  Lwr  Wing  Skin  @  Xw  429 

•  W2  -  Wing  Upr  Skin  @  Xw165 


luoEie 

CPID  LOCATION /DESCRIPTION  DAMAGE  DRIVING 

PARMETERS 


AFl""'" . 

Aft  Fuselage,  Canted  Front  Spar,  Panel 
Close  to  Longeron  #3  (5)  Yf  =  1 582.90 

Lateral  Gust  & 
Maneuver  Loads 

AF2 

Aft  Fuselage,  Pressure  Bulkhead  @  (X46.9 
and  Z304  Aft  Cap)  Y,  =  1 797.55 

;  Pressure 

CF1 

Center  Fuselage;  Longeron  #1  (intercostal 
rib)  Wing  rear  Spar  @  Yf  851 .60 

Gust  &  Maneuver 
Loads 

HS1 

Horizontal  Stabilizer;  Front  spar  Upper  Cap 
@  Xhfs  =  126.87 

Gust  &  Maneuver 
Loads 

PI 

Inboard  Pylon-Stub,  Inboard  Upper  Cap  at 
■  Hole ‘A1 9’ @  Yin  =  263.1 4 

Gust  Loads 

VS1 

Vertical  Stabilizer;  Front  Spar  Skin  Flange 
^  @Zvie  =  613.00 

Lateral  Gust  & 
Maneuver  Loads 

WF1 

;  Wing/Fuselage,  Rear  Trapezoidal  Panel 
Upper  Inner  Doubler  at  Edge  @  Yf  = 

892.00  (Xw  =  116.0) 

Gust  &  Maneuver 
Loads 

W1 

Lower  Wing  Basic  Structure,  Skin  (§>  Xw  = 
429.00  and  Stringer  #50 

Gust  &  Maneuver 
Loads 

W2 

Upper  Wing  Basic  Structure,  Skin  @  Xw 
=1 65.00  and  Stringer  #20 _ 

Taxi  Loads 
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A  gap  in  the  flight  data  being  processed  may  be  caused  by  various  factors  such  as: 
malfunctioning  equipment,  exceeding  memory  capacity,  bad  data  rejected  during  the 
EDITING  process.  The  amount  of  crack  growth  for  the  gap  period  is  determined  by 
factoring  the  baseline  or  analytical  crack  growth  curve  based  on  the  design  usage  for 
each  respective  control  point  by  the  recent  usage  or  history  of  the  individual  airplane. 
The  procedure  used  is  summarized  as  follows: 

1 .  The  edited  flight  data  is  assessed  by  the  damage  calculation  program  and  the  gap 
in  recorded  usage  is  identified  in  terms  of  the  number  of  hours  unaccounted  for. 
This  is  .  The  total  flight  hours  from  initial  delivery  to  the  end  of  the  last 
recorded^ffight  is 

2.  From  the  appropriate  baseline  curve,  determine  the  time  required  for  the  crack  to 
reach  the  size  at  the  end  of  the  last  recorded  flight  (aj,  from  the  initial  size.  This 
is  (N'Jsi 

3.  The  equivalent  of  Tg^p  on  the  baseline,  (Tg^p)BL,  is  determined  as  follows: 

V-r 

sap)  gap 

4.  The  baseline  curve  is  entered  at  a„,  and  is  followed  for  a  time  interval  equal  to 
(T.ap^BL-  The  increment  of  growth  due  to  the  unlogged  hours  of  usage. 

5 .  The  crack  size  at  the  end  of  the  unlogged  period  is  given  by: 

+ 1  = 

The  procedure  presented  above  applies  to  all  control  points.  By  using  the  ratio  of 
(^,)bl  io  (^JiATP>  difference  between  the  baseline  usage  and  lATP  usage  is 
compensated  to  some  extent. 
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Future  Damage  Projection 


Projection  of  the  time  for  a  crack  to  reach  its  limiting  size  is  made  for  each  of 
the  nine  control  points  on  each  airplane  in  the  force.  The  procedure  is  based  on 
the  assumption  that  the  relationship  between  the  recorded  usage  and  the  design 
usage  crack  growth  curve,  in  terms  of  the  ratio  between  flight  hours  to  reach  the 
same  crack  size,  will  remain  constant  out  to  limiting  crack  size,  The 
procedure  to  be  used  to  project  the  time  to  reach,  from  a  given  crack  size, 
and  time,  is  as  follows: 

1 .  The  baseline  curve  is  used  to  determine  the  time  to  reach  from  the  initial 
crack  size  based  on  baseline  usage.  This  is  (NJ^p. 

2.  The  baseline  curve  is  used  to  determine  the  time  to  reach  the  limiting  crack 
size.  Op.  This  is 

3.  The  projected  time  to  reach  the  limiting  crack  size  is  found  from: 

)iATP  =  (N^  )bl  X 
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lAT  Data  Report 

•  Fleet  usage  leaders  (Airframe  time,  flights,  landings,  etc.). 

e  Usage  by  aircraft,  duty  base,  and  fleet. 

•  Time  spent  in  segment  (ascent,  cruise,  air-drop,  etc.). 

•  Mission  type  usage  distributions. 

•  Take-off  and  landing  gross  weight,  fuel  weight,  and  cargo  weight 
distributions. 

•  Air  drop  cargo  weight  distribution. 

•  Fleet  structural  status:  Fleet  leading  current  and  projected  damage 
tolerance  and  durability  life  at  the  control  points  . 

•  Individual  aircraft  structural  status;  Current  and  projected  damage 
tolerance  and  durability  life  at  the  control  points  for  each  aircraft. 

•  Data  capture  rate 

•  Duty  base  and  horizontal  tail  assignment  dates. 

Note;  All  data  compiled  for  the  current  report  period  and  cumulative. 


The  output  of  I  ATP  is  the  I  AT  Data  Report  which  typically  is  generated  on  a 
semi-annual  basis.  This  report  contains  two  types  of  information:  one  being 
usage  statistics  and  the  other  being  aircraft  structural  status.  Usage  statistics 
give  information  such  as  airframe  time,  nximber  of  flights  and  landings,  etc. 
This  information  is  given  by  aircraft,  duty  base,  fleet,  and  fleet  leaders  (aircraft 
with  most  flights,  landings,  etc.).  Information  is  also  given  in  terms  of  time 
spent  in  the  different  segment  types  and  a  distribution  of  flights  over  the  35 
design  mission  types.  In  addition,  take-off  and  landing  gross  weight,  fuel 
weight,  and  cargo  weight  distributions  as  well  as  air-drop  cargo  weight 
distribution  are  given. 

Aircraft  structural  status  is  given  as  current  and  projected  damage  tolerance 
and  durability  life  at  the  control  points.  This  information  is  given  for  each 
aircraft  and  also  summarized  by  the  fleet  leading  aircraft  (aircraft  with  largest 
damage). 

In  order  to  disposition  the  success  of  collecting  and  processing  SFDR  data 
capture  rates  are  given  in  terms  of  flight  hours  processed  over  flight  hours 
flown.  Also  given  are  a  history  for  each  aircraft  of  duty  base  and  horizontal 
stabilizer  assignment  dates. 

The  above  information  is  compiled  for  the  current  report  period  and 
cumulative. 
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Computer  System  Characteristics 


•  Platform:  IBM  mainframe  computer 

•  Source  code  written  in  FORTRAN 

•  System  runs  interactively  or  as  a  batch  process  depending  on  the 
operations  being  performed. 

•  Employs  a  DB2  relational  database  for  ease  of  data  processing  and 
managing. 

•  Incorporates  utilities  to  permanently  backup  processed  flight  data  to  tape 
in  order  to  relieve  on-line  storage  requirements. 

•  User  interface  in  the  form  of  pull  down  menus  and  panels  for  ease  of 
operation  and  user  data  entry. 


The  lATP  system  resides  on  an  IBM  main  frame  computer  and  is  composed  of 
12  separate  FORTRAN  programs  that  perform  the  various  data  processing 
functions.  This  system  runs  interactively  for  those  functions  requiring 
extensive  user  interaction,  or  as  a  batch  process  for  those  functions  requiring  a 
large  amount  of  data  processing  time.  The  system  employs  a  DB2  relational 
database  for  storing  and  retrieving  processed  flight  data  and  user  input  data.  A 
relational  database  allows  for  efficient  data  storage  and  ease  of  data  retrieval 
and  manipulation  without  extensive  coding.  The  system  includes  a  utility  to 
permanently  backup  processed  flight  data  to  tape  in  order  to  relieve  on-line 
storage  requirements.  The  user  interface  is  in  the  form  of  pull  down  menus  and 
panels.  This  allows  for  ease  of  system  operations  and  user  data  entry  without 
extensive  knowledge  of  mainframe  computer  operations. 
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Heather  P.  Roland  TSgt  Alfred  G.  Tans 

Lockheed  Martin  Altus  AFB 

Marietta,  Georgia  Altus,  Oklahoma 


Background 
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•  Safety  concerns  expressed  due  to  lateness 
of  AFTO  451  usage  form  submittals. 


Background 
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•  Weekly  tracking  runs  required  to  insure 
safety. 


Purpose 


s  «  ^ 

03  T3 

p  3  S 

u  w 

S  o 

OS)  ^  , 

2  3  2f3 

A  P  C5  '*5 
mh  ;5  </:i  P 

1  £ 
"S  o 

^  S  =■!! 

O  <p 

;n 

1/^ 

O  “ 

b  3 

CJ  X} 

‘S  a 

2-3 

45  a> 
o  >► 
a>  ^ 

S 

^  X 

-fi  ?? 

H  ^ 


*a 

5?  c« 

M 

2 


CQ 


a> 


ejD*rt 

^  Qj 

a 

u 

Cm 

o 


u 

u 

a 

QJ 

a 


c/5  _5i 

c/5 

QJ  .5 

U  > 

o  o 
u  u 

Q^  CU 


M  0^ 

^  c« 

pi2 

'TS  ^ 

3  « 
.1'^ 

Ph  ^ 


a> 

o 


s 

o 

o 
o 

c«  ^ 

S 

^  O 

.2  o 

.P 

H 


p 

'P 

dJD 


1098 


U.S.  Air  Force  C-141B  Starlifter 
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standard  Missions 
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AD2C  SKE:  Three  Lifts 

AD3C  SKE/VFR:  Three  Lifts 


Advantages  of  the  New 
Altus  Form  1 
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entered  daily  into  the  Electronic 
451  program  for  quick  turnaroun 


The  Electronic  AFTO  451 

Main  Menu 


Data  Entry  Screen 
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Electronic  Data  Transfer 


Advantages 
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•  Flight  data  accounted  for  in  tracking 
within  a  week  of  actual  flight  date. 


Additional  Benefits 
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FAA  MSR/LSR  Flight  Inspection  Fleet 
Aircraft  Structural  Integrity  Program 


Mr.  James  M.  Marks  Jr. 

Mr.  James  J.  Abel 
Raytheon  E-Systems 
Greenville  Division 

P.O.  Box  6056  Greenville,  Texas  75403-6056 


ABSTRACT 

Raytheon  E-Systems  has  developed  and  implemented  a  comprehensive 
Aircraft  Structural  Integrity  Program  (ASIP)  that  aims  to  ensure  the  long-term 
structural  integrity  and  continued  airworthiness  of  a  broad  portion  of  the  Federal 
Aviation  Administration  (FAA)  Flight  Inspection  Fleet  of  aircraft.  Specifically,  six 
Leaijet  Model  60  and  three  Canadair  Challenger  CL601-3R  aircraft  are  modified 
by  Raytheon  E-Systems  to  perform  the  Flight  Inspection  role  for  the  FAA.  Flight 
Inspection  mission  profiles  are  characterized  by  routine  operation  below  2000  feet 
MSL  (mean  sea  level).  This  low  level,  structurally  harsh  environment  contrasts 
significantly  with  typical  corporate  jet  operational  environments.  The  Flight 
Inspection  usage  has  a  notable  impact  on  the  loading  spectra  experienced  by  the 
aircraft  structure  due  to  the  routine  operation  in  the  low-level  gust  environment, 
and,  as  a  result,  an  appreciable  impact  on  the  service  life  of  the  aircraft. 

The  ASIP  is  based  on  a  tailored  implementation  of  the  requirements  of  the 
U.S.  Ar  Force  MIL-STD-1530  approach.  The  ASIP  task  elements  consist  of  a 
Durability  and  Damage  Tolerance  Assessment  (DADTA),  an  Individual  Arcraft 
Tracking  Program  (lATP),  a  Loads/Environment  Spectra  Survey  (L/ESS),  and  a 
Fleet  Structural  Maintenance  Plan  (FSMP).  The  interrelationship,  implementation 
and  functional  responsibilities  of  each  ASIP  task  element  are  identified  and 
distinguished  in  an  ASIP  Master  Plan.  This  paper  focuses  on  the  development  and 
integration  of  the  lATP  as  well  as  the  relationship  of  the  lATP  with  the  other  ASIP 
task  elements. 
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INTRODUCTION 


The  objective  of  the  I  ATP  for  the  FAA  Flight  Inspection  Fleet  is  to  provide  a 
method  for  collecting,  processing,  and  analyzing  data  from  various  sensors 
installed  onboard  each  aircraft  in  order  to  predict  crack  growth  at  potentially 
critical  areas  on  each  airframe.  The  Durability  and  Damage  Tolerance  Analysis 
(DADTA)  for  both  aircraft  types  [1]  [2]  identifies  these  critical  areas  as  well  as  the 
corresponding  predicted  crack  growth  expected  in  service,  predicated  upon  an 
assumed  initial  baseline  loading  spectrum.  Crack  growth  calculations  using 
measured  flight  data  are  then  compared  to  the  baseline  crack  growth  calculations 
contained  in  the  DADTA.  The  results  are  used  to  project  the  scheduling  of 
inspection  and/or  maintenance  activities  identified  in  the  Fleet  Structural 
Maintenance  Plan  (FSMP)  for  both  the  Leaijet  Model  60  and  the  Canadair 
Challenger  CL601-3R.  Prudent  fleet  management  decisions  may  then  be  formed 
based  upon  actual  in-service  usage,  some  of  which  may  include:  accelerating  or 
delaying  inspections  and/or  maintenance  actions,  re-assigning  aircraft  from  one 
base  to  another  in  order  to  evenly  distribute  mild  and  harsh  usage  environments, 
and  re-evaluation  of  the  DADTA  based  on  the  measured  flight  load  spectrum  after 
a  significant  amount  of  flight  data  is  collected. 

CURRENT  AIRCRAFT  IN  THE  lATP 

The  lATP  is  designed  to  track  all  Learjet  Model  60  and  Canadair  Challenger 
CL601-3R  aircraft  in  the  Flight  Inspection  Fleet.  The  Learjet  Model  60  performs 
the  Medium  Size/  Medium  Range  (MSR)  flight  inspection  role  for  the  FAA.  The 
Canadair  Challenger  CL601-3R  performs  the  Large  Size/  Long  Range  (LSR)  flight 
inspection  role  for  the  FAA.  These  aircraft  differ  in  size,  manufacturer,  fracture 
analysis  methodology,  flight  recorder  parameters,  and  original  certification  basis. 

A  depiction  of  each  aircraft  is  shown  in  Figures  1  and  2,  respectively. 

The  Leaijet  Model  60  aircraft  has  a  maximum  gross  takeoff  weight  of  23,500 
pounds  with  a  maximum  operating  altitude  of  51,000  feet.  The  Learjet  Model  60 
is  certified  to  Amendment  23  of  FAR  Part  25.571  [3]  using  fatigue  and  static  fail¬ 
safe  methodologies.  Leaijet,  however,  performed  a  fracture  mechanics  based 
damage  tolerance  analysis  of  the  modified  Model  60  MSR  Flight  Inspection 
aircraft  using  assumed  baseline  load  spectra  for  the  planned  flight  inspection 
operating  environment. 

The  Canadair  Challenger  CL601-3R  aircraft  has  a  maximum  gross  takeoff 
weight  of  45,100  pounds  with  a  maximum  operating  altitude  of  41,000  feet.  The 
CL601-3R  is  certified  to  Amendment  45  of  FAR  25.571  [3]  as  a  “damage 
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tolerant”  aircraft  using  fracture  mechanics  based  damage  tolerance  analysis 
methodologies.  In  addition,  Canadair  re-evaluated  the  structural  integrity  of  the 
CL601-3R  to  account  for  the  unique  load  spectra  encountered  in  the  flight 
inspection  environment. 

FLIGHT  DATA  RECORDER 

The  Flight  Data  Recorder  (FDR)  system  is  installed  on  each  aircraft  in  the 
fleet.  The  FDR  system  is  comprised  of  eight  major  line  replaceable  units  (LRUs):  a 
Signal  Acquisition  Unit  (SAU),  a  Data  Transfer  Interface  Unit  (DTIU),  a  Data 
Transfer  Module  (DTM),  a  Crash  Survivable  Memoiy  Unit  (CSMU),  a  Flight  Data 
Panel  (FDP),  an  Engine  Signal  Data  Converter  (ESDC,  LSR  aircraft  only),  a  Flight 
Data  Recorder  Interface  Unit  (FDRIU),  and  a  triaxial  accelerometer.  In  addition, 
eight  strain  sensors  are  permanently  bonded  to  each  MSR  airframe  and  ten  strain 
sensors  are  permanently  bonded  to  each  LSR  airframe.  The  FDR  system  diagram 
is  shown  in  Figure  3.  The  strain  sensor  locations  for  both  the  MSR  and  LSR 
aircraft  are  shown  in  Figures  4  and  5,  respectively. 

The  SAU,  DTIU,  DTM,  CSMU,  and  FDP  are  supplied  by  Smiths  Industries. 
The  ESDC  and  FDRIU  for  the  LSR  aircraft  are  supplied  by  ICE  Corporation.  The 
FDRIU  for  the  MSR  aircraft  is  supplied  by  Leaijet.  The  accelerometer  is  supplied 
by  Magnetek  Transducer  Products.  The  strain  sensors  are  supplied  by  Columbia 
Research  Labs,  Inc. 

All  aircraft  signals  monitored  by  the  FDR  system  interface  with  the  SAU, 
which  processes  the  data  using  its  internal  Operational  Flight  Program  (OFP)  and 
then  stores  it  into  the  two  separate  non-volatile  memory  devices:  the  CSMU  and 
the  DTM.  The  CSMU  is  permanently  mounted  in  each  aircraft’s  tail  section  and  is 
to  be  retrieved  in  the  event  of  a  mishap.  Flight  data  required  by  FAR  Part  135 
Appendix  B  are  stored  in  the  CSMU.  The  DTM,  which  is  inserted  into  the  face  of 
the  DTIU  (much  like  a  cassette  cartridge)  is  used  to  periodically  transfer  the  stored 
engine  and  structural  data  from  the  aircraft  to  a  ground-based  computer  system  for 
processing  and  analysis. 

The  SAU  also  receives  user  input  data  from  the  FDP.  These  data  are  entered 
by  the  crew  during  preflight  and  include  such  parameters  as  fuel  weight,  gross 
weight,  center  of  gravity,  date,  mission  base,  etc.  The  FDP  also  provides  a  built-in 
test  (BIT)  button  and  fault  annunciators,  as  well  as  memory  capacity  indications 
that  signal  DTM  memory  at  80%  (or  greater)  and  100%  full. 
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The  DTIU  provides  a  receptacle  for  the  DTM  and  transfers  the  data  to  the 
DTM  for  download.  It  is  also  used  for  uploading  the  OFF  and  aircraft  specific 
configuration  data  from  the  ground-based  sources.  The  DTIU  also  contains  a 
display  and  operator  interface  used  for  upload/download  commands,  as  well  as  for 
displaying  system  status  messages. 

The  ESDC  unit  (LSR  aircraft  only)  converts  existing  analog  engine  signals 
ITT,  Nl,  and  N2  (from  both  engines)  into  ARINC  429  high  speed  digital  data 
format.  ITT,  Nl,  and  N2  engine  signals  are  digitized  to  ARINC  429  digital  data 
format  because  these  three  parameters  require  the  highest  accuracy  and  are  the 
most  significant  contributors  in  the  engine  trend  analysis.  The  Engine  Signal  Data 
Converter  Unit  digitizes,  combines,  and  transmits  the  ITT,  Nl,  and  N2  data  to  the 
Flight  Data  Recorder  system  on  two  separate  ARINC  429  data  bus  channels  (left 
engine  and  right  engine  data  channels).  The  ESDC  incorporates  BIT  features  that 
allow  fault  isolation  to  the  LRU  level. 

The  FDRTU  is  essentially  a  relay  box  that  converts  discrete  signals  from 
various  parts  of  the  aircraft  into  a  format  that  is  acceptable  to  the  SAU.  The  S  AU 
accepts  mostly  open/28VDC  type  discrete  signals  and  only  a  few  open/ground  type 
discrete  signals;  whereas,  the  aircraft  provides  mainly  open/ground  type  discrete 
signals.  The  primary  function  of  the  FDRIU,  therefore,  is  to  convert  a  number  of 
discrete  signals  from  open/ground  format  to  open/28VDC  format.  The  FDRIU 
also  provides  diode  isolation  between  the  monitored  aircraft  signals  and  the  SAU. 
The  FDRIU  incorporates  BIT  features  that  allow  fault  isolation  to  the  LRU  level. 

The  accelerometer  provides  three-axis  acceleration  data  to  the  SAU  and  is 
used  in  determining  airframe  structural  loads  due  to  maneuvers,  gusts,  and  ground 
events.  The  strain  sensors  are  permanently  installed  (bonded)  at  critical  locations 
throughout  the  airframe  for  measuring  and  recording  structural  strains.  The  strain 
sensors  provide  their  signal  data  to  the  SAU.  These  data  are  used  to  track  and 
maintain  a  structural  loading  history  of  the  airframe. 

Stored  engine  and  structural  data  are  periodically  downloaded  from  the 
aircraft  to  a  ground-based  computer  system  via  the  DTM  for  processing  and 
analysis.  The  computer  system  is  called  the  Ground  Replay  and  Display  Unit 
(GRDU)  and  is  supplied  by  Smiths  Industries.  It  consists  of  an  IBM  type  personal 
computer  (PC)  with  an  internal  DTM  interface  card,  and  an  external  Data  Transfer 
Module  Receptacle  (DTMR).  The  GRDU  contains  several  software  programs 
used  for  engine  trending  and  structural  analysis,  including  the  Smiths  Industries 
Operational  Ground  Program  (OGP)  and  various  analysis  programs.  Engine 
performance  trend  analysis  for  the  LSR  aircraft  is  accomplished  by  processing  the 
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downloaded  engine  data  with  a  General  Electric  supplied  software  program 
named  CF34  RJ  TREND  VERSION  4.0  (RJ  stands  for  Regional  Jet;  this  program 
is  normally  used  for  trend  analysis  of  the  CF34  type  variant  engines  on  the 
Regional  Jet  Challenger  aircraft,  but  is  also  applicable  to  the  CF34  engines  on  the 
CL601-3R).  Engine  performance  trend  analysis  for  the  MSR  aircraft  is 
accomplished  by  processing  the  downloaded  engine  data  with  a  Pratt  &  Whitney 
supplied  software  program  named  ECTM IV.  Structural  Integrity  data  analysis  is 
performed  on  a  Hewlett  Packard  workstation  using  software  supplied  by  Raytheon 
E-Systems  as  part  of  the  Individual  Aircraft  Tracking  Program  (lATP). 

The  GRDU  can  also  be  used  as  a  Direct  Parameter  Display  (DPD)  by 
connecting  it  to  the  FDR  system  via  an  aircraft  mounted  test  connector  and  a  test 
cable.  The  test  connector  is  called  the  Ground  Replay  Equipment  (GRE) 
connector,  and  it  is  installed  inside  the  auxiliary  equipment  rack.  The  test  cable  is 
called  the  Ground  Replay  Equipment  (GRE)  cable.  The  DPD  allows  individual 
parameters  to  be  monitored  real  time,  but  the  FDR  system  cannot  record  data 
when  in  DPD  mode.  In  addition,  a  laptop  computer  DPD  is  available  to  monitor 
parameters  real  time  in  lieu  of  using  the  GRDU. 

TRACKING  ANALYSIS  METHODOLOGY 

Raw  flight  recorder  data  are  downloaded  periodically  (approximately  every 
two  weeks  or  thirty  flight  hours)  from  each  aircraft  in  the  Flight  Inspection  Fleet. 
The  data  are  transmitted  via  modem  to  the  FAA  Flight  Inspection  Office  in 
Oklahoma  City,  OK,  from  the  various  field  bases  located  throughout  the  United 
States  and  abroad.  The  data  are  then  forwarded  to  Raytheon  E-Systems  in 
Greenville,  TX,  via  the  internet  for  subsequent  processing  in  the  lATP. 

The  Raytheon  E-Systems  lATP  software  reads,  stores,  and  processes  raw 
flight  recorder  data  that  have  been  downloaded  and  processed  on  the  GRDU  using 
the  Smiths  Industries  supplied  Cartridge  Ground  Program  (CGP)  software.  The 
processed  decompressed  data  files  (ddf)  are  transferred  from  the  GRDU  to  a 
Hewlett  Packard  (HP)  workstation.  The  Raytheon  E-Systems  lATP  software 
resides  on  this  HP  workstation. 

The  lATP  software  reads  and  stores  each  ddf  in  a  database.  These  data  form 
the  Loads/Environment  Spectra  Survey  (L/ESS).  The  L/ESS  database  contains 
the  structural  loading  history  for  all  critical  areas  on  each  aircraft  in  the  fleet  as 
recorded  by  each  aircraft’s  individual  FDR.  The  I  ATP  software  utilizes  these  data 
to  perform  crack  growth  analyses  at  these  critical  locations. 
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The  lATP  reads  and  interprets  the  FDR  BIT  file.  The  BIT  file  indicates  the 
fault  status  of  various  FDR  parameters.  The  results  are  displayed  on  the  screen  and 
are  also  stored  in  a  separate  file  that  may  be  accessed  later  by  the  user.  Those 
parameters  critical  for  crack  growth  algorithms  are  highlighted.  Also,  all  FDP 
entered  parameters,  strain  sensors  and  accelerations  are  verified  for  plausibility  by 
determining  if  the  data  exceed  certain  predetermined  maximum  and  minimum 
values.  Raytheon  E-Systems  furnishes  a  monthly  report  to  the  FAA  customer 
detailing  any  BIT  faults  detected  by  the  lATP.  This  information  is  used  to 
investigate  and  correct  potential  hardware  or  software  problems  associated  with 
the  FDR. 

The  lATP  software  utilizes  an  object-oriented  architecture  with  a  fiilly 
integrated  graphical  user  interface.  One  significant  advantage  realized  from  the 
object-oriented  software  architecture  is  the  ability  to  add  new  aircraft  types, 
beyond  the  Learjet  Model  60  and  Canadair  Challenger  CL601-3R,  with  no 
modification  to  the  essential  software  code.  The  graphical  user  interface  enables 
the  software  to  be  user-friendly  and  dynamic. 

lATP  FATIGUE  CRACK  GROWTH 

The  lATP  fatigue  crack  growth  methodologies  are  based  on  the  Leaijet 
Model  60  MSR  Service  Life  Analysis  [1]  and  the  Canadair  Challenger  CL601-3R 
LSR  DADTA  [2].  The  theoretical  growth  of  cracks  at  each  control  point,  from  an 
assumed  initial  rogue  flaw  in  the  aircraft  structure,  is  determined  utilizing  the 
applied  stress  spectrum  as  collected  by  each  flight  data  recorder.  A  cycle-by-cycle 
analysis  approach  is  utilized  such  that  a  full  accounting  of  sequenced  peak  and 
valley  stress  data  are  incorporated.  The  lATP  utilizes  a  standard  cycle  counting 
method  known  as  Rain  Flow  Counting.  This  stress  spectrum  cycle  counting 
method  ensures  that  the  largest  delta  stress  cycles  (Omax  -  c^min)  are  not  neglected  in 
the  analysis.  Control  point  locations  for  both  the  MSR  and  LSR  aircraft  are  shown 
in  Figures  6  and  7,  respectively. 

Stress  spectra  are  determined  in  two  distinct  manners.  The  primary  method  is 
to  use  data  from  the  strain  sensors  installed  onboard  each  aircraft.  The  secondary 
alternative  is  to  calculate  the  stress  spectra  based  upon  other  parameters  collected 
by  the  FDR  such  as  acceleration  (g’s),  aircraft  weight,  airspeed  and  altitude.  The 
lATP  will  use  the  secondary  method  only  when  the  necessary  strain  sensor  data 
are  not  available  on  a  control  point-by-control  point,  decompressed  data  file-by- 
decompressed  data  file  basis.  Transfer  functions  are  established  that  define  specific 
relationships  between  strain  sensor  values  and/or  various  combinations  of  other 
parameters  collected  by  the  FDR  in  order  to  determine  the  appropriate  values  of 
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stress  at  each  control  point  location.  In  addition,  if  the  FDR  data  are  completely  or 
partially  unusable,  a  gap-filling  technique  is  employed  by  the  lATP.  This  is 
accomplished  by  copying  portions  of  decompressed  data  files  from  previously 
flown  missions  in  an  attempt  to  create  data  that  are  representative  in  terms  of  flight 
hours,  landings,  and  pressure  cycles. 

The  fatigue  crack  growth  analysis  involves  the  modeling  of  various  structural 
control  points  on  the  Leaijet  Model  60  and  Canadair  Challenger  CL601-3R.  Each 
control  point  model  consists  of  a  crack  growth  system  with  one-dimensional  and 
two-dimensional  crack  growth  combined  to  achieve  a  multi-phase  crack  growth 
capability.  Each  crack  growth  model  includes  the  relevant  material  properties, 
stress  intensity  solution,  and  crack  growth  retardation  model.  Crack  growth 
involves  an  incremental  growth  (Aa)  per  stress  cycle  (AN).  Material  properties  are 
included  in  the  form  of  da/dN  data.  These  data  are  a  function  of  applied  stress 
spectrum  (amax  and  Omin)  and  environment.  Control  point  geometry  effects  are 
included  through  the  use  of  a  factor,  P,  such  that  the  stress  intensity,  K,  is  given 
by: 


K  =  (Eqn.  1) 

Crack  growth  is  retarded  when  an  enlarged  plastic  zone  is  developed  at  the 
crack  tip  due  to  an  overload.  The  crack  growth  retardation  model  incorporated  in 
the  lATP  is  the  industry-accepted  Generalized  Willenborg  Retardation  Model. 

INSPECTION  INTERVALS 

Inspection  intervals  for  each  control  point  are  based  upon  a  linear  function  of 
time  with  respect  to  structural  failure  of  the  control  point.  Initial  as  well  as 
recurring  inspection  intervals  are  determined  by  each  aircraft  manufacturer.  The 
lATP  includes  a  comparison  of  calculated  crack  growth  based  on  actual  usage 
versus  predicted  crack  growth  in  terms  of  crack  length  versus  flight  time.  Criteria 
are  implemented  in  the  lATP  to  adjust  inspection  intervals  when  crack  growth 
based  on  actual  usage  deviates  from  predicted  crack  growth. 

Figure  8  depicts  the  initial  inspection  interval  adjustment  criterion.  Let  tip  be 
the  time  corresponding  to  the  predicted  initial  inspection.  Let  ta  be  the  time 
corresponding  to  the  current  crack  length  based  on  actual  usage.  Let  tp  be  the  time 
at  which  the  predicted  crack  growth  curve  would  yield  the  same  crack  length 
corresponding  to  the  crack  growth  curve  based  on  actual  usage.  Therefore,  the 
adjusted  initial  inspection  time  (tia)  as  a  function  of  predicted  initial  inspection  time 
is: 
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(Eqn.  2) 


tia  =  (ta/tp)  *  tip 

The  delta  time  to  next  inspection,  Ati,  is  therefore  defined  as: 

Ati  =  tia  -  ta  (Eqn.  3) 

Figure  9  depicts  the  recurring  inspection  interval  adjustment  criterion.  This 
criterion  is  identical  to  the  initial  inspection  adjustment  criterion  except  that  the 
times  ta,  tp,  tia,  and  hp  are  referred  to  the  time  of  last  inspection  rather  than  time 
zero. 


When  actual  crack  growth  deviates  sufficiently,  a  significant  change  in  usage 
is  reported  by  the  lATP.  This  occurs  when  there  is  a  10%  reduction  in  the 
predicted  inspection  time  interval  (tia/tip  <  0.90). 

The  lATP  also  has  the  ability  to  adjust  the  actual  usage  crack  growth  curve 
to  account  for  inspection  and  maintenance  results.  This  is  depicted  in  Figure  9. 
The  crack  length  is  adjusted  to  the  maximum  non-detected  length  if  no  crack  is 
found  during  inspection.  Theoretically,  the  crack  length  could  also  be  set  to  the 
length  discovered  during  inspection;  however,  in  practice,  all  cracks  found  during 
inspection  will  likely  be  repaired. 

INSPECTION/MAINTENANCE  SCHEDULING 

Semi-annually  the  lATP  data  are  compiled  and  reported  to  the  FAA 
customer  along  with  appropriate  discussion  and  recommendations.  Each  semi¬ 
annual  report  relates  the  structural  integrity  status  of  the  Flight  Inspection  Fleet. 
For  example,  the  following  information,  at  a  minimum,  is  included  in  each  semi¬ 
annual  report:  crack  growth  curve  and  stress  exceedance  curve  for  each  control 
point  for  each  aircraft  for  the  current  reporting  period  as  well  as  accumulated 
totals,  remaining  hours  to  inspection  or  maintenance  for  each  control  point  for 
each  aircraft,  a  listing  of  the  inspections  and  maintenance  actions  required  during 
the  following  two  calendar  years  for  each  aircraft;  flight  hours,  landings,  and 
pressure  cycles  grouped  by  individual  aircraft:  and  by  aircraft  type  (MSR  or  LSR) 
for  the  current  reporting  period  as  well  as  accumulated  totals,  and  a  stress 
exceedance  curve  for  each  control  point  grouped  by  aircraft  type  (MSR  or  LSR). 

Detailed  inspections  and  maintenance  procedures  are  included  in  the  FSMP 
or  maintenance  manuals  for  each  aircraft  type.  The  inspection  intervals  are  based 
upon  the  DADTA  performed  by  each  aircraft  manufacturer  using  assumed  initial 
baseline  load  spectra.  The  lATP  semi-annual  reports  will  convey  any  applicable 
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recommendations  relative  to  inspection  interval  adjustment  due  to  actual  usage. 
This  is  of  particular  importance  to  safety  in  that  each  aircraft  may  experience  a 
more  severe  loading  in  service  than  originally  anticipated. 

A  link  among  the  FAA  customer,  the  FAA  maintenance  infrastructure,  and 
Raytheon  E-Systems  is  vital  in  order  to  ensure  that  inspection/maintenance  results 
are  incorporated  into  the  lATP.  This  link  ensures  that  inspection/maintenance 
information  flows  freely  among  the  three  entities.  Each  entity  performs  a  function 
that  affects  the  overall  execution  of  the  ASEP  for  the  Flight  Inspection  Fleet. 

Inspection/maintenance  actions  are  dependent  upon  the  projected  usage  of 
the  fleet  considering  design  mission  profiles  and  the  stress  spectra  that  result  from 
these  missions.  These  baseline  stress  histories  will  either  be  validated  by  the  lATP 
or  indicate  that  an  update  is  required.  If  the  latter  is  found  to  be  true,  the  DADTA 
will  be  re-evaluated  with  more  accurate  stress  spectra  from  the  L/ESS.  This  may 
then  result  in  changes  to  the  inspection  intervals  found  in  the  FSMP  for  each 
control  point. 

CONCLUSION 

As  the  FAA  Flight  Inspection  MSR  and  LSR  aircraft  accumulate  flight 
hours  in  a  structurally  harsh  low-level  environment,  the  lATP  system  will  provide 
detailed  inspection  and  maintenance  recommendations  to  the  FAA  operator. 
Concurrently,  the  L/ESS  data  contained  in  the  lATP  will  further  define  the  loading 
spectrum  unique  to  the  operation  of  the  MSR  and  LSR  aircraft  in  the  flight 
inspection  environment.  The  careful  application  of  results  from  the  lATP  will  no 
doubt  increase  aircraft  safety  and  reduce  overall  maintenance  costs. 
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DEFINITIONS,  ACRONYMS,  ABBREVIATIONS 


ARJNC  Aeronautical  Radio  Incorporated 

ASIP  Aircraft  Structural  Integrity  Program 

BIT  Built-in  Test 

CGP  Cartridge  Ground  Program 

CSMU  Crash  Survivable  Memory  Unit 

DADTA  Durability  and  Damage  Tolerance  Analysis 

DDF  Decompressed  Data  File 

DPD  Direct  Parameter  Display 

DTIU  Data  Transfer  Interface  Unit 

DTM  Data  Transfer  Module 

DTMR  Data  Transfer  Module  Receptacle 

ESDC  Engine  Signal  Data  Converter 

FAA  Federal  Aviation  Administration 

FDP  Flight  Data  Panel 

FDR  Flight  Data  Recorder 

FDRIU  Flight  Data  Recorder  Interface  Unit 

FSMP  Fleet  Structural  Maintenance  Plan 

GRDU  Ground  Replay  and  Display  Unit 

GRE  Ground  Replay  Equipment 

HP  Hewlett  Packard 

lATP  Individual  Aircraft  Tracking  Program 

ITT  Interturbine  Temperature 

L/ESS  Loads/  Environment  Spectra  Survey 

LSR  Large  Size/  Long  Range  (Canadair  CL60 1 -3R) 

LRU  Line  Replaceable  Unit 

MSL  Mean  Sea  Level 

MSR  Medium  Size/  Medium  Range  (Lear  Model  60) 

N1  Engine  Fan  Speed 

N2  Engine  Core  Speed 

OFP  Operational  Flight  Program 

OGP  Operational  Ground  Program 

PC  Personal  Computer 

SAG  Structural  Advisory  Group 

SAU  Signal  Acquisition  Unit 
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Figure  1 

Learjet  Model  60  MSR  Flight  Inspection  Aircraft 
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Figure  2 

Canadair  Challenger  CL601-3R  LSR  Flight  Inspection  Aircraft 
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Figure  3 

Flight  Data  Recorder  System  Diagram 
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Figure  4 

FAA  Flight  Inspection  (MSR)  Lear  jet  Model  60 
Flight  Data  Recorder  Strain  Sensor  Locations 
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CP  #8:  Fwd  Fuselage  NLG  Wheel  Box  FS  235  BL  5.7 


Figure  7 

FAA  Flight  Inspection  (LSR)  Challenger  CL601-3R 
DADTA  Control  Point  Locations 


Recurring  Inspection  Adjustment  Criterion 
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